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Performance Evaluation of ERS-1 Synthetic Aperture Radar Fast Delivery 
Processor 


43070007a Kyoto Selection from INTERNATIONAL ARCHIVES OF PHOTOGRAMMETRY 
AND REMOTE SENSING in English Vol 27 1988 Part 88 pp 1-7 


{Article by J.P. Guignard, European Space Research and Technology Centre 
(ESTEC), Keplerlaan 1,2201 AZ Noordwijk, The Netherlands; P. Hasan, 
MacDonald Dettwiler, 3751 Shell Road, Richmond 8.C., Canada V6X 229, 


Commission No. II] 


[Text] Abstract 


The Factory Acceptance Test (FAT) of ESA’s ERS-1 Synthetic Aperture 
Radar Fast Delivery Processor (SAR FDP) took place recently at MacDonald 
Dettwiler in Vancouver, Canada. The process will eventually be 
installed in the principal ERS-1 ground station in Kiruna, Sweden. 


This paper describes the comprehensive set of tests performed during the 
FAT. It focusses on the tests and measurements of: 


--image quality parameters (3-dB resolution, side lobe level, etc.); 


--throughput for both image and wave mode preduct generation; 


--operational aspects of the processor demonstrating its ability to run 
in an operational environment. 


The SAR FDP shows a significant improvement in throughput compared to 
previous SEASAT processors, while preserving the integrity of the image 


quality. 
1.0 Introduction 


The development of the ESA ERS-1 satellite offers an opportunity to 
develop a fast processor for Synthetic Aperture Radar (SAR) data. 








The key requirement is to minimize the delay between data acquisition 
and delivery of processed images to users. A maximum delay of three 
hours is allowed for both the image and wave modes. In the image mode 
the satellite acquires SAR data for 100 x 100 km earth surface images. 
In the wave mode SAR data is acquired for 5 x 5 km images of the ocean 
surface at along-track intervals of 200 or 300 km. The SAR FDP 
processes wave mode data to an intermediate image, and then to a final 
inage spectrum product. Sea state information (wave spectrum product) 
will be derived from the image spectrum product. The image mode 
operates only while the satellite is visible from a ground station, 
while the wave mode can operate around the orbit, since data can be 
recorded on-board and later downlinked to the ground station (see [1]). 
















The features of a SAR processor capable of meeting such a requirement 
while preserving the quality of the product have been defined and an ESA 
contract placed with MacDonald Dettwiler for the processor’s 

development. Development based on the design presented in [2] is now 
complete, and this paper describes the processor’s actual performance. 












2.0 Methodology 






The SAR FDP was tested against all detailed requirements stated in the 
SAR FDP requirements specifications. A verification cross reference 
matrix was used to map each requirement to the test number in the 
acceptance test procedures document. 











A total of 16 main tests, each of them consisting of various subtests, 
were developed. These tests were designed to demonstrate that the SAR 
FDP meets its image quality requirements as well as its operational 

requirements. 

















The image quality tests were done by processing simulated point targets 
as well as pre-processed SEASAT data and then analysing the processed 
images. The operational aspects of the processor were tested by > 
scheduling the processor to generate all possible products from image 
mode and wave mode data. The processor was also tested against possible 
errors in commands from the control computer, or in the input raw data. 


3.0 Test Data Specifications 
3.1 Simulated Point Targets 


Point target data is generated using a specially developed simulator 
running on the ST-100 array processor (developed and manufactured by 
STAR Technologies Inc.). Both image mode and wave mode data are 
Simulated using ERS-1 parameters and copied to High Density Digital Tape 
(HDDT), which conforms to the ERS-1 downlink format. 


Parameters are chosen to simulate extreme conditions such as maximum 
squint angle and maximum and minimum pulse repetition frequencies. 





e 


Point target positions are chosen to allow image quality measurements 
across the whole image. In particular, point targets are placed at the 
boundaries of the processor range invariance regions to demonstrate that 
image quality parameters are not degraded by any implementation specific 


restrictions. 
3.2 Simulated Wave Data 


The simulated wave data is used to test generation of the image spectrum 
product, which is the final SAR FDP wave mode product. 


The wave data is approximated as a sinusoidal function characterized by 
its wavelength and direction. 


Several data sets are generated using different wavelengths and 
Girections so the averaging process and the rectangular to polar 
coordinates transformation during spectrum product generation can be 
tested extensively. 


3.3 Pre-processed SEASAT Data 


Pre-processed SEASAT data is used to demonstrate that the SAR FDP can 
process real SAR data with image quality parameters comparable to those 
produced by MacDonald Dettwiler’s SEASAT processor. 


The pre-processing includes shifting the SEASAT data to baseband 
frequency so that the raw data is represented as complex numbers (real 
and imaginary parts) rather than as real numbers. Furthermore, the data 
is converted to S-bit to resemble ERS-1 data and then copied to HDDT to 
conform with the ERS-1 downlink format. 


4.0 Performance Test Results 


The performance tests include image quality parameter measurements such 
as spatial resolution, side lobe level, Integrated Side Lobe Ratio 
(ISLR), and also throughput measurements. Typical results of the 
measurements are shown below. These results are compared with the 
values stated in the processor’s requirements specification. 


4.1 Spatial Resolution and Side Lobe Level 


The spatial resolution is the 3-dB width of the point target response 
measured in both range and azimuth directions. The side lobe level is 
the highest side lobe level in the point target response relative to the 
main peak level. Results from the simulated point target data are as 


follows: 











Requirement Measured 


3-db width in slant range: < 1.44 pixels 1.29 pixels 
3-db width in azimuth: < 1.66 pixels 1.35 pixels 
side lobe level in range: < -20.0 GB -24.5 dB 
side lobe level in azimuth: < -20.0 GB -24.1 cB 


Spatial resolution of the SEASAT image is also measured and compared 
with the image processed by the original SEASAT processor. Two SEASAT 
scenes, Goldstone and Vancouver, are used for the test. The result is 
as follows: 


SEASAT SAR FDP 
3-db width in ground range: 24.8 metres 25.2 metres 
3-db width in azimuth: : 22.0 metres 22.35 metres 


4.2. Integrated Side Lobe Ratio (ISLR) 


The integrated side lobe ratio is the ratio expressed in dB between the 
Side lobe power and the main lobe power of the response from a point 
target. The main lobe power is the total power inside an area centered 
about the main peak and extending three resolution lengths in the range 
and azimuth directions. The side lobe power is the total power outside 
the main lobe but within an area bounded by ten resolution lengths in 
the range and azimuth directions. 


The result of this ISLR measurement is -20.0 dB .while the stated ISLR 
requirement is that the value be less than -12.0 GB. 


4.3. Point Target Misregistration 

The point target misregistration is the misregistration of a point 

target’s energy from its expected zero doppler location caused by 

approximations in the SAR processor. Misregistration due to orbit data 

and terrain height error is not included. 

The results of the misregistration measurements are as follows: 
Requirement Measured 


Range misregistration + 2.88 pixels 0.47 pixel 
Azimuth misregistration + 3.32 pixels 1.56 pixels 

















4.4. Radiometric Linearity 


Radiometric linearity is measured by the correlation coefficient between 
output point target peak power regression and input point target power. 


The simulated point target area is generated with varying power so that 
the data cover the dynamic range of the ERS-1 input data. The 
radiometric linearity value is required to be greater than 0.95, and the 


measured value is 0.98. 


4.5. Product Localisation 


The SAR FDP is required to compute the geographical position of the four 
corner points and the centre point of a 100 x 100 km image. The 
geographical positions are expressed in terms of geodetic latitude and 


longitude. 


The computation is checked by comparing results with those obtained from 
an independent software program developed by ESA. Typical results are 


as follows: 
Latitude Longitude 


ESA program result 36.369 degrees 239.338 degrees 
36.369 degrees 239.337 degrees 


4.6 Throughput 


The SAR FDP is required to generate the following image mode products 
within one orbit: 


--three fast delivery image mode products (100 x 100 km scene); 
--two image mode noise statistic and drift calibration pulse products; 


--three image mode chirp replica products. 


If the SAR FDP is configured to process wave mode data it is required to 
generate the following products within one orbit: 


--150 wave mode intermediate products (5 x 5 km scene); 
--150 fast Gslivery wave mode products (image spectrum); 


--ten wave mode noise statistic and drift calibration pulse products; 


--150 wave mode chirp replica products. 


The throughput has to be sustained from orbit to orbit for up to ten 
consecutive orbits. Considering all ground station operational 

















restrictions such as data acquisition and HDDT positioning, the SAR FDP 
is required to generate one fast delivery image mode product in less 
than 28 minutes. The tests show that one 100 x 100 km image mode scene 


can be generated in 27.3 minutes. 


Similar consideration of the ground station operational rastrictions 
during low rate data processing shows that one chirp replica, one 
intermediate and one image spectrum wave mode product must be generated 
in less than 28 seconds. The throughput measurement during the test 
shows these products can be generated in 26 seconds. 


5.0 Operational fspects 


The operational aspects of the SAR FDP are tested by scheduling the 
processor to generate all possible combinations of products. These 
inciude all the products listed for the throughput tests plus extraction 
products such as raw data products, general header products, and 
instrument header products. All products are checked against the output 
product specifications in terms of their formats and expected values in 


the product headers. 


In addition, several tests are designed to test the error handling 
strategy of the SAR FDP. The simulated errors include the following: 


--errors in the command parameters (inconsistent parameters) 


--invalid commands 


--errors in the system software configuration (e.g., a missing parameter 
table) 


--missing lines in the input data 
--insufficient input data 

--overlapping and obsolete wave mode source packets 

--error in spacecraft time code 

--range line counter wrap around (caused by the 16-bit representation) 


The SAR FDF is verified to handle these possible errors correctly. The 
verification was done by examining the SAR FDP responses and also the 


generated products. 


6.0 Conclusion 


The SAR FDP as developed by MacDonald Dettwiler under the ESA contract 
represents the newest implementation of the well proven SAR range- 


doppler algorithm. 








The development’s most significant achievements are its processing speed 
improvement (27 minutes per image compared to 2 few hours) and its 
preservation of image quality in a ground station operational 
environment. 
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[Text] ABSTRACT 


This paper describes the Brazilian experience in acquiring, storing, 
processing and disseminating Remote Sensing data collected from 
satellites since the first Landsat in 1973 up to SPOT I in 1987/1988. 
INPE’s roles in operating, upgrading, specifying, designing and 
integrating the related systems are covered, giving an evolution history 
which can be useful to other developing countries in the field of Remote 


Sensing. 


1. INTRODUCTION 


Fifteen years ago a decisive action taken by the Brazilian Government 
Placed ocr country in the promising position of becoming one of the 
leaders .n the Remote Sensing field of activities in the Southern 
Hemisphere. Thus, at that time, presicential authorization was granted 
to install and operate in Brazil a LANDSAT (ERTS-1) Ground Station, the 
third in the world to become operational, after only Canada and United 


States. 


The initial availability of data from the LANDSAT series of satellites 
since 1972 and the operation of its own ground station from 1973 on gave 
to Brazil the required conditions to develop their Remote Sensing plans 


in a more effective manner. 














Today, fifteen years have passed and the activities and perspectives in 
Remote Sensing have undergone an unequivocal evolution giving rise to a 


SOlid Brazilian experience in the field. 


This experience is represented both by the practical know-how and 
maturity of specialists in the use of orbital remotely sensed data for a 
wide variety of applications, and by their consolidated experience in 
the operation of satellite ground stations as well as in developing and 
establishing new image generating and processing systems for satellite 


data applications. 


2. HISTORIC RETROSPECT 


It seems proper to review now the results obtained during this long 
range effort and their main consequences in conjunction with the 
governmental priorities and the perspectives at internal level. 


During this time two different (but complementary) actions were 
developed in Brazil on the Government Level: the RADAMBRASIL Project 


and the Satellite Remote Sensing Program. 


The RADAMBRASIL Project, initially established to perform the Land 
Resources Integrated Survey of an area of 1,500,000 km2 in the Amazonic 
Region, was gradually increased until, by 1975, to cover about the 
entire 8,500,000 km2 of the Brazilian territory, thus becoming the 
world’s largest airborne radar (SAR) data coverage project, 38 reports 
with cxcellent thematic maps in the 1:1,000,000 scale were generated and 
became importrant assets for Government planners both on the Federal and 
State levels. This event gave rise to a very important by-product which 
was the training of large number of technicians on the handling and 


assessment of non-conventional images. Until then, the Brazilian 
technical people were used almost exclusively to handle aerophotographs. 


Today, both the data generated by RADAMBRASIL project and a fairly large 
number of specialists can be found in the IBGE “oundation, the Brazilian 
Institute for Geography and Statistics, in Rio de Janeiro. 


A different strategy had to be pursued, however, by tha Satellite Remote 
Sensing Program. First, it was soon recognized the priority need of 


preparing the peopls for the then unknown techniques (even at 
international level) for generating and handling earth resources 
satellite data. Opportunities were then created to send Grazilian 
technicians abroad to be submitted to different levels of specialization 
(i.e., from short on-job-training courses all the way to MSo’s and 


PhD’s). 











Further, it was enough just to look at the characteristics of Brazil to 
infer that very shortly the country would turn into a very important 

user of the data generated by these earth resources satellites. The 
authorization for purchasing a complete LANDSAT Satellite Ground Station 
was then issued by the Government in 1972, guided by the urgent need of — 
better knowing the large country’s territorial area and its renewable 
and non-renewable natural resources, the need of checking on the 
progress of large projects ** remote areas, the need to access, ina 
global basis, the results of an agricultural policy, and the cost 
involved in obtaining data for theses objectives using traditional 


methods. 


Benefiting from the strategic location of Cuiaba Receiving Station in 
the continent, a huge amount of data of our own territory and also of 
several other countries in South America became then directly accessible 
to a growing and potentially active user community. 


During these past years a considerable experience was gathered by INPE 
technicians in the process of developing techniques and methodologies 
and also in the identification of problems which could be either solved 
or mitigated with the use of satellite images. 


In this period a very large number of specialists from governmental 
agencies, universities and private companies were trained at INPE, 
through the offering of specialization courses in several levels and 


graduate courses, MSc level. 


Such an effort was intended primarily to spread the Remote Sensing 
disciplines and techniques through the user community in order to awake 
them for this exciting new technology and its potential applications in 
many areas of earth resources management and monitoring of the 
environment. A well trained user community would surely promote the 
country’s economic growth using satellite technology. 


Using the experience gathered from monitoring, operating, adapting and 
enhancing sophisticated receiving and processing equipments and 
associated software to meet satellite specifications of the LANDSAT 
series and the requirements of a growing user community, INPE gradually 
became a competent R&D center able to develop fairly complex units 
mainly in the areas of receiving and dats processing stations as well as 
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image analysis systems. It is worth to mention that this experience 
enabled INPE to reach, in the early 80’s, the status of co-participant 
in the design and implementation of the new LANDSAT 4/5 Brazilian ground 


processing system. In the mid 80’s INPE, reaching a deeper knowledge 
and experience among its technical personnel, focused its attention to 
the French Satellite SPOT getting prepared to answer alone the 
requirements for the full integration of a ground processing system to 
handle SPOT data. 


Nowadays the industry, strengthened by the government policies, is able. 
to make and sell abroad, competitively, not only items but also systems 


developed by INPE. 


Unfortunately, the economic crisis which hit Brazil in the last years 
reached very severely the R&D Centers, INPE among them. Suffering from 
lack of financial support and struggling for survival from lack of 
financial support and struggling for survival INPE (and some other 
important government R&D Centers too) was forced to pursue application 
target areas which naturally should have been transferred to the private 
companies, for example, biding for the execution of projects using 
technologies already established. 


However, with the new Government installed in March 1985, opportunities 
were granted to review every policy in the country, including the Remote 


Sensing subject. 


Meanwhile, a simple inspection of what is happen*ng in the western 
countries led the observer to conclude that the satellite remote sensing 
activities are going now through the transition from a well proven 
technology to a commercial venture, with a very high probability of 
success. The new policy from the Remote Sensing in Brazil should 
certainly take these facts in consideration. 


Thus, in harmony with the aims and priorities of the new Brazilian 
Administration and guided also by international trends, INPE, from 1985 
on, revised its own policy and structure in the area of Remote Sensing 
with the firm objective of intimately interact with the user community, 
universities and private industry in order to attain a highly desirable 
technological union. This will end up positioning the satellite 
remotely sensed information as the decisive tool for the investigators, 
specialists, planners and managers to achieve daily, weekly, monthly, 
seasonal, annual or long term monitoring of changes in the country’s 
environment and resources. In particular, consolidating its union with 
the private sector, INPE is in process to promote, as said, the 
country’s economic growth using satellite technology. 


3. SIGNIFICANT STEPPING STONES 
3.1 Ground Stations and Instrumentation 


Once acquired its first Receiving and Processing Ground Station in 
1972/1973, few months after launching of LANDSAT-1, INPE, in a gradual 
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way, became capable of modifying and upgrading the station to handle 
data from new satellites of the initial series (LANDSATsS 2 and 3), to 
fulfill new operational needs and to improve overall efiiciency. The 
knowledge absorbed in this process allowed INPE the effective 
participation in the system design and analysis as well as software 
development and system integration of the new ground station, capable of 
receiving and processing the Thematic Mapper data from the LANDSAT-D 
series of satellites, to come, sharing these efforts with the supplier. 










These efforts were managed within what was called "INPE LANDSAT-D 

Project" which set up, from late 1980 on, several phased activities 
aiming the establishment of the new facility which would allow the 
Brazilian Remote Sensing user community to access the markedly improved 
capabilities of the new instrument. Not only that but it would also 

fulfill the increasing requirements of this fast growing community in 
terms of timely, accurate and reliable data on earth resources. The new 
ground station for TM then came to reality, starting in early 1984 the 


generation of Thematic Mapper products. 


This joint development program was fundamental in giving INPE the 
ability to modify, improve and expand its station as needed. New 
products have been placed available to the user community these last 
years and several extensions were made in the overall processing 
software to allow easier management, higher throughput and accuracy, and 


better quality of products. 


This ability is being also put to test, in the last two years, through 
the efforts for upgrading the receiving and processing systems to handle 
the data from the French satellite SPOT. The development, 
implementation and tests of SPOT software as well as full system 
integration activities are being totally carried out by INPE 


specialists. 


At this point it is most appropriate to mention the full support INPE 
has received from the Government in the procurement of the hardware 
subsystems necessary to upgrade their station for receiving and 
processing SPOT’ data. This financial support was provided despite the 
severe economic crisis the country is lately suffering, in clear 
recognition from the governmental planners of the priorities and 
importance of INPE’s role in fulfilling the country’s needs. 


INPE has also been developing equipments for meteorological satellite 
Stations and for Data Collection Platforms (DCP’s). In order to 
stimulate the use and fabrication of PCDs in Brazil long range efforts 
were engaged by INPE, reaching finally the successful development and 
testing of an ARGOS PCD prototype. This prototype was certified by 
ARGOS Service in France, end of 1983. Once certified the ARGOS PCD, 
INPE transferred to the Brazilian industry the associated technology 
endowing the country with this very effective tool for collecting data 


under remote conditions. 
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Efforts are also being conducted to test, certify and transfer to the 
industry PCDs compatible with GOES satellites, and a mini-station for 
directly receiving and processing ARGOS PCD data. The developments of 
both GOES PCD prototype and the ARGOS mini-stations have been concluded 
and the test and certification phases are in progress now. Worth to 
mention is also the development of the GOES-AAA format synchronizer 
capable of processing the new data format pertinent to the GOES-6 and 7 
satellites.-. A prototype has been already developed and transferred to 
the industry which is now delivering the units to the national agencies 
which operate similar meteorological GOES stations. 


Concerning the TIROS/NOAA Satellites, a software package has been 
developed to correct, in the digital domain, radiometrically and 
geometricaliy, the high resolution image data provided by the AVHRR 


instrument. 
3.2 - Image Generation and Distribution 


Table 1 and figures 1.1 and 1.2 summarize the volume of data (in 
photographic media and digital tapes) distributed to the users during 
these fifteen years. 


The data above indicates an unprecedented utilization of photographic 
images by the user community during the years of °77, *°78 and ’79. 
However, during these years a big percentage (over 40 percent) of 
dGelivereuU ocroducts were processed in the smaller scales 1:3,704,000 and 
1:1,000,000 (black and white). Only about 15 percent of the products 
were requested in the 1:250,000 scales. The profile of the requests. 
indicated a user community mainly oriented to the use of LANDSAT 
products in global studies or analysis in smaller scales. 


A very significant reduction in the user demand occurred particularly in 
the years of ’83 and ’84 as a result of the very restraining financial 
actions taken by the Brazilian private companies and governmental 
agencies in general great difficulties in reaching their desired annual 


budget. 
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TABLE 1 - IMAGES DISTRIBUTED TO USERS 


YEAR | IMAGES (PHOTO) j; ccTs | TOTALS 














The situation changed and in the last three years a clear recovery is 
being observed. Not only that but the user tendency is now markedly 
towards the larger scales 1:250,000 and 1:100,000 which account for more 
than 50 percent of the total volume of photo products distributed in the 
year. The interest for color products has significantly increased being 
now around 20 percent of the total. These reflect the growing maturity 
and knowledge of the user community which is accessing more deeply the 
greater potential of the remote sensing products and is using more 
sophisticated ways for monitoring the environment and managing of earth 
resources, taking advantage of the finer resolution and accuracy of the 


new sensors. 


Another aspect which deserves notice is the increase of the number of 
CCTs distributed to the users. This is a clear answer from the 
community to the efforts INPE is conducting to make more known and 
accessible to the users the image analysis systems called SITIMs, 


manufactured by the Brazilia.. industry. 


Brazil has an outstanding position among the countries operating the 
LANDSAT system. The data base includes nowadays nearly 200,000 LANDSAT 
images. The number of users (mostly institutions) is steadily growing 
and today exceeds 1,500, many of them being from outside the country. 
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Figure 1.1: Photographic Images Distributed to Users 


MomrdibitP ie 





24 9 6%OCUT?tCO OO OkCOktCkkS GC? 
Year 
Figure 1.2: Digital Images Distributed to Users 


During these fifteen years INPE has taken under its own responsibility 
the activities of marketing, selling and distributing the remote sensing 
images. In the last two years, a preparation is in progress to 
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gradually transfer to the private sector, particularly the Aeserosurvey 
Companies, the activities of marketing and commercialization of the 
data. This surely can improve the overall efficiency besides yielding 
better service to the user community. This is indeed most opportune, 
once SPOT data, which are being already handled in a commercial basis, 
are being acquired directly at the INPE Ground Station, and are in 
process to be placed available to the user community. 


3.3 - Cartography 


It is significant to mention the efforts developed by INPE since 1980 
towards the implementation of techniques for geometric correction of 
orbital imagery, mainly LANDSAT-MSS data, in order to transform the 
images into means more efficiently usable by the cartographic community. 


After TM data became available it was soon recognized their enormous 
potential for cartographic applications. The advanced characteristics 
of the orbital platform, the new geometry of the instrument and the 
precise knowledge of its behaviour and associated modelling, made T 
data of unprecedented value for cartography from space. Taking 
advantage of all these benefits, INPE specialists developed and 
implemented a processing station software capable of generating TM fully 
system-corrected products of very encouraging performance: internal 
geometric errors of the order of one pixel “rms,” with no use of ground 


control points. 


Benefitting, therefore, from special image treatments and from the 
existence of more advanced sensors, INPE has been directing its 
cartographic activities towards: a) production of planimetric satellite 
image-maps for future inclusion of altimetric information, 6) production 
of thematic maps in different scales, and c) updating of conventional 
topographic maps in 1:250,000 and 1:100,000 scales. An achievement 
which also deserves notice is the successful experiment of digital 
mosaicking of two adjacent TM fully system-corrected scenes (Barbacena 
image-map), in preparation for the creation of a 1:250,000 satellite 
image-map. This product is particularly important in poorly mapped 


areas of the country. 


Presently, INPE specialists are refining the operational procedures for 
the production of such image maps at 1:1,000,000 and larger scales. 


3.4 - Image Analysis Systems 


Intensive efforts havo been made by INPE to develop, in the last few 

years, a family of image analysis systems called SITIMs. These systems, 
based on Brazilian-made microcomputers, are used mainly in the analysis 
of natural resources, although their frames allow the use in other areas 


such as Metallurgy, Meteorology and Medicine. 


Four of these systems are presently installed at Image Processing Labs 
of INPE, supporting directly research and application projects in Remote 
Sensing and Meteorology. Twelve units are installed in other agencies. 
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Engineers and system analytts specialized in digital processing, 
modelling and simulation techniques have, since 1975, been doing 
researoh in enhancing, filtering, classification, registration, 
geometric correction and image analysis techniques. 


INPE is also working in the field of Geographic Information Systems, 
conducting a project where satellite imagery is being incorporated and 
associated to other related information plans (maps, digital terrain 
models, etc.) with a wide variety of image analysis techniques, 
everything manipulated through a user-friendly software package. The 
hardware concept is based on Brazilian-made microcomputers and the first 
version is finished and being submitted to an evaluation from the user 
standpoint. Finished this phase it will be also transferred to the 


private industry. 
3.5 - Environmental Analysis 


Qualified people of INPE in the field of Environmental Analysis have 
been conducting projects on techniques and methods for gathering remote 
sensing environmental data through the use of orbital and airborne 
imagery. These projects have been oriented to the following areas: 


a) Environmental Geomorphology with studies of erosion and equilibrium 
of inland water systems, | 


b) Hydric Systems with studies on water basin managing and water 
quality and availability analysis, 


c) Land use with studies on the soil occupation, and 
qd) Urban Areas with studies of their growth and use. 


3.6 - Geology 


In Geology, studies have been conducted for mining and energy (i.e., 
©1il, gas) surveys through the analysis of the land surface from not only 
satellite imagery but also airborne data and field conditions. It is 
known that to extract‘ geological information from remote sensing data is 
necessary a careful work by experienced geological interpreters 
knowledgeable of both remote sensing system characteristics and field 
conditions. The reason is that most of the observed surface phenomena 
are often related to geologic parameters in only a tenuous manner. 
Despite these difficulties success has been verified on various cases 
being worth to mention projects on: a) Mineral and Energetic Resources 

with mineral and oil deposit studies through remote sensing for the 
establishment of prospection models, 6b) Geological Mapping with 
structural analysis and identification of lithostratigraphic units in 
different scales seeking the definition of evolutionary models, c) 
Engineering and Hydrological Geology with remote sensing applications on 
geotechnical and hydrogeological problems. 
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3.7 - Forestry 


In the field of forestry INPE efforts have had the objectives of 
research and development of methodologies for utilization of remote 
sensing data in the effective monitoring and management of the forest 
resources. The main project lines in this area have been: vegetation 
mapping, forest management and inventory, estimates of timber volumes, 
detection and assessment of damages due to burns, diseases or droughts. 
A number of forestry agencies are also becoming increasingly involved in 
the development of techniques for applying remote sensing methods to 
environmental monitoring problems. 


3.8 - Agriculture 


The use of remote sensing techniques .to carry out projects on crop 
survey and potential yield of areas is of utmost importance for the 
national economy. In this framework INPE has conducted research and 
studies on wheat, soybean, corn, rice and sugar-cane. Estimates of 
planted and/or tilled areas are done through visual and/or automatic 
interpretation of satellite or airborne imagery, thus performing the 
integral or sampled mapping of regions of interest. These techniques 
can also be used to estimate the losses caused by plagues, diseases and 


periodic events such as frost or drought. 


Yield estimates are obtained through modelling with data from orbital 
remote sensing, eventually supported from airborne data and conventional 


meteorological studies. 
3.9 - Basic Research in Remote Sensing 


The efforts in this area have been conducted primarily to increase the 
knowledge about the remote sensing data and its possible correlations 
with terrain features. Research activities encompass the study and 
analysis of the parameters which affect data acquisition and 
registration, such as characteristics and conditions of terrain 
features, spectral and radiometric properties, atmospheric effects, 
sensor specifications and image geometry, among others. 


The efforts also deal with evaluation studies of the capability of 
certain sensors to well characterize certain target areas of the 
country. Under this point of view efforts are also being conducted to 
access the potential of SPOT data over the Brazilian territory. 


3.10 - Technology Transfer and Technical Orientation 


As said, the Brazilian remote sensing user community is steadily growing 
and INPE is vitally interested in promoting the establishment of a large 
team of specialists throughout the country in order to maximize the use 
of the advanced remote sensing tools and to join efforts with the 
universities and private companies in such way that a strategic bridge 
will be established between the space technology and the social and 


economical development of our country. 
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The efforts developed towards these objectives in the last years, 
particularly from 1985 on, can be summarized below: 


a) Graduate courses in remote sensing and image analysis, with regular 
courses beginning in February each year and leading to a MSc Degree. 


b) On-the-job training for professionals needing a familiarization with 
the remote sensing techniques. 


c) User assistance and orientation service with a specialized team, 
when applicable, INPE’s user-aid services units spread throughout the 
country are used to facilitate user access to the data or to the 
technique (presently there are centers in the following Brazilian 
cities: Natal, Manaus, Brasilia, Rio de Janeiro, Sao Jose dos Campos, 


Cachoeira Paulista and Campina Grande), 


d) International cooperation in remote sensing with other countries and 
international agencies. To the present INPE develops cooperation 

programs with Argentina, Bolivia, Chile, Colombia, Equador, Canada, GDR, 
France, FRG, Great Britain, Iraq, Japan, Panama, Paraguay, People’s Rep. 
of China, Peru, Uruguay, USA, FAO, UN, World Bank and the interamerican 


Development Bank. 


e) International training courses, for Latin American Specialists 
(1985-1986) and for African Specialists (1987-1988), being offered in 
cooperation with the U.N. Outer Space Affairs Division and supported by 
other agencies such as SELPER (Latin American Society of Remote Sensing 
Specialists), ESA (European Space Agency) and CNPg (National Council of 
Scientific and Technological Development of Brazil). 


f) Dissemination and establishment of regional remote sensing 
laboratories, organization of symposia, seminars, lectures, workshops, 
demonstrations, etc., on remote sensing have been very active in 

INPE in the last years. The goal is to familiarize the different 
professionals with the most recent developments in. the field. 


INPE is also working on another powerful way of transferring to the user 
community the remote sensing knowledge absorbed during this long range 
experience, whioh is the planning and installation of regional remote 
sensing laboratories. They are oriented, in general, to work guided by 
the characteristics and needs of each region. Regional Remote Sensing 
Laboratories are already in operation in Campina Grande, at Northeast of 
the country, and Belem, at the North. Eight new laboratories are 


already approved for the period ’88-’89. 


4. EXAMPLES OF RECENT APPLICATIONS AND ACTIVITIES ON REMOTE SENSING IN 
BRAZIL 


4.1 - Land Reform Program 


The Brazilian Land Reform Program which is being conducted by the 














Ministry of Land Reform and Agrarian Development includes among its 
activities the identification and selection of nonproductive land 
parcels for redistribution among small farmers. These activities 
include the identification of land suitable to establish new 
settlements. This process tackles with technical, social, political and 
legal issues. The way it has been conducted in the past has been 
subject of criticism and laid on political grounds. INPE developed a 
methodology using satellite remote sensing techniques and is helping the 
government in the process of land identification and selection of land 


areas for allocation. 


To date remote sensing techniques have been used to examine 1,000,000 
km2 in Brazil for the Land Reform Program. The following major 
conclusions may be mentioned regarding the application of remote sensing 


techniques in land reform issues: 


R.S. techniques have allowed to establish a technical approach to 
lessen the political aspects related with land redistribution. 


--R.S. techniques have provided for lowering the cost of land selection 
and inspection and have expedited the process of land identification and 


selection. 


--Satellite data have allowed multitemporal analyses which permitted the 
study of the historical evolution of the land use practice or degree 
during a time interval. This has provided a strong argument concerning 


the history of the occupation of the land. 


--R.S. techniques have permitted to unveil “make up” practices. Some 
farmers, knowing that their lands were subject to land act decrees, 
would fastly plow their lands and question in judge government land 
misqualifications or errors. Through the historical archive of 
satellite images it has been possible to show the status of the land use 
before and at the time of the land reform acts. 


--Government lawyers have been briefed upon remote sensing techniques, 
study areas, and provided with a kit of satellite images and reports on 
the historical land use. These space data have been present to the 
judges. Thus, remote sensing has been used now also as a legal tool. 


4.2 - National Irrigation Program 


Brazil has devoted considerable efforts in the development of its 
agriculture. In 1986 the federal government decided to promote an 
increase in grain production of the order of 20 percent, based on a 
National Irrigation Program. In order to reach that goal government has 
to apply the resources of the irrigation program in areas which have an 
agriculture vocation or a return potential. This required to know in a 
short period of time the actual distribution of the agriculture 
activities, a task that by conventional methods could take years. Based 
on this, INPE has developed a six-month remote sensing project which has 
covered an area of 3,500,000 km2. The resulting maps provide updated 
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information on the distribution of crops. Based on this information, 
government may now establish priority project areas and apply the 
resources more effectively. 


4.3 - Water Resources 


The Government of the Pernambuco State has been concerned with the 
drought in the semiarid Northeast of 6razil and has asked INPE to 
present alternative plans in 45 days. First it was identified by the 
governor himself that the first priority should be the man, that is, 
water for drinking. Drought periods in Northeastern Brazil causes 
migration of population from rural to urban areas or to more developed 
cities in Southern Brazil. During drought periods the government has to 
hire water trucks to distribute water to the population. INPE has 
developed a forty-five-days remote sensing project for the whole State 
of Pernambuco (91,000 km2) which has included 1:100,000-scale maps 
containing information on the distribution of alluvial-sediment areas, 
updating of surface water resources (dams, ponds, drainage, etc.), human 
activity concentrations, actual drainage (including intermittent 
drainage) and road networks, and existing water wells. 


Based on the maps provided to the government, a strategic and more 
efficient plan for water distribution might be made. Planning of the 
distribution of the water trucks and schedules (e.9., population versus 
pond distribution), planning for the establishment of public fountains 


and other actions may be taken. 


5. FUTURE PERSPECTIVES 


A challenging step has been taken by the Brazilian Government with the 
official approval of the so-called Brazilian Complete Space Mission, 
MECB. MECB has the objective of designing, building, integrating, 
testing, launching and operating four national application satellites: 
two for Data Collection and two for Remote Sensing. The launcher is 
being developed in Brazil by the IAE--Institute of Space Activities 
(Ministry of Aeronautics). The expected date for the first satellite 
launching is mid 1989. The first two satellites, i.e., the data 
collection ones, will use the Brazilian made DCP’s (Data Collecting 
Platforms) which are an economically viable option for national 
applications. The greatest advantage of the telemetry via satellite is 
that costs are not too high and are not related to the distances 


involved. 


The remote sensing satellites, whose launching is foreseen to occur 
early 1993 will relay information on natural resources evaluation, 
through moderate resolution and a highly repetitive global coverage, 
using an onboard Multispectral Earth Observation Camera. This will 
enable the country to achieve its own data bank at low cost and on 


nearly real time. 


The INPE’s opportunity to work as the designer and system integrator in 
the SPOT ground station turned out to be a necessary and very important 
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step towards the preparation for the Remote Sensing Mission of the MECB 
where both ground and space segments will be designed and built in 
Brazil by INPE with the help of the Brazilian industry. There are 
Similarities on both missions. For example, the experience gathered in 
the calibration, correction and processing of the charge-coupled- 
detector (CCDs) arrays featured by the SPOT imaging process will be 
shared by the MECB-Remote Sensing which is using the same process in its 


multispectral camera. 


Nowadays, aS it has been already mentioned, the policy among the western 
countries is the transference to the private sector of the space 
technologies, which have been fully tested and approved after years of 
operation. This has happened to the telecommunication satellites, and 
is happening now with the French Satellite SPOT which the French 
Government operates under commercial basis through the SPOT Image Co. 
LANDSATS 6 and 7, to be launched in the early ’90s, will also be 
operated in a commercial basis by the EOSAT Co. in USA. 


Drawn this picture, it becomes imperative the need for Brazil to possess 
its own satellites. 


A satellite configuration is chosen to attend the aims of its owners. 
In a domestic satellite, the spectral bands and the channels of the 
observation camera are established as a function of parameters the 
country wants to observe, measure or collect. These actions will allow 
the establishment of long range economical and social programs, without 
incurring in the heavy satellite fees often imposed by the satellite 
operators abroad as well as gaps in the data availability due to 


political or economical reasons. 


The definition of the Remote Sensing Mission of the MECB was then 
conducted in accordance with the above aspects. Aware of the present 
constraints of the satellite and taking them into account the remote 
sensing specialists of INPE have come to a proposal of certain mission 
characteristics in order to better fit it, as much as possible, to known 
needs in the areas of research and application of the Remote Sensing in 


Brazil. 


A technically oriented inquiry was also conducted by these specialists 
to evaluate closely the needs of Brazilian user community in terms of 


remote sensing data. 


As a result of this joint analysis it came up a mission with the 
following main characteristics: moderate spatial resolution, high 
temporal resolution with larger distance swath and a sensitivity in both 
the red and near-infrared regions of spectrum, most oriented to studies 
where the target phenomena are characterized by a high dynamics 
especially in the areas of Agriculture, Environmental Control and 


Oceanography. 


Still in the international scenario, an ever growing number of new 
satellites and data acquisition systems will be placed in orbit, so that 
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satellites and data acquisition systems will be placed in orbit, so that 
natural resources surveying processes and environmental monitoring 
techniques will become more and more simple, precise and less expensive. 


So it is time to plan ahead, thinking in terms of the future. 


Brazil has today thousands of Landsat images of its territory and a 
complete radar coverage (SAR) is available. Besides, it has data 
gathering and processing facilities which may be upgraded at reduced 
costs to handle future satellite data. More important, it has quite a 
reasonable number of skilled technicians who can make adequate use of 


the available technologies. 


However, it is high time to bring in the private industry. Our 
electronics and informatics industries are capable of producing the 
items required by the remote sensing specialist. Besides, the 
aerosurvey companies can also be present in the data commercialization, 
thus improving the overall efficiency and yielding better service to the 
users community. The Government would stay in the area of the research 
of new technologies, improving the existing facilities (e.g., installing 
Regional Remote Sensing Laboratories), training personnel and warranting 
the data access to new satellites and systems of interest. 


This way the base which will allow the Brazilian society to have an 
effective participation in defining future Space Missions, either 
National or international, will be established. 


6. REFERENCES 


Barbosa, M.N., Bastos Netto, D., Cunha, R.P.: Overview of the Brazilian 
Satellite Remote Sensing Program and Selected Examples of Recent 
Applications. 38th International Astronautical Congress of the IAF, 


Brighton, United Kingdom, Oct. 11-17, 1987. 


Barbosa, M.N., Bastos Netto, D.: Brazilian Satellite Remote Sensing 
Program. International Conference Remote Sensing for Development, 


Berlin, East Germany, Sept. 1-7, 1986. 


Parada, N.J.: Brazilian Remote Sensing Receiving, Recording and 
Processing Ground Systems in the 1980’s. 18th International Symposium 


of the Environment, Paris, France, Oct. 1-5, 1984. 


/ 08309 


as 
es 


23 

















JAPAN 


ESA’s Activities and Plans for Remote Sensing Data Systems 
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[Article by Dr. Livio Marelli, ESA/ESRIN, Frascati, Italy] 


{Text] 1. INTRODUCTION 


The European Space Agency (ESA) is involved in Earth Observation from 
space since the mid 1970s. This has been through the Meteosat 
programme, operated now on behalf of Eumetsat, and through the Earthnet 
programme, which up until now have acquired, archived, processed and 
disseminated data from non ESA remote sensing satellites like Landsat, 
SPOT, MOS-1, TIROS-AVHRR (Advanced Very High Resolution Radiometer ) 
plus, in the past, Seasat, HCMM (Heat Capacity Mapping Mission) and 


Nimbus-7. 


At present ESA is building ERS-1 which is an ocean and ice monitoring 
satellite due for launch in 1990. The Earthnet programme is upgrading 
its facilities to meet this new challenging mission. It will guarantee 
service not only to the science and research community, but also to 
global and regional weather and sea state forecasting organisations, 


offering near real time sea state information. 


Recently, the Long Term Plan of ESA has been approved (up to the year 
2000) and within the plan, a comprehensive Earth Observation programme 
is included. The plan involves, in particular, a second flight unit of 
the ERS series and the participation of Europe to the Polar Platform 
complex developed in the framework of the International Space Station 


(ISS) programme. 
Within each one of the phases of ESA’s activities in remote sensing 


mentioned above, the problem of data systems have, and will be, central 
and will involve major technical and operational challenges. 
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This paper outlines the programmatic approach followed so far, and 
describes the evolution planned for the coming generation of ERS 
satellites as well as the preliminary end-to-end system concept 
envisaged for the Polar Platforms later in the 1980s. 


2. THE PAST EXPERIENCE 


Up until very recently remote sensing has been mostly an area of 
research and proof of concept: missions like Seasat, HCMM, Nimbus-7, 
and to some extent Landsat Multi Spectral Scanner (MSS) and Return Beam 
Vidicon (RBV), were aimed at demonstrating the feasibility of specific 
sensing techniques for a number of research/application domains. 


The satellites mentioned above were not developed by ESA and had no 
assurance of continuity: in fact, Landsat is the only service which has 


maintained practical data continuity over the years. 


Another aspect of the series of satellites launched towards the end of 
the 70s was the fact that they relied mostly on regional ground stations 
to collect data since only some had on-board recorders, the vast 


majority of which proved rather unreliable. 


ESA’s philosophy for the data systems associated with the above missions 
can be summarised as follows: 


--Decentralised network of ground stations covering the major areas of 
interest for the European users (e.g., Europe, Scandinavia/Greenland and 
the Arctic region, North Africa and the Sahel). 


--Each regional station charged with data acquisition, archiving and 
pre-processing (bulk correction) on request. 


--Centralised network management and users interface service including 
catalogue, quick-looks, order handling and quality control. 


-~-Dissemination system, exclusively off-line delivering digital and 
photographic products through mail and/or special delivery to nominated 
centres in participating countries or to individual users. 


The archiving task is restricted in preserving the original raw data 
recorded digitally on High Density Digital Tape (HDDT) (one should bear 
in mind that Seasat Synthetic Aperture Radar (SAR) and Landsat RBV were 
transmitting in analog form) or on Computer Compatible Tape (CCT) for 
HCMM/Nimbus-7. 


¢ 
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Data received was screened systematically to generate quick-looks and 
catalogue entries primarily to obtain reliable cloud cover information, 
although the cloud cover assessment was assigned by operators with 


computer support. 


The catalogue of such missions as Landsat was assessable on-line from 
low speed land lines or packet switching networks. 


A full set of quick look prints of data available was set up at the 
centralised user services facility (ESRIN, Frascati) as well as at most 
of the national centres which act as points of contact for their 


national user community. 


Products were generated exclusively upon user request in view of the 
large difference between data collected and data actually exploited: it 
should be mentioned that, in the case of Europe, up to 70 percent of 
acquired data are cloud covered for a good part of the year. As regards 
acquisition strategy, two approaches were adopted: 


a) Landsat data was acquired systematically over the coverage of member 
states and to a lesser extent over the Arctic due to swath overlap. 


b) The experimental satellites were handled instead according to an 
agreed science plan defined between Principal Investigators (PIs), NASA 


and Earthnet. 


So far the system has previded an acceptable service and the large 
archives have proved accessible longer than expected though the old 
Landsat data (’75-’78) cannot be guaranteed any longer as the rest. 


All attempts of purging the data sets of cloud covered passes have 
demonstrated that the exercise doesn’t pay off because, transcriptions 
are labor intensive and strain the old HDDTs which are difficult to 


maintain after so many years. 


On the other hand, user demand for very old data is so modest that it 
doesn’t justify major investments. 


For several years, quick-looks for Landsat MSS were compacted on 
microfiches but also this exercise proved costly and of limited value to 
users due to long delays between sensing and the availability of the 


microfishes. 


In the specific case of Meteosat, the handling of the archive on HDDT 
proved infeasible and it was decided to transfer data on 6250 bpi CCTs: 
the volume of data involved and the importance of the data set for such 
projects as FGGE made the exercise worthwhile. Figure 1 gives some 
general indication on volume of data handled so far by Earthnet. 
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3S. THE PRESENT GENERATION OF REMOTE SENSING MISSIONS 


Starting with the launch of Landsat 4/5 in the early ’80s and later with 
the launch of Spot 1, a new generation of pre-operational/commercial 


missions have become available. 


Furthermore, the remote sensing user community became more fully aware 
of the value of the AVHRR of the TIROS-N series which though designed 
originally for polar meteorology, proves very valuable for a number of 
other applications (vegetation index, etc.). 


Lastly NASDA, in Japan, launched their first remote sensing satellite 
MOS-1 which Earthnet is acquiring with its network of ground stations. 


The peculiarities of the mission mentioned dictated some evolution in 
the system used with ESA which will be briefly reviewed here after. 


The overall philosophy adopted for Landsat Thematic Mapper (TM) and MOS 
is not significantly different from the one quoted in the previous 


chapter with few exceptions: 


_=--The quick-look system is undergoing modifications to allow generating 
of colour masters which are expected to improve users ability to 
determine suitability of images for their applications. 


--The network of ground stations is permanently linked via leased lines 
to the Central User Service (CUS) facilities which allows daily updates 
of the catalogue and reception of stations logs including information 
about data acquired, quality assessment reports on products generated, 


etc. 


--Experiments were successfully carried out disseminating within 6 to 8 
hours from acquisition MSS full resolution scenes over areas of interest 
to a number of pilot/demonstration projects using high speed satellite 


link. 


Spot operations in Maspalomas (Canary Islands) are carried out with the 
same basic philosophy, however, since the access fee is linked to the 
number of scenes requested acquisitions are not systematic but scheduled 


to match resources and user demands. 


At present, Maspalomas can only acquire and record data while 
processing, when requested by users, is carried out by Spotimage in 
Toulouse. Later on this year a standard Spot processing chain will be 
installed at the station which will then serve users directly. 
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The TIROS-N/AVHRR and TOVS (Tiros Operational Vertical Sounder) data are 
part of a newly started activity of Earthnet aimed at preserving those 
data which are widely used in Europe for meteorological purposes but are 
thereafter not usually kept for long term exploitation. 


The approach adopted here is somewhat different: having recognised that 
many stations are already in existence, Earthnet has negotiated 
agreements with national High Resolution Picture Transmission (HRPT) 
operators to host at their premises a data archiving subsystem which 
extracts the data to be preserved, derives catalogue parameters 
including cloud cover, sunglint, etc., generates a more advanced quick- 
look and evaluates, based on specific criteria, whether or not data 
should be preserved in a long term archive. In view of the more 
manageable volume of data involved, archiving is carried out on optical 
disks (12", 2Gbytes disks). The optical disk generated at the various 
HRPT stations will eventually converge onto a centralized archives where 
historical products requested by users will be generated when required. 


A central computerized catalogue connected with the archival stations is 
available on-line from ESRIN. Another feature of interest in the TIROS 
system is that the archival strategy is decided centrally and the 
system, based upon the planned HRPT operations and on actual cloud cover 
evaluation, results will automatically direct station to or not to 
archive data so that specific cloud free data sets over area of interest 


be secured. 


The option of generating automatically at the stage of data archiving, 


level 2 products like vegetation index over land, sea surface 
temperature over sea, etc., is being evaluated in cooperation with major 


user groups. 
4. THE ERS-1 PAYLOAD DATA SYSTEM 


ERS-1 is the first major remote sensing mission developed by the 
European Space Agency. 


ERS-1 has scientific and application objectives: for this reason ESA is 
committed to provide, within three hours from sensing, geophysical 

parameters associated with global sea state (wind fields, wave high and 
wave image spectra, etc.) plus high resolution regional data over Europe 


derived from the SAR. 


ERS-1 payload data system has been designed as an end-to-end complex 
which includes the following main components: 


--a number of Real Time Acquisition stations charged with the task of 
acquiring from the on-board recording system over an orbit and 
processing it to geophysical level and make it available for 


dissemination. 
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--a number of SAR acquisition stations capable of generating few SAR 
scene in near Real Time. This service is exclusively available for the 


European coverage. 


--a centralized user services centre (based at ESRIN, Frascati) in 
charge of handling all interfaces with the users, cataloguing, order 
handling, quality assurance network management and, upon delegation, 


mission management. 


--a number of archive and processing facilities in charge of preserving 
data generated by ERS and generate products for the users either on 
request or by retrieving relevant products which had been generated by 


systematic processing. 


--a wide band data dissemination system ensuring the distribution of 
near real time data products to users in nominated centres of member 


states. 


--a number of national or foreign acquisition/processing centres having 
signed agreement with ESA to have access to ERS-1 data. | 


The salient feature of the ERS payload data system are: 


--all ERS raw data are going to be preserved for at least 10 years after 
the lifetime of the satellite. Figure 2 summarizes the anticipated 


volumes involved. 


--global Low Bit Rate (LBR) data will be systematically processed and 
disseminated in near real time while SAR products will be generated and 


disseminated exclusively upon user request. 


--the archival task will be carried out in a decentralized way 
separately from the data acquisition facilities. Processing and 
Archiving Facilities (PAFs) have been selected as centres of expertise 
in different application domains of ERS. These PAFs shall, in 
particular, support the Agency in the calibration and validation of the 
mission, its sensors and products, in the definition of novel or 
improved or additional products to be extracted from ERS either off-line 
or in near real time as well as in the definition and validation of 


quality assessment methods and procedures. 


--the central catalogue shall provide information about data already 

acquired as well as on acquisition plans: likewise the users shall have 
the capability of placing orders for products from the archive but also 
for data yet to be sensed and for products to be made available in near 


real time as well as off-line. 


The catalogue shall contain reference to all data collected for ERS-1 
including those acquired by foreign stations and users will not be 
bothered to know where physically data they need is kept or production 


is carried out. 
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At present it is not anticipated to generate the equivalent of Landsat 
or Spot quick looks: it is possible that such a service might be 
developed for the Along Track Scanning Radiometer (ATSR) which is an 
experimental sensor provided as an instrument of opportunity but it is 


not planned for the SAR. 


ERS-1 will represent a major challenge as regards the payload data 
system and will allow development and validation of a number of novel 
features which will be applied for the Polar Platform complex. 


5. THE POLAR PLATFORMS 


The polar platform complex is a part of the international space station 
programme which is being finalized between ESA, Japan, Canada and 


Europe. 


The polar platforms will carry on-board a very large number of sensors 
of scientific and application nature and they are aimed at a long term 
earth monitoring system lasting well into the next century. 


The polar platforms will rely on a complex space infrastructure 
developed as part of ISS which will greatly enhance the system and will 
also impose significant changes to the payload data handling design and 
implementation: typical examples are the Data Relay satellites (TDRSS, 
European DRS and Japanese DRS) which will allow real time world wide 


access to data from the platforms. 


At present ESA is only studying the Polar Platforms’ (PPFs) end-to-end 
data systems. Therefore, only very preliminary considerations on its 
final set up can be provided to date. We anticipate that in a number of 
aspects the PPF payload data system will be similar to the one presently 


in place, or that foreseen for ERS-1 e.g.: 


--the system will be decentralized making use of several regional 
facilities as well as of specialized centres in charge of specific 
services along the lines of ERS PAFs. Such facilities however, will not 
be operated necessarily just off-line but could receive global data 


through data Relay channels. 


--specialized user community like meteorology could and probably shall 
set up dedicated facilities to handle sensors data of specific interest. 


--the centralized user services will offer information about data 
available for the PPF complex rather than for individual missions. 


--near real time as well as off-line services will be ensured with 
systematic and “on request” operations according to the product type and 


associated user requirements. 


Conversely, a number of new aspects will have to be addressed and solved 
with novel approaches: 
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--the central user service for Europe will have to integrate user 
requests for instrument operated on non-European platforms and liaise 
with the relevent mission management centres to ensure their feasibility 


and actual implementations. 


--data collected in the United States or in Japan needs to be 
transferred to Europe when required, either in near real time or in an 
off-line mode, be processed and delivered to the requesting users. This 
data traffic may well prove to be complex and expensive and needs 


therefore to be optimized. 


--the data dissemination systems of the major partners covering NORTH 
AMERICA, EUROPE and Japan might need to be interconnected. 


--standardization and intercalibration of products will become an 
essential requirement to allow users to handle indifferently data 
generated by different sensors operating on different platforms and 
received and processed by different partners in the system. 


--the different nature of the sensors operated on PPFs being either 
scientific/application of commercially orientated will imply complex 
problems of data allocation, protection, accounting, etc. 


The challenge is major but this path appears the only viable one to 
address the complex issues of environmental monitoring or global change 


studies. 


6. CONCLUSIONS 


Remote sensing systems seems to change generation every 5 to 6 years: 
Landsat-1 in 1972, Landsat-3, Seasat, Nimbus, HCMM in 1978, Landsat- 
5/Spot in 1984/86, ERS, JERS, TOPEX 1990-92, PPFs in 1995-97. 


Within such a time span two major aspects evolve: 


1) User requirements 


2) Technology available 


We have attempted to describe how in Europe the transition from 
experimental to pre-operational space remote sensing system is taking 
place, and its impact on data systems. The technology evolution is 
obviously playing an equally important role but unfortunately, each 
generation of remote sensing systems operates approximately with the 
technology prevailing during the previous generation and little can be 


done to modify this pattern. 


So far the most glamorous evolution in remote sensing payload data 
handiing is in the area of available computing power and on-line 


storage. 
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High speed recording technology has not yet evolved at the same pace. 


Archival technology based on optical recording is reaching maturity only 
in these last years and has a long way to go before adequate 
standardization is enforced. 


High speed data dissemination and networking are becoming increasingly 
available and are likely to impact on traditional remote sensing 
services ina significant way. 


Remote sensing is on the move to become a major resource management tool 
in the next decade and beyond. 


Future is indeed a big challenge in this area. 
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FUCINO 247366 63447 8600 - 7 
KIRUNA 177657 68472 17300 ~ - 
MASPALOMAS 29089 3896 13600 1998 2134 
TROMSOE - - vile 465 - 
LANNION - ~ - - 1291 
TOTAL 454112 135815 46800 2463 3425 





Total Gbytes 
One MSS scene equals 35 Mbytes MSS 15894 
One TM scene equals 256 Mbytes T™ 34760 
One MOS-1 scene equals 16 Mbytes MOS-1 748 
One TIROS pass equals 70 Mbytes TIROS 172 
One nimbus pass equals 90 Mbytes NIMBUS 308 
TOTAL 51882 


Fig. 1. Scenes Acquired by the Earthnet Stations up to Spring '88 
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EXPECTED GLOBAL VOLUME PER LBR INSTRUMENT PER YEAR 


Instrument Raw Data FDP + IP 
Wave 870 Gbyte 302 Gbyte 
Wind 472 Gbyte 3 Gbyte 

RA 72 Gbyte 2 Gbyte 
ATSR 600 Gbyte ~ 
TOTAL LBR 2014 Gbyte 306 Ghyte 

A Rie He Fils Se te te a oe S2saeene 








POSSIBLE TOTAL VOLUME FOR SAR PER YEAR FOR ALL 
REGIONAL STATIONS 





Instrument Raw Data © FDP 
SAR 28000 Gbyte 990 Gbyte 


Fig. 2 ERS-1 Yearly Data Volumes 
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The Possibilities of Application of Satellite Data for Monitoring Ice 
Condition in the Antarctic Region 


43070007d Kyoto Selections from INTERNATIONAL ARCHIVES OF PHOTOGRAMMETRY - 
AND REMOTE SENSING in English Vol 27 1988 Part 68 pp 1-6 


{Article by N.N. Shvyrkov, All Union Research Institute for Marine 
Fisheries and Oceanography; V. Krasnoselskaya 17a, 107140 Moscow 8-140 


USSR; Commission Iv] 
[Text] ABSTRACT 


Based on satellite observation in visible and microwave ranges the 
possibility of investigation of ice formation process and ice melting in 
high latitude zones of the World Ocean is considered. The modern 
methods of data acquisition of the borders of drift ice, size of ice 
fields, age of ice formations are analysed. The methods of construction 
of weekly ice charts of the Arctic and Antactic regions resulting from 
satellites, aircrafts and vessel observations are considered. The 
analysis of ice distribution in the Antarctic region in light period of 
1980-1984 seasons is made. Annual and seasonal dynamics of ice cover is 
estimated. The results obtained are of practical and theoretical 
interest for ecological studies and for the investigations of 
possibilities of exploitation of live resources of the Antarctic region. 
The data available can be used for the development of some methods for 
prediction of ice conditions in the South Ocean. 


The processes of ice formation and ice destruction are the main elements | 
of oceanographtec regime in the high latitudes of the World Ocean and in 

the seas of temperate latitudes. Ice is the major factor which produces 
the effect on the security of navigation and performance of various sea 
operations. The successful fishery in freezing seas is possible only on. 
the basis of good understanding of the ice situation. The knowledge of | 
the peculiarities of ice regime is necessary in plannins and performance 

of the scientific research in different almost inaccessible region. 
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Sea ice can be used as a tracer 1n studying large-scale oceanic 
circulations, streams, rotations anc other sea water dynamic features. 
The quantity of ice in partially freezing seas is the mair index of 
winter type and good predictor of some productive processes in 


forecoming spring period. 


The process of melting ice in the Antarctic region provides the complex 
of hydrometeorological processes at the edge of drifting ice and 
icebergs. The production of phytoplankton in these regions can be 
compared with some highly produstive zones of the World Ocean. High 
abundance of phytoplankton at melting edges of drift ice is responsible 
for considerable formations of macrozooplankton and various fish 
species. Therefore, the investigations carried out aboard the R/V 
"“Gizhiga” in the Atlantic sector of Antarctic region revealed a group of 
icebergs around which at the distance of 50 miles a concentration of 
macrozooplankton was observed. The total biomass of thjs concentration 
was estimated at about several hundreds of thousand tons. 


Sea ice serves as a habitat for numerous marine animals such as penguins 
and crabeater seals. In high latitudes the ice distribution and shift 
produces a considerable impact on marine animal migrations, promoting 
the distribution of some species and restricting the distribution of the 


others. 


Due to its physical properties ice is considered to be one of few 
parameters of marine environment which can be easily detected by remote 
sensing. Ice formations can be observed in visible, infra-red and 
radio-range waves of the electromagnetic spectrum. Each method 
mentioned above has its own advantages and limitations. 


At present TV-images, obtained from meteorologic satellites or ‘Meteor’ 
are used in the Soviet Union for operative deciphering. The application 
of these images allows for the determination of large ice-floes, 
icebergs, ice-fields, the locations of channels and polynias, ice unity 
and age with a high degree of accuracy. On the basis of this 
information the Hydrometeorological Service of the USSR regularly 
constructs 5-day charts of ice situation in Arctic and Antarctic region. 


After geographical bridging of some TV-images from satellite the 
composite photographs in the polar stereographic projection (the scale 


1:15,000,000) are made. 


By deciphering the borders of drift ice, the location of fast shore ice, 
the position of ice-fields, borders of various density and age, the 
location of channels and polynias are determined. The data obtained 
from the analysis of satellite information are refined and supplied by 
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some information from shore stations, vessels and aircrafts. The 
refined information is transferred to the blank of the charts of Arctic 
and Antarctic region in the composite photograph scale and projection. 
Weekly charts serve as the synthesis of satellite information, 
observation of aircraft, vessels and shore stations. Since 1985 for the 
construction of charts the radar data obtained from satellite “Kosmos- 
1500" was used. The satellite equipment allowed for the acquisition of 
data on motion and ice distribution irrespective of the condition and 
synoptical situation. On the radar images the ice fields and icebergs, 
vast extensive zones, shelf and continental glaciers are easily 
differentiated. The series of several radar images of the same aquatic 
area allows for the detection of general drift of sea ice, for 
estimation of the dynamics of large polynias and for the determination 
of the tendencies of ice compression and rarity. Figure 2 demonstrates 
the composite of the Antarctic ice cover. This figure is made on the 
basis of radar images of satellite “Kosmos-1500." Various ice 
formations are perfectly distinguished on the images. 


Weekly charts are used for the estimation of seasonal dynamics and 
annual changes of ice conditions and for the determination of the 
extreme borders of the ice cover development. Ice charts serve as the 
basis for the estimation of ice drift rate and can be used for various 


predictions. 


Recently the Soviet Union as well as the other countries performs the 
studies on biological resources of the South Ocean. The investigation 
of the biological structure of the South Ocean demonstrated that 
existence of the whole complex of sea organisms functioning on the basis 
of macrozooplankton abundance occurring in the high.’ productive pelagic 
zone combined with the natural zone of drift ice. Vise viability of the 
whole complex of these organisms are found to depend on ice conditions, 


which becomes more intensive in summer. 


Based on satellite data the monitoring of ice conditions in high 
latitude zones was carried out. The charts broadcasted by 
Hydrometeorological Service of the USSR are received by the vessels and 
by different scientific institutions. The charts thus obtained are 
resistered by graphical table on the magnetic disks and demand future 
processing. In particular for the performance of monthly analysis of 
the ice processes in Antarctic region the charts of the average monthly 
and extremal borders of ice cover development both in the whole South 


Ocean and in some region are made. 
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Fig. 3 ayy of ho ~y A of Antarctic 
on covere t ice in 
1980-1984 r : 


o—« 1980-81, o—« 1981-82, o—« 1982-83, —— 1983-84 


The location of the edges of the drift ice is compared with the average 
year-to-year data. The aquatic areas of waters covered by ice are 
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calculated. The degree of accessibility of some regions is also 
evaluated 1n seasonal and year-to-vear aspects. Figure 3 demonstrates 
the graphical changes of the areas covered by drift ice in 1980-1984. 
As a rule the maximal development of ice-cover is observed in October. 
The edge of the drift ice 1s approximately located at latitude of 60 
degrees S. The area covered by ice was 17-18 mln.km2. When the 
Antarctic spring begins the ice area gradually reduces. The most 
intensive destruction of ice cover is observed in late December-early 
January. The intermediate maximum observed in late November and 
December 1n all years are of interest. S8y early February the area of 
the drift-ice reaches its minimal values and composes 2.5-3 mln. km2. 
In the whole the edge of fast shore and drift-ice almost coincides with 
the slope of slightly mobile ice of the Antarctic regions. By late 
March a new period of ice formation begins. By the end of observation 
(in mid-April) the ice area doubles and reaches 6-7 mln. km2. In spite 
of the common tendencies in the process of ice formation in the 
Antarctic region in some years the several variations in the rate of 
ice-formation and ice-melting were observed. Thus the presented diagram 
illustrates the considerable variations in 1983-1984 as compared with 
previous seasons. It is known, that the first half of the season the 
ice conditions were more severe. In January and February the area of 
water was covered by ice by 2-3 mln. km2 more than usual. In the second 
part of the season the tendency changed. Figure 4 illustrates the 
extremal boundaries of drift ice in light period of 1982-1984. 
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Pig. 4 Extremal boundaries of drift-ice in s 
wow 1982-83 and —— 1983-84 ererw 


In 1983-1984 the maximal development of the ice cover was observed in 
October-November. From the beginning of October till late November the 
ice area reduced from 18 to 13 min. km2, which was by 2-3 mln. km2 less 
than observed last year. The most significant changes in the location 
of the bordes of maximal development of ice cover were observed in 
central parts of the Indian and Pacific Ocean sectors. The average 
monthly edge of October 1983 was located 2-4 degrees southward than that 
in last year. From the beginning of December 1983 till middle of 
February 1°84 the ice area reduced from 13 to 4 min. km2. During these 
months more severe ice condition were registered. In December 1983 the 
ice area was by 2-4 mln. km2 more than last year. By the end of March- 
early April 1984 the ice area extended to 5.5 min. km2. In this period 
in 1983 the process of ice formations was more intensive. By the 
beginning of April 1983 the area of drifting ice was 7.5 min. km2. 


All the above mentioned examples show that the satellites provide the 
valuable information on distant regions. The results of satellite 
observations are of great practical value for researches. The material 
obtained can serve as a basis for improving of the prediction of ice 
conditions in the Antarctic region. 
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Sea Spectra as Determined From Low-Altitude SLAR Images 


43070007e Kyoto Selections From INTERNATIONAL ARCHIVES OF PHOTOGRAMMETRY 
AND REMOTE SENSING in English Vol 27 1988 Part 88 pp 58-67 


{Article by Konstantin Voliak, General Physics Institute of USSR Ac. 
$ci.38 Vavilov st., 117942 Moscow, USSR, Commission VII] 


[Text] ABSTRACT 


The analytic composite theory for the surface of Gaussian slope 
statistics has revealed that backscattering of vertically (V) polarized 
microwaves from the sea at large angles of incidence might have an 
azimuthal characteristic, which corresponds to the angular dependence of 
the small-scale Bragg spectrum, with an accuracy to within the slope 
variance of long waves. For horizontally (H) polarized radiowaves, the 
model relates the image speckles to the statistics of spikes in a random 
surface-wave field. Fourier analysis and photometering of radiocoimages 
taken over various seas have allowed us to estimate the spectra of 
large- and small-scale sea surface waves as well as to study other 
phenomena imprinted on surface waves. 


1. INTRODUCTION 


Imaging radars yield principally novel data on the physical state of the 
sea. The aim of this paper is to develop simple models for image 
formation by an airborne real-aperture side-looking radar (SLAR) at 
different polarizations as well as to elaborate experimental methods and 
data processing techniques to effectively extract sea parameters, 
particularly wave spectra, from radicimages. 


2. THEORY OF MICROWAVE SCATTERING FROM THE SEA SURFACE 


Microwave backscatter from highly conducting surfaces, including the 
sea, depends strictly on polarization. At fairly large angles of 
incidence, 8, the H-polarized backscattered Gragg power is notably 
weaker than the V-polarized one. Generally the V-polarized power (Fig. 
la) shows up in a scattered field of the first order, which is rather 
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smoothly modulated by sea surface features. On the other hand, an H- 
polarized backscatter (Fig. 1b) produces an image formed by randomly 
spaced bright spots,. or by the so-called speckles, whereas the Bragg 
return is virtually at the noise level. A similar temporal behavior of 
polarized scattered fields was first established by Long. (1) 


The composite theory allows us to decode V-polarized SLAR images of the 
sea. Assuming that the conventional specular component and the second- 
order scatter, i.e., depolarization, are negligible at @ >30 degrees, we 
can write a local cross-section of V-backscattering[2): 


oCe.a, (2 = kfeot*e [o|* weak gine’ ak cose’ sing, C12. 


where ke * 2nn~* cosé ; = cosé’cosf?, 6°= Ota, a and #2 are the 
angles between the smoothed surface and the horizon in the plane of 
incidence and in the plane normal to the latter, € is the dielectric 
constant, W (K,K) is the roughness spectrum, the x-axis lies in the 


plane of incidence. 
2 ea a Z2ae 
6= (eCi + sin'@) - sin*@ IT sin*e’cos"p + T sin’, C2) 
ose -2 a _ 
A Ce 1)Cecos@ + {e-sin*é > a Ce- 13¢cosq 7 Je-sin"6) a 


v 
The average cross-section, %, , is determined by an integral over all 
a 


the slopes within an irradiated area: 
oscar ={foce.a,mopctan a, tan MdCtan wectan 9, 9 


where, neglecting the shadowing factor, 


Ktan a, tan = p Ctan a, tan "cose, /eone@ , c4)> 

‘ 

actan a, tan 7’) is the anisotropic probability density of sea slopes 
along its principal axes, the coefficient cose / cos@ is related to the 
"area defect’ in the slope distribution, observed at an angle of 
incidence @ [3]. The averaged value of (3) is conventionally computed 
for an isotropic normal distribution without any areal correction and 
for a certain (usually isotropic) spectral model. However, o” can be 
also evaluated analytically avoiding numerical computations. - 


2.1 Radar Gross-Section Evaluation. 


The actual long-wave slopes are small, [@]- |b] = 9 Cw << 1 | ang 
their distribution is approximately Gaussian [4]. Within a fairly good 
accuracy Oc u*> for slightly rough seas, i.@e., neglecting the slope 
skewness, we can replace the tangents of small angles in (3), (4) by 
their arguments. Further assuming Je]= 100 for cm and dm microwaves, 
the calculation of og reduces to a serigg expansion of |e} cot *e and 
Wek sind’, 2k code’) in w and je| , followed by their 
proatet avdéraging over slope distribution (4). Here, we can restrict 
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ava lt 
ourselves to the terms oC je] > MO . The result, accurate to a 


few per cent, is 
0,9) = Bk*ci+sin®6>"Cak sine>”*|i+c{e coseS*|*rcpoM . ced 


where @ is the azimuthal angle along the pricnipal axes, the typical 
roughness spectrum is given by the power-law [5] 


WCK.p) = BK “tcp , ca 


the modulation is expressed by 
pd A = an 2 -a 8 ~ 
pice.o, 8.0, 201m fr af mat m0 + m (+O Jo] 7M Ja andp 


are the average slopes, the tilt the variance 


Cm om, and mm m9 modul ation transfer functi 
Catt. Pia> depita SbiaSAtiy on @. jel **.e, and tte. 


2.2 Tilt and Variance Modulation 


The m.t.f.’s behave quite differently with respect to 6. For instance, 
at a saturated spectrum (6) of the actual azimuthal factor [3], assuming 


s=4, 140.8, |f°-f[¢1. the functions aC@> are similar to thoes 
in Fig.@. They are calculated for A = @.6 cm, i.e. for Cef.[6}5 


when @ € 7 the power scattered from an anisotropic sea is modulated in 
a complex way by all the slope components. In this case, the 
backscatter from isotropic waves [2], [7] is Slightly sensitive to 
slopes owing to the small m.t.f.’s, when @» 74°., however, 


j= | >> j= ,I- je sal >> 1... |. in... | We and the 
modulation arises primarily from the slopes of long waves oriented in 
the plane of incidence. At grazing angles, [=,|- , jm,,| are 
proportional to tan @, tan @ , respectively. Hence, the modulation by 
such slopes should predominate over other [7] scatter modulation types, 
which are somewhat limited. The number of S has a very slight bearing 
on the m.t.f’s and, as it changes from 3.5 to 4.5, the considered angles 
are shifted only by 1.5 degrees. These m.t.f’s occur over a wavelength 
range up toA= im, while the angles shift slightly towards @ = 90 


degrees. 
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Fig.2. M.t.f’s versus @ for A= 2.5 cn, 
|t’t [<0.6, je"-t ia 


2.3 Radar Measurement of Bragg Resonant Wave Spectra 


Another form of (6) yields 
_¥ 2... 8 Se ~-4 24 ~4 
o 60. go Ci +sin @> ja +c ecoadé) j Wk sind, di +0Cy, je] >.¢7> 


Thus, by measuring the backscattered power at various azimuths and 
angles of incidence and by averaging it over the longest sea wavelength, 
we can define, accurate to an order of long-wave slope variance, i.ée., 
to a few per cent, the directional spectrum wok, 9) of short resonant 
waves of i = aCk,sin@) . For isotropic long waves, (7) agrees with 
the numerical | findings [2]. It is also clear that (7) implies an 
accordance between the azimuthal anisotropy of @%¢g) and the angular 
spreading f€g) (6). For large averaging scales, all the conclusions 
drawn are also valid, sometimes at the expense of accuracy, for any kind 
of scatter modulation [7]. For directional spectrum measurements. 

SLARS and synthetic aperture radars (SARs) are most suitable. If their 
antenna vertical aperture admits angles 30° ¢ @< 90° , the 
measured spectral range can be as large as AK = @KC@ =30°> . 
high spatial resolution of SARs and SLARS permits, in principle, a more 
accurate spectral measurement than is achievable with conventional 
scatterometers. The azimuth can be varied by changing the aircraft 


heading. 
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2.4. Long-Wave Spectra From an Image Fourier Analysis 


The problem of how an image Fourier-transform can conform to the actual 
sea is still disputable. It can be resolved, on the basis of the 
foregoing examination, for angles of 70-80 degrees and moderately 
developed seas. A SLAR operates as a low-pass filter with a cut-off at 


= ani~*, SN being resolution. Then, in the first approximation, 
the 1mages recorded at ¥ = const correspond for 7>L_ to a spatial 
pattern of the smoothed slope a@%@ in the azimuthal plane. A 


Fourier transform of the image yields distribution of « yy 1n wave 
numbers and directions. 


Following [8], we examine the transformation of a coherent plane light 
wave passing through a radioimage modulated in amplitude as 


TR = Ti + a9 5), ce | 
where ® = {X,Y} is the coordinate in the image, 7, is its average 
transparency; ® @ is the mean value of m Ccos ‘nD ; h is the 


SLAR altitude; Y Stands for the range. The'light intensity in an off- 
axis region in the focal plane of a converging lens is proportional to 


1c = | [faci eon 18 ak |* . ce 


where K is the radius vector in the plane of transform. Measuring 
intensity (9), we determine the directional spectrum of long-wave slopes 


Ay- 
3. MODEL OF H-RADIOIMAGE FORMATION 


At oblique angles of incidence, the Bragg H-scatter may be a few dozens 
of decibels lower than that of specular reflection and the so-called 
’wedge’ [10] scatter. The two mechanisms occur only at the surface of 
breaking or nearbreaking waves. This point is corroborated by numerous 
experiments, in which a noticeable H-backscattering appears 
Simultaneously with the onset of white capping. A spike of the 
reflection, even averaged over radar resolution may prove an order of | 
magnitude more intense than the Bragg background and, hence, may produce 
a saturated density of recording films. Thus, each bright spot should 
be referred to a site at the sea that has undergone breaking during 
irradiation. To develop a simple model of the phenomenon, we assume 
that there exists a certain limiting amplitude of waves. If the 
amplitude is exceeded, waves are no longer stable and start to break 
{9}. Then, the speckle statistics in images reduces, in general, to the 
statistics of overthrowings of a two-dimensional random function above a 


specified height level. 
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3.1. Statistics of a Random Sea Surface 


To simplify the solution to the problem, we consider the Gaussian field 
of sea waves because their actual statistics differs just slightly from 
the normal case [4]. Most important for us is the mean number 

of overthrowings of the random function &¢x, yd above some appropriate 
level & = a per unit area. If the level a is fairly high, the number 
of overthrowings for an anisotropic random Gaussian field is determined 


by [11]. 


2.4728 ~ 
N <a m ~a? can, > athe oqXt -a" am.) ° ¢€10) 
tj 
wher @ at =f fren nH a a stands fcr spectral momenta. 





Take simple power-law spectrum (6) in the range K¢ K€ K, as a 
Saturated sea model, where K. corresponds to the high-energy component, 
is the high-frequency cut-off. This yields the momenta and (10) 


. 
chahges over to 


2.2 ‘ x a** 
N nf 95 2,8 pa 
oe Ny Sate Mean Gd Ke we EP 9, cas 
oO 


wher e@ ti Lad J fC pecce pe’ sin ede ‘ 


Hence, we see that the speckle density depends strongly on the typical 


width f - of an azimuthal spreading. Another important feature is the 
exponential dependence of N on the square of breaking wave steepness 
Cak 52 . Itis noteworthy that measurement of speckle density in 


images*is a very special procedure. a more common practice is 
photometering of an image density. A description of this operation 
shows that, at Gaussian seas, the mean signal obeys largely the same 


rules as does the number of spikes. 
3.2. Current-Induced Variations in the Spike Statistics 


If a surface wave runs to an inhomogeneous surface current U = UC x, >» 
the former and the latter undergo an energy exchange. Waves over an 
opposite-directed current are supplied by the energy from it, while, in 
the case of a divergent current, the reverse effect occurs. To analyze 
the spike statistics of the sea on currents which are due, for example, 
to internal waves, we use the conservation law of kinematic wave density 


9Coe + KID = o.. C12 


where o = (gkK*“* is the wave frequency at KU = KUcosep = 0. 
Using (129, we obtain variations in the wave number 


K’= 4K(1+¢1+4UfK/g coap)*)"*= xci-eufk/g coap)+ocufK7g> .¢130 
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For saturated spectra on the current, the high-energy and the cut-off 
components shift as K,-> Kj and K > K,. respectively. Substituting 
(13) into (10), we derive variations in W on the current y 

a 


AN 7 N.e< -Uccep [1 + oca"K> 3 ‘ c14) 


Thus, the speckle density in H -radicimages increases on currents moving 
against the wind and decreases in the opposite case. In addition, we 
evaluate the nonlinear effects in the speckle statistics and estimate 
essential corrections OaK 9 to the speckle density due to the second- 
order nonresonant nonlinearity of gravity waves. 


4. EXPERIMENTS 
4.1. Experimental Procedures 


Circular flights, which are conventional in azimuthal measurements [3], 
have been undertaken. A SLAR operated at A™®8. 28 cm was installed 
aboard the Antonov-24 aircraft of the ’Aeroflot’ Co. Due to the high 
spatial resolution of the instrument we have taken azimuthal scattering 
characteristics in detail. The flights were performed in the coastal 
areas of the Pacific near Kamchatka and of the Barents Sea near Kola 
Peninsula as well as in the Black Sea. The other flights were various 
’poxes’” and ’stars’, similar to that of Fig. 3. 
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Fig. 3. Aircraft 12-tack flight of the ’rotating box’-type 


4.2. Azimuthal Measurements 


With stable SLAR parameters. an average image density is proportional to 
the scattering cross section. Fig. 4 shows GO} for V-polarized (crosses) 
and H-polarized (dots) scatter versus. azimuth The Pine allege wae 
averaged over an area of S800 x a80 s* C3 x ~ % at @ » and 
then normalized to its maximum. The ratios of the Pete R _ <b D 
presented of Figs. 4 a,b to that in Fig.4 c are 0.1 : 1 and 0.72:1, 
respectively. The weather conditions shown in the caption were reported 
by a ship drifting 20 km away from the center of the experiment area. 


The azimuthal characteristics have conventionally an axis of symmetry 
coinciding with the wind direction, so Fig.4a suggests, for instance, 
the direction of a calm breeze. Figs. 4b and c differ very strongly in 
spite of close wind velocities. Thus unlike [3], this velocity is not a 
single parameter determining the backscatter. Its nature is quite 
different for different polarizations. Evidently, at V-polarization 
only Fig. 4¢ corresponds within a fairly good accuracy to the actual 
spreading fCgd for capillary-gravity waves of A®1.@em. According to the 
theory, this accuracy is provided by a deep azimuthal dependence of 

and by virtually symmetric scatter related to the wind direction. The 
other two V-characteristics reproduce f€@ , involving some distortions 
due to the skewness of long-wave slopes. While H-scattering is, in 
practice, isotropic under the calm (Fig. 4a), it is highly anisotropic 
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even at moderate winds (Figs. 46, c). Processing the azimuthal 
characteristics for the H-backscatter in the sector @ = 45-60 yields 
the second spectral momenta and estimates nonlinear skewness and 
curtosis of long-wave slopes. It is remarkable that the estimates are 
in a very good accordance with {12], while those made according to [4] 
under our experimental conditions, are considerably higher. 





| 
a 
U, 
s 
Fig.4. Azimuthal etarestersotite oe at 6 = 87°: 
crosses - ‘Ton Ge ee ~ hg py po 
Ca> calm, U = 0, handy, wind, Us 6 ms 
Cc> gusty wind, us x! S m= 


4.3. Fourier Analysis of Sea Radiocimages 


Two-dimensional Fourier transform of radioimages is observed in the 
focal plane of a converging lens in a coherent optical system [8]. v- 
images, taken over each of the 12 tacks of the flight and shown in Fig. 
3, were processed in the analyzer. A series of optical spectra is 
presented in Fig. 5, where each spectrum is placed at the azimuth of an 
original image. In the center, there is a spectrum of reference grating 


with A. = 50m. 


A feature of importance to us is the disappearance of waves at azimuths 
orthogonal to the wind, namely, at f = 150, 180, 330, and O degrees. 
Conversely, in directions almost collinear to the wind, the waves are 
Clearly pronounced. This confirms the theory asserting the predominant 


modulation of V- backscatter by long-wave slopes. 
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fig.6. Optical spectra of the 12 SLAR ‘rotating box’ tacks 








y, 


a a ee 


Fig.6. Resonant-wave spectrum KIO versus aziauth ¢ 
as determined from the 12 SLAR tacks 
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4.4 Spectra of Short-Wave Resonant Roughnesses 


The radar cross-section for each of the 1i2 azimuths was determined by 
measuring the image density versus slant range. Using a transform 
inverse to (7), the spectral density ®% @k, siné> was computed. 
Fig. 6. shows 12 of such spectra cut- off at @€ 70° to avoid strong 
modulations and covering the total azimuthal circle. At 42<¢ 6¢ 70 
all the spectra decline with increasing O to obey the universal power- 
law WhOec x~8-*t0-2 . Note that these angles cover the wavelengths 


1.2¢ %A¢ 1.7 cm. ; 


4.5. Nonlinear Wave Interaction as Observed by an Imaging Radar 


The typical Pacific radioimages were analyzed in a spectral domain. 

Part of the obtained spectra displayed several wave systems, some of 
which originated probably from the nonlinear interaction. Let us 
consider, for instance, the smoothed spectrum presented in Fig. 8., 
where one can see three wave systems. Two of them, making angles 55 and 
-25 degrees with the horizon, are found in other spectra. Hence, they 
are fundamentals. The third (additional) wave system makes an angle of 
around 30 degrees with the aircraft heading, and it is just interesting 
as a probable result of the nonlinear interaction between the 


fundamentals. 


a 
a 
Bb, 





Fig. 7 Smoothed radioimage spectrum and corrected wave vectors near the 
spatial resonance 


To verify the hypothesis, Fig. 8. plots the spatial resonance for the 
vectors of the additional (central) wave K. and of the fundamental 
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KK g, This plot 1s superimposed on the equal spectral 


(side) waves “sg 
density lines after corrections, to account for the relative motion of 
It results in the fact 


the platform referred to the sea peapserers.. 
that the frequency resonance xi73, Ki? oxi? is fulfilled within 
20. : 


a very : accuracy, while the considerable ae detuning 


ax - 
It a from theory that such nonsynchronous waves interact 
quasistationarily. Thus, the central wave component is expected to have 
some spatial beats. To check this claim, the spectrum of Fig. 8 was 
multiplied by a binary filter, which is totally transmitting over a 
sector 20 degrees orierted along . and nontransparent in other 
directions. The product was Fourier-transformed. The density of the 
reconstructed image is presented in Fig. 9. Here one can see the 
spatial beats of a period of 450-700m marked with broken lines. The 
theoretical spatial period was estimated to be about 500 m, which is 


close to the beats in the filtered-out image. 


- & —<, 1s measured from Fig. 10. 

















Fig. 8 Spatial Beats of the Central-wave Amplitude as a Result of the 
Image Processing 


S. CONCLUSIONS 

(1) The modulation of the sea radioimages of V-polarization, taken at 
oblique angles, is due to the long-wave slopes in the plane of antenna 
look. 

(ii) The speckles of H-polarized sea radioimages are explained by the 


statistics of specular reflections and wedge scatter from steep waves. 
The speckle structure depends strongly on the directional spectrum and 











(i1) The speckles of H-polarized sea radioimages are explained by the 
statistics of specular reflections and wedge scatter from steep waves. 
The speckle structure depends strongly on the directional spectrum and 
nonlinearity of the sea as well as on the inhomogeneous surface 


currents. 


(111) The developed models provide imaging radar methods to be used for 
experimental measurement of sea spectra of different scales and for 
observation of the effects of various atmospheric and oceanic processes. 
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[Text] ABSTRACT 


A lidar method is developed to measure the sea wave spectrum and other 
sea characteristics by using a narrow amplitude-modulated beam. Some 
techniques are suggested for determining the slope and curvature 
variance from a backscattered (reflected) power. The model field 
experiment was carried out at an oceanographic platform. The 
experimental data demonstrate that the introduced method is highly 


effective. 


1. INTRODUCTION 


To study spatial-time structures of sea waves is urgent for solution of 
various scientific and applied problems arising in the field of oceanic 
physics. Remote sensing methods seem to be very appropriate for this 
purpose. The most highly developed technique of oceanic remote sensing 
is associated with the use of air-space imaging radars. Now in 
radiococeanography the so-called synthetic aperture radars [1], recording 
the coherent backscatter from the sea, or the noncoherent side-looking 


radars [2] enjoy great favor. 


Since the early 1970s laser methods of sea probing have been also 
developing. They are superior to radiococeanography in several respects. 
First, lidars can provide higher spatial-time resolution, which is 
Suitable for measuring small-scale waves. For instance, a sharp-focused 
laser beam permits us to measure waves of a horizontal scale up to l mn, 
while the microwave radar resolution is limited just to several meters. 


The duration of a radar pulse is typically from 1077 to 10% {3}; 
accordingly, the spatial resolution is 15215 Qn. . In the optical 
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range one can produce pulses as short as 5.10713, [4] enabling a 
resolution of 0.08 mm to be obtained. We observe the principal 
difference between radar and lidar probings in the reflection of signals 
by the sea. The radar wave length exceeds often those of the capillary 
waves, and therefore, various diffraction effects are of importance 
here. Meanwhile the laser wave length is considerably shorter than any 
surface one, and the mechanism of laser reflection is reasonably 
explained in terms of geometric optics. This aspect gives rise to much 
Simpler relations between the statistics of the lidar return and that of 
the sea surface, compared to radars. One of the important parameters of 
the lidar return is its average power. As shown theoretically [5], [6], 
lidar measurement of an average backscattered power at different 
uncidence angles establishes variances be and & y Of the downwind and 
crosswind components of the surface slope. 


Presently there are two principal ways of oceanic lidar remote sensing, 
namely, pulse ranging and phase profilometry. The former [7] involves 
determination of the time delay between sounding and reflected pulses as 
well as analysis of their shapes and durations. Depending on the 
experimental facilities and position (for example, on the beam 
divergence and vertical or slant incidence), we can estimate the 
significant wave height, the slope probability density, or measure 
roughly the sea surface profile. Laser measurement of the profile can 
be also carried out with airborne phase altimeters. Continuous flight 
altitude measurements permit reconstruction of the surface profile along 


the flight ground line [8]. 


An analysis of the state of the art in these methods [9] reveals that 
the two techniques are effectively applicable to investigation of sea 
waves as long as 10mm and more. But in many field experiments it is 
important to study both large-scale and comparatively short waves. A 
good example of this situation is the interaction between internal and 
surface waves, where most changes are observed just in the short surface 
wave spectrum. Further development of laser investigations of cm and dm 
sea waves requires a finer spatial resolution. In the field of pulse 
ranging one can decrease the sounding pulse duration. The use of 
subnanosecond laser techniques calls for faster and, hence, more 


sophisticated receivers. 
2. NARROW-BEAM LIDAR RETURN FROM THE SEA SURFACE 


Let a laser beam of an intensity I(r,0), where r is the distance from 
the beam axis and O is the azimuthal angle, be normally incident upon 
the sea surface. In the horizontal plane at above-surface elevation H, 
there 1s placed a receiving aperture D (Fig. 1). The distance between 
the axes of the beam and the aperture equals 1. The characteristic 
Giameter of a sounding beam is equal to d. Then introduce the 
coordinates (x,y,z) with their origin at the receiver axis. The 
coordinates x,y lie in the horizontal plane of an undisturbed surface, 
and z is directed upwards. Any point of the surface is defined by its 


elevation g¢x,y) and the slope vector. 
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P Cxry2 = CE = dtvdx, g © aedy> . : c1> 


where K = (K1,K2) is the wavevector of surface random disturbance. 


the momenta of the directional surface spectrum Eck, .k,2 can be 
Witten in the form : 


Jj 
a,° ff Ck ok.> kd dk dk. 4.J00,4.8..., C2 
here Rack ok, > is the wavevecter of surface random disturbance. 


To evaluate the statistics of the return, it is necessary to assume the 
following: (i) the surface yr nege are. obviously gentle, %4™ << 3; 


(11) the altitude is great, iD><e oe ; a the sea returiv Ss 
homogenous within the beam spot, «<< w ; (iv) the transmitter 
and the receiver are nearly coaxial 1 44 sg 4e8 >; (v) the beam is 
narrow, nS << 1, where n is the density of = feouhh-’ points at surface 
Slopes (#2206. , and S$ stands for the spot square; and (vi) the 
surface is ¥ii1uminated from the ’far’ range, X<jRI>A1 << n “*? where 
<{R|> is the average modulus of the radius of surface curvature. Under 
all these conditions the lidar return from the sea surface is provided 
by the specular points within the lidar spot. The power returned from a 


Single specular point is 
P_= 0. 0688" ja) “cD1D"p Icr.6d. ¢3D 


Here 2” 0.02 is the sea reflectivity at the normal incidence, (9) ls 
the total curvature at the specular point. Let us normalize the point 
return as 


{ 1, if the specular point is placed inside the beam spot, 
* 


0, if it is outside the spot. 


Assuming that the sea surface is a random Gaussian field, one can obtain 
some statistical quantities: the average power 


* 
Ps 


<P> 2s a6 $ C4) 


the power variance 


a<jat) * 
“ I 


and the average normalized return 


(p> = < [apes = ns. ; ce 





SS t%cr,@) r dr de ; : rq .») 
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Fig. 1. Horizontal cross-section of the probing geometry. 


Here ae . ~ 2** the radiated power, 6« ~pel 
sad nt, - m* is the surface invariant, ch 
<b « on mae wave ea and divergence. . 





Fig. y: Discretization of the process S (t) by a random pulse sequence 


4. 

Thus, the measurements of ‘P,>. <B>, <P> make feasible to 
determine the sea cross-section @& , the invariant , and the mean 
moduli of the curvature and of its reciprocal: 


o = <P> 4m “CHDO®, aX? = avenP, ; C7) ,C@. 
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a <P <p*) p* 
(jap = 5 » <fa°* pe #2 
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a <P> 


respectively. 
3. PHASE PROFILOMETRY WITH A NARROW BEAM 


In the case of phase profilometry, the laser beam is amplitude-modulated 
with some frequency #, . The return is formed by specular points 
mentioned above. To monitor the profile continuously, the lidar surface 
spot should be widened, thereby degenerating spatial resolution and 
averaging surface elevation over the spot. Oceanographers are primarily 
concerned in the sea spectrum or other statistical quantities rather 
than in isolated sea elevation records. The problem is, therefore, to 
reconstruct the statistics of an unknown random sea from an experimental 
discreet sample corresponding to the specular points. 


Let - (t) be a stationary random function of f e surface elevation in 
time (or in space). The discretization of €*% (cf. Fig. 2) by a 
random sequence y(t) produces a new random n(t). If n(t) is of th 
Poisson distribution with a mean pulse period Te, the spectrum S, ww of 


n(t) depends on the spectrum S,(W) of £ct> as (10): 





4 er* cw’>*s Cw’d 
sco are () ° fn 2.3 
ite Tf . 4 +<Cw’>" 2, 


It is evident from (11) that the discretization of a random process by a 
pulse Poisson sequence is equivalent to incorporation of a high-pass 
filter of a cut-off frequency G), = ~£ , which is connected in 
parallel with the source of ‘white’ Poise of the spectrum S. 
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ded nok 
Sea eurface 
Fig. 3. Experimental set-up. 
For an initially limited sea spectrum, where 36-0 © at W>%&%, 


we reconstruct %S,Cw by the following: (i) Sg©@,9 is evaluated at’ 
47m ; (ii) the ‘white’ noise spectrumbe $cwcie we > 


is found; and (iii) the initial spectrum (wo . SCw¢1+ ws 2- § 

is reconstructed. When probing is carried ouc from a stable platform, 
we can show that the sea specular glitter is described rather well by 
the Poisson statistics, that allows us to apply the suggested algorithm 
for reconstructing the sea spectrum. 


4. EXPERIMENTAL SYSTEM AND SITE 


To study the statistics of the laser return from the sea surface, we 
have designed a lidar shown schematically in Fig. 3. The source was a 
15mw CW He-Ne red (633nm) laser LG-79/1. The sounding beam was 
amplitude-modulated by an electrooptical modulator. The modulation 
frequency (40 MHz) corresponded to the 2n-uncertainty range of 3.75m. 
The modulated beam was focused by a system of lenses, Ll, L2, and 
directed vertically downwards through prisms, Pl, P2, and an output 
mirror, Ml. 


The return reflected by the sea surface and transmitted through mirrors. 
M1, M2, was focused by a lens, L3, onto a diaphragm, di, inserted in the 
image plane to limit, for the sake of better sensivity, the receiver 
field of sight. The input aperture, D, was defined by the L3 diameter 
of 4.7cm, while the dl diaphragm diameter was from 0.04 to 0.3mm. So 
the receiver visual angle was 10’, i.@e., the linear dimension of the sea 
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Site visible through the diaphragm was from 0.7 to 5 cm. The return 
beam cut-off by the field diaphragm was made to pass through a 
collimator, L4, and an interference filter, IF, then it was registered 
by a photomultiplier, PM. The radiopulses of frequency 40 MHz from the 
PM output were fed to a detector of envelope and phase. The analog 
envelope and phase signals from the detector were recorded by a Bruel 
and Kjaer 7005S tape recorder. The signals were also processed using a 
DISA 56 Nli-type mean-value meter, a Bruel and Kjaer 7005S spectrum 
analyzer, and a CLANP 0270 computer. 


The lidar experiments were conducted during May-June 1987 at an 
oceanographic platform located in the Black Sea SOOm off the Crimea 
coast. Under various hydrometeorological conditions, we carried out 70 
series of 20-minute measurements. The distance between receiving 
telescope and the sea was 16m. Some of the experimental series were 
accompanied by in situ measurements with a standard wire waverecorder. 


5. DATA ANALYSIS AND DISCUSSION 


After processing the geturn eyelope in accordance ath (7)-(10), we 
evaluated o,A,* 6A,A,), = <jaj>*”*, and R =< jal” 

then anelysed” their “dependence on the wind velocity, We Spectral 
processing of the phase signal and its correction (11) allowed us to 
determine the sea spectra under various winds. The results were checked 
partly by wire waverecords and compared to the data reported elsewhere. 
The sea cross-section versus the wind velocity for different sounding 
spots 1s plotted in Fig. 4. The value of @ for vertical incidence 
falls rapidly with the wind strength and it varies from 0.5 to 0.0l, 
when W changes from 1 to 12 m/s. For comparison, we show here the data 
[S] obtained with an airborne green (530 nm) Nd: TAG lidar. 
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Fig. 7 R versus WwW. 


Assuming the Gaussian sea statistics, we can readily evaluate the 
product ft, fy of slope variances from OO. Fig. ge Pketones the 
relationsnip between # and W, as well as a*’"cW, . which is 
derived from semiempirical {11}) relations (xc. a*cw and 
evaluated from a model [12] spectrum. As seen, out data aré ina 
satisfactory agreement with the results cited. Figs. 6,7 show the 
experimental graphs aCW> and RCW. respectively. Fig. 6 also gives 
the semiempirical [13] relationship @CW . The errors present in Figs. 
4 through 7 are related to the imperfect lidar graduation. The 
variations in experimental data may be partly attributed to several 
uncontrollable factors, i.@., to duration and length of wind fetch, wind 
turbidity, differing temperatures of the air and sea, etc. Notice also 
that, at the widest sounding spot of 5 cm dia, correct measurements of X 
and R are possible only under rather slight (W < 4m/s) winds because 
under stronger winds the specular point density increases sharply and 
the condition (v) of narrow spot (Sect. 2) is no longer fulfilled. 

Figs. 8a,b exemplify different sea spectra, where we demonstrate: (i) 
the spectrum of an initial phase signal; (11) the phase spectrum 
corrected in accordance with (11); and (iii) the waverecorder surface 
elevation spectra; respectively. A comparison of all the spectra 
exhibits good coincidence of the rough lidar and recorder data up to 
frequency of 1,25 Hz. After correction the upper limit of the 
coincidence shifts to around 2.5 Hz. 
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Fig. 8 Examples of the spectrum reconstruction. 


Thus, the suggested method enables a noticeable progress to be made in 
lidar remote measurements of the high-frequency sea waves and of the 
Slope variance and the curvature of a disturbed sea. 
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[Text] 1. SATELLITE PROGRAMS: 


In India, the National Remote Sensing Agency (NRSA) of the Department of 
Space, Government of India had been given the responsibility to receive, 
process and disseminate the satellite data products to various users in 
the country. Accordingly, NRSA set up the satellite data reception and 
processing facilities initially for the MSS data from Landsat-2 and 3 in 
1979. Subsequently, the station was upgraded to handle data from 
Landsat-4 and 5. When the more powerful sensor, viz., Thematic Mapper 
(T™) was introduced and the TM data was available to foreign ground 
stations (on Landsat-5), NRSA set up the facility to receive and process 
this data. The system integration, hardware and software development 
for T data processing were done inhouse. The supply of TM data 
products to the user community commenced during end i985. 


During 1979-85 the growth in the use of remote sensing techniques had 
been substantial and this led to increased expectations regarding the 
potential uses of remote sensing. The SPOT data reception and 
processing facilities were established and declared operational from May 


1987. 


Since March 1988 the Earth Station is receiving the LISS-I & LISS-II A&B 
data from the Indian Remote Sensing Satellite. 


Since 1979, NRSA has also been receiving the data from meteorological 
satellite series (NOAA 6 to 11). The comparison of satellite/sensors 


volume of data handled is given in Table 1. 


Being an active sensor which can ’see through’ clouds, Synthetic 
Aperture Radar (SAR) has been an all weather sensor with invaluable 
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applications in tropical countries such as India where perennial cloud 
cover in many parts of the country precludes the utility of optical/IR 
sensors. Owing to these advantages India also is planning at present 
for its own space borne microwave remote sensing mission (with the 
Synthetic Aperture Radar (SAR) forming an important component by 1993- 
94). 


2. DATA PRODUCTS GENERATION AND DEMAND: 


When the distribution of data products began in early 1980, the remote 
sensing technology in India was in its initial promotional phase. 
Satellite data products generation and supply were a few hundreds in 
1980. 


In the last few years, the number of users requiring the data products 
increased at a fast pace. Initially, the major users were from within 
the Department of Space, which is the nodal agency in India for 
developing and promoting remote sensing technology. Gradually, the 
other Central Government agencies also started using data products. 
During this period, to bring out the user awareness, seminars were 
organised, presentation made, workshops held, courses conducted, pilot 
projects completed, results presented to decision makers, etc. These 
promotional efforts have had a telling effect on the growth of demand. 
Now Central and State Government agencies are the major users. Many of 
the universities have started teaching remote sensing and have post- 
graduate diploma courses on the subject. The growth in the number of 
users also has been substantial. 


With this large growth in users, the demand for data products also 
increased substantially. While in 1980, about 30-40 products requests 
per month were processed, the number increased to more than 900 products 
per month in 1987. The demand in the last few years and the 
corresponding supply is shown in Fig. 1. 


3. SATELLITE DATA PROCESSING: 


With such a multiple satellite data processing program, hundreds of 
gigabytes of data are to be handled every day. The processing steps 
involved in handling satellite data are reception, archival, storage of 
"information on data’, and its retrieval, processing of the data to 
convert it to a usable form using computers, photo processing 
equipments; distribution, application/analysis of such data, maintenance 
of such facilities, etc. Allied steps are training of personnel, 
conception and execution of national level projects to use the data. 
(Fig. 2 and 3) 











TABLE 1. COMPARISON OF SATELLITE/SENSORS VOLUMES 
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LISS-1i 36.25 87 . 
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To handle the large increase in products supply, a multi billion dollar 
program is under execution at NRSA, which envisages going in for newer, 
technologically advanced, faster systems. Firstly, large volume data 

calls for a well thought out archival policy where only useful data are 


archived. 


Next, an integrated information management system in the form of an 
electronic data base is planned to be implemented which serves as an 
external user interface while at the same time controlling the internal 
data flow and production. This system is designed such that it serves 
as the master assigning jobs to all the slave work centres. 


The data processing systems so far have been using the ‘Conventional’ 
approach where in the data is first transferred from high density tapes 
to disk and then the data is corrected in Array Processors. This poses 
limitations in thruput and utilizes system resources less efficiently. 
To overcome this, using an I/O computer connected to the conventional 
processor, the input/output speed is enhanced and more efficient 
utilization of resources achieved. In yet another innovation, a Master 
Flexible System (MFS) is being designed and implemented, which shall 
have the capability to process all types of sensor data. This is 
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essential in our environment, where the demand for different types of 
products from different sensors may vary over time and such a MFS can 
act as a ‘load’ balancer. 
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2] Demand SQ Supply 
Fig. 1 Demand and Supply of Data Products 


Large percentage of the data products required are photo products; hence 
the program envisages going for large volume production using roll film 
recorders, automatic photo-processing machines and the building of a 
modern, controlled, dust free photo processing laboratory. 


3.1 POLICY FOR THE ARCHIVAL OF DATA: 
The archival of satellite data has been in the form of raw data on High 
Density Tapes, quick look Video Cassettes, Quick Look Films, User 


Computer Compatible Tapes (CCTs), 240mm Negative/Positive 8/W or Color 
master films. 


The High Density Tapes (HDTs) contain data from Landsat MSS, RBV, TM, 
SPOT, Metsat as well as IRS data. 
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Fig. 2 Master Data Products Generation: Data Processing, Filming and 


Photoprocessing 


Enormous amount of data had been archived in the last 5 years. (Table 2 
and 3) The policy so far, had been to acquire and archive as much data 
as possible. With the increasing volume of data, it became apparent 
that a suitable archival policy had to be worked out and implemented. 
Accordingly, after wide consultations, a draft archival policy has been 


worked out, which when implemented will 


--ensure that an archive of representative and useful satellite imagery 
of India exists, is maintained and is accessible; 


--ensure that selected scenes from the existing data set are not lost; 


--reduce the tape storage and tape consumption to manageable 
proportions. 
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Fig. 3 Data Products Generation Flow 


The three tier criteria for archival is as follows: 


Tier 1: Within a week of acquisition, the HDTs having technical 
reproduction problems will be identified. These will not be considered 


at all for archival. 
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Tier 2: After a year of acquisition, data to be retained will be 
subject to the criteria that only when 5 or more scenes in a pass have 
less than 40 percent cloud cover, the data of that pass shall be 
retained in the archives. Before destroying the weeded out HDTs 
adequate notice will be given to the user community to indicate whether 
they have any need for the data identified to be removed. During such 
weeding process it will also be ensured that sufficient data for Target 


Sites are available. 


Tier 3: After 5S years of acquisition the following data only will be 
retained: 


i. One pre-monsoon coverage of India 

1i. One post-monsoon coverage of India 

1ii. Disaster affected areas 

iv. Data at the required frequency for Target Sites 


Table 2. ARCHIVED LANDSAT QL AND CCT DATA 








Satellite Archive Medisa Date Range of Approz.Mo. of 
heme @ Date Format Archive Data Frames 

Lenest 465 70 an Sept.1962 to 33,500 frames 
os Film rolls Mar. 1965 

CLendent 408 Quick Look April 1965 to 5,900 ecenes 
—- Video cassettes August 1965 

Lemdeat 4a4 Quick Look Feb.1986 to 20,800 scenes 
cad Video cassettes Oct. 1986 

Lendest 2,3 160@ bpd Howv.1979 to 650 ecenes 
- Digital Tape Dec. 1962 





3.2 INTEGRATED INFORMATION MANAGEMENT SYSTEM: 


NRSA Data Centre is the focal point of data products generation 
activity. It is the interface between the user agencies and NRSA. It 
handles the data products requirements of all the users both internal 
and external. The primary function of NDC are User assistance, data 
archival, User order processing and data dissemination in response to 
user requests. The user interface with NDC typically includes both 
general enquiries about the data availability and ordering procedures a 
well as requests for data covering specific areas. NDC also provides t 
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any user, on request, the information on the status of his order or 
account. It enables the visitors to view Quick Look data and make use 
of reference facilities such as catalogues, orbital calendars, etc. 


Table 3. ARCHIVED LANDSAT HDDTs 








Satellite Archived Medium Date Range of © Approx.lo. of 
6 Sensor & Data Format Archive Data Frames 

cadens 45 HDDT Aug. '82-Oct.'86 1,50,40@ scenes 
lendeat 5 Ty HDDT June'64-Nov. '84 6,250 scenes 
landest 5 Ty HDDT Mar. '85-Nov.85 15,600 scenes 
landest 5 Ty HDOT Jan. '66-Oct. '86 17,000 scenes 





The volume of information being handled is enormous. Currently, NDC is 
operating in a semi-automatic, semi-annual mode. Because of large 
growth 1n demand, it is planned that the entire activities of NDC is 
computerised and handled through the data base. This can be achieved 
through a multi satellite, Integrated Information Management System 
(IIMS). IIMS is a system to store, manage and retrieve information of 
several Users of the system as a resource for their mutual benefit and 
for several purposes. IIMS will also be used as a production 
controlling, monitoring system work centre wise and also create job 
orders for various work centres. The system will be based around two 
VAX processors. NDC data base operations are very critical. The 
information to be dealt with is voluminous and prompt information 
retrieval is desirable to answer the Users queries within reasonable 
time. To meet this criticality, it is necessary that system should have 
no ‘down time’ virtually. To ensure this, the system has a redundant 
processor, redundant intelligent I/O server and sufficient number of 
disk drives. Data processing computer system will be connected to this 
system through Eathernet Local Area Network. The Gateway will connect 
_Decnet nodes on the Eathernet network to Remote Users for information 
retrieval. This system will store and retrieve data kept in by several 
Users such as DATA RECEPTION, DATA PROCESSING UNITS, PHOTOLAB and DATA 
CENTRE. Using the Eathernet link; information such as type and details 
of the products to be generated, HDT No., priority of products, products 
generated, quality of products, acquisition and data processing 
information, system status, etc., can be exchanged between this system 
and the various Data Processing Computers, Photoprocessing, Filming 
System, Quality Control, Data Acquisition, Data Archival, etc. 


4. CONCLUSION: 


With the completion of the activities relating to the augmentation 
programme the Data Products supply is expected to be streamlined and 
will have a Turn Around Time (TAT) of within 7 days for 50 percent of 
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the products and within 21 davs for the balance. It will also be 
possible to cater to increase in user demand for any one type of sensor 
products with the corresponding decrease in any other sensor products. 
[t will also ensure guaranteed supply of products in case of emergency 
and the data products quality 1s also to be better. While, going 
through this programme it will also be ensured that future expansion can 


be easily done. 


Till about a couple of years back, the facility available for digital 
analysis of the Remote Sensing data was practically not existent in the 
country. To overcome this situation, S Regional Remote Sensing Service 
Centres along with two Associate Centres which can perform high thruput 
Gigital analysis have been established. It 1s likely that in the next 
few years the digital analysis of satellite data will increase manifold 
and Remote Sensing will be used in an operational mode regularly. 
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Real Time SAR Processing Techniques 


43070008b Kyoto Selections from INTERNATIONAL ARCHIVES OF PHOTOGRAMMETRY 
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[Article by Roland Schotter, DORNIER SYSTEM GmbH, P.0.8. 1360 D-7990 
Friedricnschafen, Federal Republic of Germany] 


[Text] Abstract 


Spaceborne Synthetic Aperture Radar systems provide raw radar data 
information at high data rates of 10 Mwords per second. Real time SAR 
processors, therefore, must perform more than 1 giga-operations 
(multiplications, additions) per second in order to generate images from 
the raw data. Since conventional computer systems are not able to cope 
with these requirements DORNIER has developed real time SAR processors 
on the basis of a modular pipeline concept. 


The processing pipeline is built up to standardized hardware modules 
which are required for digital signal processing of two-dimensional 
data. These modules show identical electrical and mechanical interfaces 
so that each hardware module can be used in any place of the pipeline. 


The paper presents the basic principles of the real time modular 
Pipeline concept and its processing rate requirements for real time SAR 
processing applications. The implementation of some of the most 
important modules like Fast Fourier Transformation, correlation, 
interpolation and data memory is described. 


Finally, a short outlook on further applications of the pipeline 
processing concept is given. 
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Introduction 


Synthetic Aperture Radar images (SAR images) of the earth are an 
important tool for many scientific ranges and application related 
problems. Due to the fact that radar waves are able to penetrate even 
though strong clouds at very low attenuation, SAR images are available 


Gay and night independent of weather conditions. 


SAR data can be applied in a multiple way. In oceanography theories on 
the origin and extension of ocean waves can be investigated. Limits 
between iced and iceless areas can be detected by means of SAR images as 
well. This facilitates for example finding optimum ship navigation 
courses and the operation of ocean platforms in polar areas. Since 
icebergs are detectable in SAR images, SAR information is also 
predestinated to be applied in iceberg warning. 


The scattering coefficient for radar waves mainly depends on the water 
content of the illuminated objects. Therefore, conclusions on the 
growth conditions of vast agricultural areas can be derived from the 
radar images (wheat and tobacco fields, detection of forest damages). 
Geological structures are often clearly visible in SAR images so that 
radar images can be used in this field, too. This applies mainly for 
the exploration of hardly accessible areas. 


Further applications can be found in overall ship control (200 miles 
zones, environmental oil pollution monitoring on sea), in archeology and 


glacier research. 


Most of the above mentioned applications require real-time conversion cf 
SAR raw radar data into images for its efficient use. For this reason 
DORNIER has implemented a digital real-time SAR processor for earth 
exploration in the framework of a BMFT/DFVLR contract. The missions ERS 
1 (C-band) and X-SAR (X-band) are intended to be the first applications 


of the developed SAR-processor. 


The SAR System 


The complete SAR system consists of a sidelooking radar and a SAR 
processor which converts raw radar data into SAR images. The basic 
geometry of the SAR is shown in Fig. 1. The radar antenna travelling at 
uniform speed illuminates a certain swath on ground which returns radar 
echos to the antenna according to its reflectivity. These signals are 
converted in several steps from the carrier frequency band into the 
baseband by the receiver. The SAR processor transforms the raw radar 
data into images which have to be analysed according to the relevant 
application. The system control data (orbit altitude and velocity, 
disturbances of the antenna trajectory, etc.) are also used by the 
processor during the generation of SAR images. 


Fig. 2 shows a block diagram of the most important components of a SAR 
system. These components are basically used both for airborne as well 
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as for onboard SAR systems. However, there are some differences 
concerning the processing of the raw radar data which can either take 
place onboard the carrying aircraft or ina ground station. In the 
latter case the radar data have to be transmitted by means of a 
telemetry system, as it will be done for ERS 1 and X-SAR. 


The Real Time SAR Processor 


The signal processing task of a real-time SAR processor is characterised 
both by complex algorithms to be implemented as well as by high input 
Gata rates. Even the most advanced computer systems are presently not 
able to solve this task by means of software programs. 


Fig. 3 shows the block diagram of a SAR processor. Ina first step the 
Gemodulated radar signals are converted into digital information. 
Usually, real and imaginery I/@ data are used in the baseband for the 
representation of the amplitude and phase of the radar signal. 
Afterwards, the range compression of the raw data is performed. The 
linear range migration effect of the signal energy has to be compensated 
before the azimuth compression and look summation can take place. The 
“corner turn memory” converts the processing direction from range into 
azimuth. The azimuth processing is succeeded by the “post-processing.” 
In this unit the processing is changed to range direction, again. 
Additionally, the correct pixel mapping of the compressed radar data is 


done in this unit. 


The “Doppler Estimator” and “FM Estimator” units derive the 
characteristics of the azimuth reference function either from the radar 
data itselves or from the orbit and attitude data of the antenna. 


Concerning processing rate requirements in terms of millions of 
arithmetic operations per seconds (MOPS) under real time conditions, 
range and azimuth compression (see Fig. 3) are the most demanding units 
within the processor. Given the ERS 1 SAR parameters, roughly 900 MOPS 
are necessary for range compression in the frequency domain, whereas 350 
MOPS must be performed for azimuth compression using the SPECAN 


algorithm. 


Processor Implementation 


A pipeline architecture as shown in Fig. 4 is most suitable for realtime 
digital SAR processing. The raw radar data enter the pipeline at the 
left hand side and are transformed--step by step--into the focussed 
image. This structure can be well compared to a manufacturing line in 
which the final product is assembled at different stations working at 
the same time on products which show different degrees of completeness. 


The processing pipeline consists of digital hardware modules like 
Fourier transformation, interpolation, lowpass filter, twodimensional 
Gata memory, complex multiplication etc., the sequence and number of 
which are determined by the algorithm to be implemented by the 
processor. Several identical pipelines can be used in parallel if the 
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requirements in terms of processing speed cannot be met by a single 
channel. 


The processing parameters of each pipeline module are provided by the 
assigned pipeline controller. Tnese controllers compute local module 
control parameters from some global control information which is 
received and distributed by the supervisor controller. 


The processor nardware was implemented under the aspect that each 
specifically designed hardware module must include enough commonallity 
and flexibility to make it a part of a more general signal processing 
hardware family. This concept allows the functional modules to be used 
at several places of the processor or even supports the implementation 
of different image processing tasks (e.9g., other SAR algorithms, pattern 
recognition, etc.) at very low development cost by using the same 
hardware modules in different arrangements. 


All hardware modules incorporate identical electrical and mechanical 
interfaces in order to guarantee full modularity within the pipeline. 
Power saving CMOS technology has been chosen as far as possible to cover 
applications which suffer from lack of electrical energy (e.g., mobile 
processing in aircrafts, etc.). 16 bit word length of both the real and 
imaginary data channel is offered by all modules allowing the hardware 
to handle the full dynamic range of present SAR systems. 


Some features of the more important modules of the pipeline SAR 
processor are described in detail in the subsections below. 


High speed storage unit: 


Intermediate storage of data frames is a most important function of any 
kind of image processing. In a SAR processor these memories are used 
Guring corner turning, look summation, recorner turning and test pattern 
generation. In order to cover the above applications a universal memory 
board has been developed which allows for writing and reading line based 
data in both azimuth and range direction. The length of a line can be 
programmed within the limit of 256k words. The memory boards can easily 
be cascaded if more than 256k words of memory size 1s required for a 


specific application. 


Dynamic NMOS memory chips and Advanced Low Power TTL control circuits 
are used for the board the mechanical size of which conforms to double 
eurocards. A read/write speed of up to 3 Mwords/sec could be achieved 
at a power consumption of 3.2 Watt only. Refresh circuits for the 
dynamic NMOS chips are included on the board. 


A complete memory system (e.9g., corner turn memory) for data frames of 
1024 x 256 words is built up of three memory boards and a simple 
controller board which provides the correct input/output interface and 
some trigger signals, which tell the memory boards when to read or write 


Gata frames. 
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Finite impulse response (FIR) filter: 


FIR filters are widely used for low pass filtering, bandpass filters in 
combination with a complex premultiplication, digital interpolation and 
correlation. In a SAR processor all functions mentioned above are 
required at-different places of the pipeline. 


A FIR filter offering programmable length of up to 256 samples has been 
developec using CMOS multiplier/accumulator circuits and CMOS/TTL 
control logic. The board 1s able to perform all functions described 
above. In case of low pass filtering the output data can be subsampled 
by noninteger factors. The filter coefficients can be downloaded from 
the assigned pipeline controller. Up to 256 different filter sets can 
be stored at a time. 


The maximum input data rate is given by the internal 20 MHz computation 
clock and the selected filter length. 


Fast Fourier transformation (FFT) module: 


A radix 2 decimation in time algorithm has been chosen for the high 
speed FFT module. This unit shows the following features: 


--complex input/output 16 bit I / 16 bit @ 
--programmable FFT length (2 up to 4096) 

--FFT or inverse FFT selectable 

--data weighting function loadable from external source 
--programmable gain 

--complex butterfly operating at 15 MHz 

--full CMOS technology 


A 4k FFT can be performed at a rate of 600 Hz whereas a shorter length 
of 32 runs at almost 140 kHz. 


There 1s a wide range of applications for this module. Within the SAR 
processor both range and azimuth compression rely on this function in 
case of longer replica lengths of the respective reference function. 
For shorter replica lengths also a time correlation solution might be 
selected as the proper hardware choice. 


A prototype real time SAR processor for 1024 range samples of range 
compressed SAR data has been built on the basis of the pipeline 
architecture as described above. This processor is able to handle both 
ERS 1 and X-SAR data Dy simply initializing the processor according to 
the specific requirements of the actual SAR. The power consumption of 
the complete processor amounts to 120 watt thus demonstrating the 











efficiency of the hardware concept also in terms of electrical power 
requirements. 


















Further SAR Processor Activities 


The developed SAR processor has preven the feasibility of digital real- 
time SAR signal processing. An operational SAR processor including 
range compression and the Mean Doppler and FM Estimator units 1¢ 
presently being developed for the ground station of the X-SAR project. 


A Quicklook processor for airborne SAR processing (DOSAR) 1s also being 
designed and developed within 1988. This processor will include 
additional functions for motion compensation. 
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Fig. 2. Biock Diagram of SAR System 




















Fig. 3. SAR Processor Functions 


81 














Fig. 4. Block Diagram of Pipeline Architecture 
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17-22. Akasaka, Minato-ku, Tokyo Japan, Commission Number VII] 


{Text} 1. Introduction 


Several sets of OPS (Optical Sensor) bands to be mounted on Japanese 
ERS-1 (Earth Resources Satellite-1) were simulated and evaluated. 
Simulated data were produced by the use of ASR (Airborne 
SpectroRadiometer) of GER (Geophysical Environmental Research, Inc., 
USA). Data acquisition was conducted in Death Valley, California and in 
Comb Ridge, Utah. Evaluation of these simulated data contributes to the 
finalization of OPS specification. In this study, ground truth using 
field spectrometer IRIS was also conducted to check ASR data. 


2. Acquisition of ASR data and ground truth by IRIS 


(1) Description of ASR 


ASR that we used in this study is capable to measure radiance, along the 
track line (Image data is not available for this sensor). Major 
parameters of ASR are listed in Tab. l. 


(2) Description of data acquisition 


Aircraft is AZTEC-E, and the two test sites, Death Valley and Comb 
Ridge, were selected, considering the degree of vegetation cover and the 
types of soil/rock. Major parameters of data acquisition are listed in 
Tab. 2. Airborne data acquisition was conducted from May 24 to May 27 
for Death Valley, on May 31 and June list for Comb Ridge, both in 1986. 
During the flight radiance, data on the ground surface were measured bv 
the use of IRIS with equal time interval in order to correct the 
incident solar energy. After preliminary check of ASR data, ground 
truth using IRIS was conducted for the targets which are located at the 
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cross points of two track lines and also for the targets including 
interesting spectral curves. 


(3) Analysis of rock sample 


Some rock samples, against which we measured on site-spectra by IRIS, 
were sent back to Japan. Samples including clay minerals or carbonates 
rock were analysed by XRD; carbonate and evaporites, by X-ray 
fluorescence analysis; and volcanics, micro scope using thin section. 


(4) Comparison between ASR data and field spectra by IRIS 


For ASR data, only radiance data were available. so, all the data were 
Givided by the spectra that 1s supposed to have no spectral absorption. 
They will be called pseudo-reflectance data hereafter, and they are 
compared with a set of IRIS spectra for the corresponding site. One 
example 1s shown in Fig. 1. This site is basically an outcrop of 
limestone partially covered by low vegetation. However, absorption at 
2.35 um 1S well expressed both on ASR and IRIS spectra. 


(S) Comparison between IRIS spectra and analysis of rock samples 


Example for limestone is shown in Fig. 2. As shown in this figure, IRIS 
spectra, XRD and X-ray fluorescence analysis indicate the existence of 
calcite. In almost all the results, IRIS spectra and analysis of rock 


samples are consistent. 
3. Simulation of ERS-1/0PS 
(1) Flow chart for simulation 


Flow chart 1s shown in Fig. 3. In this flow, the attenuation of the 
solar energy from the aircraft to ERS-1 is supposed to be negligibly 
small, based on the results by the LOWTRAN S simulation. 


(2) Analysis of ASR data 
1) Check of level of radiance data 


DN values of ASR data were converted to radiance data by the conversion 
table supplied by GER. For the short wavelength infrared (SWIR) region, 
radiance spectra by ASR, IRIS and LOWTRANS 5 are shown in Fig. 4. ASR 

and IRIS spectra show Similar curves and LOWTRAN S is also very similar 


except the regions of water absorption. 


2) Estimation of noise 


GER measured the dark current of ASR by shutting out incident light. 
From these results, noise for ASR is estimated negligibly small in 


comparison witn ERS-1. 








3) Correction of solar energy difference caused by data acquisition 
time 

Data acquisition of ASR were conducted between 10 a.m. and 3 p.m. in 
local time. This caused the difference in incident energy to ASR. so, 


by the use of IRIS measurement of solar energy during the flight, ASR 
radiance data were fitted to that at 10 a.m. in mean solar time. 


(3) Execution of simulation 


After above mentioned procedure, simulation was executed by the flow 
chart in Fig. 3. The 256 (2 X 2X 2X 4X 4X 2) sets of simulation 
parameters were tried for 6 categories listed in Tab. 8. The 
supervisors were selected as they show respectively characteristic 
spectra. Classification was executed for the data showing typical 
spectral features, and also for those showing intermediate features. 
Two examples of the results are shown in Tab 4-a and 4-b. 


(4) Evaluation of the results 


1) Number of bands 


Rate of correctly classified data significantly decreases when the 
number of band changes from 5 to 3. 


2) Band width and noise 


Rate of correct answer is higher for narrower band width for the cases 
of no noise data. However, if the noise for ERS-1 level is added, this 
rate shows the highest value for the band width of 100nm through 150nm. 
This tendency becomes unclear for the case of 3 bands or for the 6 bits 


data. 
3) Data acquisition time 


For the data with noise, the rate of correct answer is lower at 10 a.m. 
than noon. This tendency is more significant for 5 band case than for 


3 band case. 

4) Data with typical spectra and miscellaneous data 
For miscellaneous data, the above-mentioned tendencies become ambiguous. 
5) Number of bits 


Decrease 1N bit number from 8 to 6 is equivalent to noise addition 


(5S) Recommendation for ERS-1/0PS 
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From the above evaluation, following specification 1s recommended: 
1) Band width: 100mn (for the actual S/N) 

2) Central wavelength 

S bands (2.08, 2.15, 2.20, 2.33, 2.38um] 

4 bands [2.08, 2.20, 2.33, 2.38pm) 


3 bands {2.08, 2.20, 2.33pm) 


5 bands are the most recommendable, because the rate of correct answer 
decreases significantly in other cases. 


3) Bait number 
8 bit data are preferrable 


(6) Discussion 


1) In this study, we have not discussed visible and near infrared 
(VNIR) 


2) ASR data were limited to VNIR and 2um region. However, this region 
1s also important. 


3) Actual specification of ERS-1/0PS is slightly different from the 
recommendable specification. So, continuous evaluation effort of this 


specification will be needed. 
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‘able l. Major Parameters of ASR 


VNIR O.4um - 1.00pm 
number of »dands $12 
band width 1.5nm 
detector 31 

SWIR 1.9S52pm - 2.494ym 
numper of bands 64 
band width 8. opm 
detector PDS 

IFOV 
Along track 9.1L72deg 
Cros¢ track 1.720de9 


Interval of Data Acquisition 
O.033sec 10 data are integrated to be 1 record for along track 
(O0.33sec interval for 1 record) 


Aerial Pnoto for Tracking 
One shot for every 10 records by 35mm film 
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Table 2. Specification of Flianht L 
Acquisition " ines and Altitude for ASR Data 
Death Valley 
Abs Alt Alt. fm Spatial No. of Total Line 
(ft) ground resolution. Lines 
(ft) (mxm] (km) 
4,000 2,000 20.1 = 18.3 33 408.0 
10,000 8,000 25.6 x 73.2 15 129.6 
18,000 16,000 32.9%1464 14 185.6 
Total 62 732.2 
Comb Ridge 
Ats Alt Alt. fm Spatial No. of Total Line 
[ft] ground resolution Lines 
(ft) (mxm]) (km) 
6,500 2,000 20.1 =x 18.3 23 259.2 
12,500 8,000 25.6 x 73.2 5 64.0 
20,500 16,000 32.9 = 146.4 12 19.2 
Total 40 342.4 
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Table 


3. 


Parameters of Simulation 


1. Supervisors only (195) , Data including misc data (396) 
2. with / without noise 
3. 8 bits / 6 bits 
4. bandcombination (central wavelength, in ym) 

casel 2.05, 2.15, 2.20, 2.33, 2.38 

case2 2.08, 2.20, 2.33, 2.38 

case3 2.08, 2.20, 2.33, 

case4 2.06, 2.19, 2.34 (RRSS plan) 
5. band width 

50/75 /100/150 (nm) 


6. data acquisition time 
10 A.M./ 12 A.M. 
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Table 4-a. Example of Simulation 


Central wavelength 2.08 / 2.15 / 2.22 / 2.33 / 2.38 


Band width 50 

Bit number 8 

Time 12 AM 
(Supervisors only, without moise) 




















A : Actual Class B: After Classification 
9* : unclassified 

rate of correct answer for all data 100% 

rate of correct answer for limestone 100% 
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Table 4-b. Example of Simulation 


Tab.4-b Example of Simulation 
Central wavelength 2.06 / 2.19 / 2.34 


Band width 100 

Bit number 6 

Time 10 AM 
(Misc data, with noise) 























A : Actual Class B: After Classification 
9* : unclassified 

rate of correct answer for all data 53.5% 
rate of correct answer for limestone 28.0% 
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Fig. 1. ASR Pseudo-Reflectance Data and IRIS Reflectance Data for the 
Site Including Limestone 
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Fig. 2. IRIS Reflectance Data of Limestone and Results of Rock 
Analysis 
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Fig. 3 Flow Chart for the Evaluation of ERS-1/0PS Specification 
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Noise Reduction From Synthetic Aperture Radar Imagery 


430700080 Kyoto Selections From INTERNATIONAL ARCHIVES OF PHOTOGRAMMETRY 
AND REMOTE SENSING in English Vol 27 1988 Part 89 pp 329-336 


(Article by Olli Sirkia, Helsinki University of ,echnology, Institute of 
Photogrammetry and Remote Sensing , Otakaari 1 02150 Tempe Finland, 


Commission VII) 


[Text] ABSTRACT 


The characteristic of the noise in SAR images differs quite much from 
that of more conventional digital spaceborne or airborne images. 
Removal of such usually high-amplitude speckle noise without losing too 
much high-frequency information requires at least some kind of 
adaptiveness from the filter. In this pape. some adaptive spatial 
filters are described and results of their use are compared. 


1. INTRODUCTION 


Speckled noise in SAR (Synthetic Aperture Radar) images is usually found 
as a problem when thinking of interpretation or segmentation and 
classification The noise appears as heavy granularity even in 
homogeneous areas containirg no natural texture. Adaptiveness of the 
filter is required to reduce this granularity sufficiently while 
preserving meaningful edges and objects. These objects can be very 
bright (strong reflectors) because the dynamic range in SAR images is 
ucually much wider than in images in visible or infrared regions. 
Adaptiveness is usually carried out by taking into account in different 
ways statistical or spatial properties of the image and objects on it. 


2. FADING AND SPECKLE 


Inside a ground element (ground cell, resolution ceil) there are usually 
many individual scatterers contributing to the signal received by the 
radar antenna. The signals from these individual scatterers have 
different phases and the summation of the individual signals varies 
because of the relative motion between the radar and the scatterers. 
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Fading 1$ a concept connected with this variation. Fading frequency 
means the interval on the frequency scale after which the total signal 
decreases down to null value. The fading frequency of the imaging 
system depends on the spacing between individual scatterers, speed of 
the radar platform, used wavelength and the distance from object to 
radar. In homogeneous areas containing no dominant individual 
scatterers certain noiselike statistics apply for fading. 


Speckle is the visual effect of those random fluctuations in the return 
Signal caused by the summation of the randomly phased phasors inside one 
ground element. Actually speckle is not really noise but usually 
disturbing and useless information related to particular SAR imaging 
system and target. Speckle can carry some interpretable information 
about the object surface but Gaussianly distributed speckle doesn’t do 
it. In Figure 1 two grey level histograms of different homogeneous 
areas are shown as an example. The areas are from the Seasat-A test 


data used for making this article. 





Fig. 1. Two histograms of different homogeneous areas. 


Multiplicative speckle (noise) model is mostly used because the 
intensity value of a pixel is the sum of many independent random 
variables. It has been found that the signal voltage due to 
backscattering from a distributed target may be described in terms of 
Rayleigh statistics. Wence in homogeneous areas holds theoretically 


dev = 0.523 * mean / N’’?, for emplitude images § (la) 
= mean / N/2 , for intensity imeges (1b) 


where mean is the mean voltage corresponding to the average scattering 
coefficient of a ground cell and dev is the standard deviation 
associated with the measurement (estimate) of that cell. This knowledge 
has been made use of in some of the adaptive filtering methods described 


later. 


The test data used in this article was originally an intensity image 
(square law detection during observing followed afterwards in SAR 
processing stage by incoherent averaging). Later it had been square- 
rooted thus becoming a square-root intensity image. In this case the 
normalized standard deviation of a ground cell measurement should be 
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approximately the same as in amplitude images in which the observations 
are made in linear detection mode before averaging. Here the incoherent 
averaging means averaging of smaller resolution cells to form a final 
pixel in the image and it 1s done to reduce the worst speckle effect. 
Unfortunately this incoherent averaging (increasing the number of 
independent samples per pixel) reduces the spatial resolution 


respectively. 


In the test data here the pixel size is 12.5x12.5m2 each pixel 
consisting of four (N=4) independent samples about 6.6m2 each. From 
this data eight d.fferent homogeneous areas were chosen and the standarc 
deviation to mean value ratios (dev/mean) were computed. The average 
value for this ratio was 0.261 while minimum was 0.230 and maximum 
0.293. The theoretical value would be according to Rayleigh statistic: 
for amplitude images 0.2615 (0.523/41/2) so the ratio values with this 
data seem to fit well to the theory. In homogeneous areas with 
multiplicative noise model the mean to standard deviation ratio (m/s) 
can be used directly to approximate signal-to-noise ratio. This has 
been done later in this paper to compare different filtering methods. 


3. DIFFERENT METHODS 


Quite many studies have been presented about the filtering of 
multiplicative noise both in frequency domain (@.g9., Nomomorphic 
filtering, Wiener filtering) and in spatial domain. In this paper only 
spatial domain methods are described. Often in practice signal-to-noise 
ratio in SAR images is low, too low, @.g., for Wiener filtering to 
achieve good results. The window $size used in the following methods has 
been either S or 7. In adaptive filtering the window size is not very 


critical factor. 


In reference [3] 10 different filters, four of which adaptive, have bee 
compared for noise reduction and classification. There the non-adaptive 
methods were simple averaging and its logarithmized variant to transfor 
the multiplicative noise into an additive one (and exponentiation 
after), median fiitering and iterative median filtering. To some exten 
adaptive but without any speckle model was a method of choosing the 
average value of the most homogeneous subwindow (the one with smallest 


variance among four subwindows). 


One adaptive method was Frost-filtering which assumes that the useful 
information in a small neighbourhood has an exponential autocorrelation 
function. Here the image is filtered with an exponentially shaped 
weighting function and the midpixel value of this window is estimated 


(the weighted sum) 


y’ sk*a* e 3! ix-xoi+ly-yo!l) , 2(x,y) 


where z(x,y) is the image pixel value and (xo,yo) is the location of 
the window’s central pixel. So the weighting. of the window is 
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wix,y) =k*a* Pali 


where k is a normalizing constant (filter weights sum to unity) and 
e = (4/(n*s*?)) * dev? /mean* (n x n - window) 


where dev is the standard deviation of the image pixels inside the 
window and mean is the mean value respectively. So the weighting of the 
midpixel increases when the standard deviation to mean ratio inside the 
window increases. The value of s (reference, assumed real speckle 
¢tandard deviation) can be computed from (1) by setting the mean value 
to 1 for multiplicative noise. It can also be estimated from 
homogeneous areas by collecting (normalized) z/mean values. 


Another method was a maximum likelihood filtering where the most likelv 
estimate y’ for N-look intensity images is the actual positive root of 


the equation. 
y"? - mean*y '2 + N*g, 2 #y' ‘i N*g, 2 #z 20 


where z is the window’s midpixel value in the image to be filtered and 
a,’ = (dev* -mean* *g*) / (1+s? ). 


The third method aims at minimizing the mean square error for each pixe: 
value assuming that the noise is additive with a mean value equal to 
zero. Thus logarithmic transformation is needed before and exponential 
after filtering. Now the estimate y’ for a pixel value is 


_y' log * ™Miog * (2) 5-mean, .)*(dev, 7 “Biog” (dav 


where the speckle variance %, 2 can be estimated from homogeneous 
areas as described earlier or be set values 1.645/N or 0.465/N for 


intensity and amplitude images respectively. 


The fourth method is like the previous one but assuming the 
multiplicative noise model with the mean value of 1. Now the filtering 


midpixel value estimate 


y' = mean + (z-mean)*(1-mean? *s* /dev’ )/(1+s* ). 


The last method was found best in [3] and it is tested in this work too 
named here method 1. Here the value s for speckle was first comruted 
from (la). With N=4 it gives $-:0.523/2. The filtered image (window 
$1ze n=7) and the greyvalue profile of the marked image row are shown in 
Fig. 2 middle row while the original image and the row profile are 

shown on the upper row. The m/s ratio of the whole image increased 
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from the original 1.53 to 1.78 and in homogeneous areas averagely from 
3.26 up to 6.26. So the improvement was much larger in homogeneous 
areas. Also the row profile shows that edges and details are well 


preserved. 


Method 2 in this study 1s introduced by the writer and 1s very much like 
the previous one. Here (la) again 1s used to conclude whether and how 
homogeneous the area is inside the window. From (la) the reference 
speckel standard deviation 1s 0.261. The filtered midpixel value 


y' = mean + k*(z-mean) , where the factor k is simply 
k = | 1-dev/(0.261*mean) | cut into interval [0,1]. 


In homogeneous areas the dev/mean ratio should be near the value 0.261 
and then k 1S about zero and averaging is done inside the window. The 
more the dev/me n differs from 0.261 the more obviously window area 
contains an edsz, details or some texture and less averaging is done. 
With this method the m/s ratio of the whole image was a little smaller 
(1.71) than with method 1 and in homogeneous areas it increased from 
6.26 to 6.61. Also the row profile shows that some details are better 
preserved while e.g., the water area is a little more flattened (lower 


row in Fig. 2). (Figure 2 omitted) 


A variant of method 2 was also tested here where the reftrence noise 
standard deviation was computed from a big water area to get a 
statistically reliable practical estimate. The new reference value was 
then 0.310. The result of filtering is shown in Fig. 3 upper row. The 
m/S ratio in homogeiieous areas increased up to 6.97 and in the whole 
image the ratio increased slightly. So the filter was more effective 
but seemed to blur a little originally unsharp edges. 


A second variant of this method was to use window sizes 7, 5 and 3 for 
every pixel choosing the window with smallest dev/mean ratio. The aim 
was to filter effectively areas also quite beside edges. Unfortunately 
the unreliability of small window statistics caused the loss cf some 
details although, e6.g., edges between water and land areas were 
sharpened. There were no significant changes in m/s ratios. The result 
is shown in Fig. 3 middle row. (Figure 3 omitted) 


Method 3 was presented in reference [4] and is much like the Frost- 
filter described earlier. Here the weighting function is also 
exponential and made narrower in areas with higher dev/mean ratio; 


w(x,y) ask @® qe @'x.y! 


where !X,yi is the distance from a window element to window’s central 
pixel znd the decay constant a is approximated a = dev/mean. In this 
case a sufficiently effective filtering seems to cause more blurring 
than the previous methods. The result with n=S is shown in Fig. 3 lower 


row. 
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Method 4 (5S) is in some sense also adaptive but without noise model. 
This Symmetric Nearest Neighbour filter takes pixel pairs symmetrical to 
the window’s central pixel and chooses the pixel with greyvalue nearer 
to central pixel’s value. Finally mean or median value is computed from 
chosen pixels. This filter tends to flatten local variations (including 
small details) while preserving and even emphasizing more significant 
edge structures. The filter can be used iteratively as shown in Fia. 4 
(n=5). In Fag. 4 lower row is shown the result of applying Symmetric 
Nearest Neighbour filter with n=5 to the filtered image in Fig. 3 upper 
row. The effects of area flattening, small detail weakening and edge 
sharpening can be seen especially from the profile image. (Figure 4 omitted) 


4. SUMMARY 


In this paper some filtering methods for noise reduction in SAR images 
were described. Adaptive spatial filters using loca} first order 
statistics showed fairly good performance. The refersnce standard 
deviation for speckle can be theoretically evaluated or it can be 
estimated from homgeneous areas. In the case of N (number of looks used 
to generate the image) being unknown N can be infered from the 
comparison of theoretical and actual speckle distribution in homogeneous 


areas. 
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(Text) ABSTRACT 


This paper presents the main characteristics of the new generation SPOT 
4 and 5 satellites, the improved version of a s@ries of high resolution 
imaging remote sensing satellites constituting the SPOT family. The 
main objectives of the SPOT 4/5 mission concern: 


i) the continuity of the SPOT 1/2/3 services and products, 


ii) the extension of the high resolution mission, in particular, through 
the addition of a new band in the middle infrared region, 


iii) the development of a new optical instrument called Vegetation, 


characterized by a wide field of view and a high radiometric resolution, 
for the purposes of agricultural forecasts and environmental studies, 


iv) the extension of the satellite life duration, 
v) the capability to offer an embarkment to new passengers, and 
vi) the reorganization and new development of the ground segment. 


1. INTRODUCTION 


The SPOT 1 satellite, the first of a series of high resolution imaging 
remote sensing satellites, has been successfully launched by the ARIANE 
launcher on 22 February 1986. After more than 2 years of operation, the 
SPOT 1 operational behaviour appears quite satisfactory, at the 
exception of some anomalies in the tape recorder operations which 
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resulted ina definitive failure of one of the two tape recorders in 
July 1986. Moreover, the SPOT 1 image quality fulfills the original 
system specifications as explained in Ref. 1. General information about 
the SPOT 1 system may be found in Refs. 2 to 5. 


The basic mission of the SPOT program being to establish a continuity of 
the SPOT data and services, it has been decided: 


--to produce two recurrent models named SPOT 2 and SPO7 3; 


at the time of writing (May 1988), SPOT 2 has completed its integration 
phase and is planned to be launched in early 1989. The industrial phase 
of SPOT 3 has already started and SPOT 3 should be ready to be launched 


in early 1990. 


--to develop two improved satellites named SPOT 4 and SPOT 5S. 


The desicn phase (phase A) of SPOT 4/5 is over, and its preliminary 
industrial definition phase (phase 8) will be completed at the end of 
1988. The industrial development and realization phases (phases C and 
D) will start in early 1989, the goal being to be ready to launch SPOT 4 
at the beginning of 1993 and SPOT 5 in 1995. 


The main objectives of the SPOT 4/5 mission are the following: 
--provide the continuity of the SPOT 1/2/3 services, 

--ensure the extension of the high resolution mission, essentially 
through the addition of a new spectral band in the middle-infrared 


region and the on-board registration of the 10m/20m resolution channels 
in the high resolution instruments HRVIR, 


--add a new payload named Vegetation, a wide field of view, low spatial 
resolution instrument for high repetitivity observations in the visible 
and middle infrared range, 


--increase the life duration of the satellite, 


--offer new possibilities of embarkment of passengers, 


--improve the ground segment. 


This paper describes these objectives, the main technical choices 
adopted for the space and ground segments, and their evolution relative 


to the SPOT 1 system. 
2. CONTINUITY OF THE SPOT PROGRAM SERVICES 


The SPOT 4/5 system must provide the continuity of the SPOT services up 
to the year 2000 at least. 
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The SPOT 4/5 system appears as an extension of the SPOT 1/2/3/ system, 
hence the main system parameters have remained unchanged. Thus the 

orbit parameters, the local time and the main high resolution instrument 
characteristics (spatial resolution, swath, spectral bands 1n the 
visible and near infrared region, etc...) have not been modified. 





The continuity of the SPOT 4/5 system is also apparent at the ground 
Station interface level. The ground stations belonging to the SPOT 
1/2/3 system will be able to receive and process the SPOT 4/5 telemetry 
emission without any other modification than proper equipment 
adjustment. In other words, it will be possible to program the SPOT 4 
satellite in a mode compatible with the SPOT 1 modes with a telemetry 
format identical to that of SPOT 1 (see also section 3.3.). 


Moreover the preprocessing of raw data as well as the use of system 
parameters to elaborate higher level products will be made possible with 
minor modifications with the type of hardware and software currently 
used in the ground stations belonging to the SPOT 1/2/3 system. 


3. EXTENSION OF THE HIGH RESOLUTION MISSION HRVIR PAYLOAD 


The SPOT 4/5 system will allow an improvement of detailed thematic 
inventories through the addition of a new band in the middle infrared 
region (this is why the HRV instruments have been renamed HRVIR, 
standing for High Resolution in the Visible and Infrared Range) and the 
ability to register on-board the 10 meter resolution and the 20 meter 


resolution channels. 
3.1. Use of a middle infrared band (MIR) 1.58-1.75 microns 


In agronomy and natural vegetation studies, the MIR band will permit an 
estimation of the ground humidity and will improve agricultural 
inventories, particularly those of cereals. The objective in this 
latter domain is to increase the performance of cultures identification 
with a monotemporal of multitemporal set of data, first, by bringing an 
information on the differences of water content in the vegetation cover 
either in time or between different parcels of a given agricultural 
zone, and second, by offering a better sensitivity than that of the near 
infrared channel 83 to the covering rate of the ground by the 


vegetation. 


In hydrography, it will be possible with the MIR band to perform a 
Getailed cartography of the snow mantle (since this band discriminates 
@asSily snow and clouds) and to inventory the various types of snow. 


104 





In geology, the 1.58-1.75 microns band can improve significantly the 
possibilities of discrimination between various types of minerals; it 
can also give an access to the notion of ground humidity, and indirectly 


to anomalies of metal content of the ground. 
3.2. Registration of 10m resolution and 20m resolution channels 


The panchromatic 0.49-0.73 microns SPOT 1 band operating at 10 meter 
resolution has been replaced on SPOT 4 by the 82 band (0.61-0.68 
macrons). The 10m resolution channel is of strong interest in urbanism, 
forestry, small parcels inventory, texture analysis, etc.... Many 
applications of the 10m resolution channel require its combination with 
the 20m data. The adoption of the same band 82 for the 10m and 20m 
resolution channels has made possible an instrument concept where bpcth 
channels are registered on-board, which is most convenient in the 


production of data for such applications. 
3.3... Main characteristics of the HRVIR payload 


The HRVIR payload contains two identical HRVIR instruments. As those of | - 
SPOT 1, the instruments are pointable in the across track direction up ae 
to plus over minus 27 degrees (plus over minus 460 km on ground) to 
provide accessibility and stereo coverage capability. The HRVIRs 
operate independently with the capability to produce monospectral (M), 
multispectral (xX), and/or SPOT 1 compatible multispectral (XS) images 


whose swath width is 60 km at Nadir. 


In the monospectral mode, the HRVIR generates a black and white image in 
the B2 band with a ground resolution of 10 meters by 10 meters. 


In the multispectral mode, the HRVIR generates a multispectral image in 
4 bands 61, 82, 83, MIR whose ranges are 0.50-0.59, 0.61-0.68, 0.79- 
0.89, 1.58-1.75 microns respectively at 20 meter resolution. 


Q@ 


In the SPOT 1 compatible multispectral mode XS, the HRVIR generates a 
multispectral image in 3 bands 81, 82, 83 at 20m resolution. 


The X and M data are compressed on-board using Delta Pulse Code 
Modulation techniques (8-5-5-8), while the XS mode is not compressed on- § 
board. The pixel is restituted on ground on 8 bits. ™ 


105 








Each HRVIR can deliver two 25 mesabits/sec. data streams corresponding 
to the monospectral (M) and multispectral (xX or XS) modes. Only two 
among the 4 bit streams available are transmitted to the ground or 
recorded on-board, the total telemetry data rate being 50 megabits/sec. 
The recording capability of each of the two tape recorders has been 
increased from 29mn (SPOT 1) to 40 mn on SPOT 4. The transmitted data 


are ciphered ori-board. 


The radiometric and geometric image quality criteria are listed in Table 
1. Note that the instrumental radiometric resolution and dynamic range 
are obtained from the Table 1 ground reflectance data and from the 
assertion that the solar incidence angle be in the range 15-60 degrees. 
Note also that while tne localization and length alteration criterie2 
apply to rectified images using system corrections only (no ground 
control point), the photogrammetric restitution criterium applies to the 
processing of stereo pairs with the use of ground control points. 


We mention finally two interesting improvements of the SPOT 4/5 system 
relative to SPOT 1. First, it will be possible to program each HRVIR 
independently regardless of any possible motion of the pointable mirror 
of the other HRVIR, whereas on SPOT 1 such motion on one HRV invalidates 
the image acquisition of the other HRV due to pointing inaccuracies, 
resulting in significant programming difficulties. Second, there will 
exist at any time 3 operational analog gains per spectral band and per 
HRVIR; the SPOT 4/5 mission center will program the appropriate gain out 
of the 3 gains available on different segments of a given orbit from 
criteria elaborated in advance with an histogram data base made of all 
histograms of cloud free scenes acquired by SPOT 1/2/3. Such a measure 
should be most appropriate in bands 81 and 82 where the scene 
radiometric dynamic range may vary significantly with season, latitude 
and the type of scene considered (sand, vegetation, etc...). 


4. ®& NEW PAYLOAD: VEGETATION 
4.1. The VEGETATION mission 
The HRVIR instruments provide very detailec images of the Earth surface 


but their repetitivity of observations is sometimes insufficient for 
applications relative to renewable resources or environmental studies. 
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Table 1. HRVIR Image Quality 
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The Vegetation instrumert, named also VGT, is conceived as a high 
radiometric resolution instrument with a wide field of view, with 
consequently a high repetitivity observation capability allowing a 
continuous monitoring of evolutive phenomena. The first priority 
objective of VGT is the worldwide operational monitoring of evolutions 
of crops and natural vegetation for the purposes of agricultural 
forecasts and environmental studies. The second priority objective of 
VGT is the observation of oceans mainly for scientific purposes. These 
two objectives are detailed hereafter. 


4.1.1 Estimate of crops and natural vegetation monitoring 

The possibility to estimate crop yields relies on the capability of vVGT 
to monitor temporal evolutions of the ground reflectance. Hence the VGT 
instrument will deliver data which will eventually permit: 


i) to retrace at any point the evolution of biomass and of the vegetal 
cover structure to detect possible growth anomalies or hydrous stresses, 
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11) to cartography the disparities of evolution, ina given 
agricultural or climatic zone, by a multidate measurement over this 


entire zone, 


211) to provide an indication of the amount of biomass at its maximum. 


The VGT data will also permit, besides the biomass evaluation, the 
estimation of the proportion of given species in large homogeneous 
regions. Such estimate, wherever possible, should be performed through 
a spatial extension of detailed classifications elaborated from HRVIR 


data. 


These various informations will be obtained with a combination of the 
VGT informations at any point in the bands 61, 82, 83, MIR, the same 


bands as those of HRVIR. 


4.1.2. Oceanography 





The addition of a “blue” band BO (0.43-0.47 microns), centered on the 
main absorbtion band of phytoplankton pigments, makes possible the 
monitoring of the amount of vegetal biomass in water, with the following 


potential applications: 


i) cartography of spatio-temporal variations of the distribution of 
phytoplankton and archive, 


ii) quantitative estimation of worldwide primary production of 
phytoplankton, 


111) climatologic and biogeochemical studies (carbon fixation, CO2 
cycle), 


iv) use of biomass as a tracer of the dynamics of superficial water 
masses, 


v) assistance to fishing by fast localization of zones of high 
phytoplankton concentration. 


The bands B81 and 82 can also be used, in addition to BO, to perform a 
mesoscale monitoring of water turbidity, in particular along the coast 
lines. All these applications require careful atmosphere corrections 


which will be taken care of by the system. 


4.2. The VGT operations 


The VGT instrument can generate 2 channels at $00 kbits/sec. data rate 
each. 
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The first channel is called "direct" or "regional’ channel and transmits 
data in real time with a spatial resolution at nadir of about 1 km. 

This operating mode is called “Regional Observation" and provides 
information from the region in which the station (Toulouse or another 
ground station) is located. The working zone for the Regional 
Observation is subject to the condition that the sun incidence angle be 


lower than 60 degrees. 


The second channel is called “recording” or "worldwide" channel and 
transmits data which have been recorded at 4 km resolution at nadir. 
The 4 km resolution is obtained by agglomerating onboard 4x4 
neighbouring image elements at 1 km resolution. The production of 4 km 
resolution recorded data is called “Worldwide Observation’ and can be 
obtained only at the Toulouse station. The working zone for the 
Worldwide Observation is nominally located between 60 degrees North and 


40 degrees South. 


A very high repetitivity of observations is achieved, taking into 
account the orbit characteristics, by a wide field of view of the 
instrument (+ over - 50 degrees) corresponding to a swath width of 2200 
km. With this configuration, a given point on ground 1s seen every day 
by VGT above a latitude of 35 degrees, and seen 3 days out of 4 at the 


Equator. 
4.3. The VGT instrument 


The VGT instrument contains a radiometer made of 5 objectives (one for 
each spectral band 80, 81, 82, 83 and MIR) and an electronic chain which 


amplifies and codes on 10 bits the signal. 


An array of 1728 CCD detectors is located in the focal plane of each 
ob jective. 


The main characteristics of the VGT payload are displayed in Table 2, 
while the VGT geometric and radiometric image quality criteria are 
listed in Table 3. 


Note in Table 3 that the local distortion on an image remains 
sufficiently small that 2 images taken at different dates be 
superimposable to a 1 pixel rms accuracy and to within a translation of 
an integer number of pixels (less than 5 by virtue of the localization 


performance). 


S. EXTENSION OF LIFE DURATION AND MODIFICATIONS OF THE SATELLITE 


A number of improvements have been inserted in the satellite design, in 
order to improve its performances and increase its life duration from 2 


years (SPOT 1 concept) to 5 years. 
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Table 2. VGT Main Characteristics 





Weight 160 kg 
Volume 0.7xn121 23 
Power 200 watts 
Data rate 2 x 500 kbits/sec. 
Pixel codiixy 10 bits 
Telemetry frequency 8.2 Giz 
Recording memory 256 Megabits 














Table 3. VGT Image Quality 
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The increase of life duration will be ensured by an appropriate choice 
of elementary components and by a new thermal design in order to lower 
the average operating temperature and hence reduce the failure rate of 


the “lectronic components. 


The flexible solar arrays of the solar generator have been replaced by a 
solar array with rigid panels, a concept both simpler and more reliable, 
and the available power end of life has been extended from 1100 watts 
(SPOT 1) to 1800 watts. Moreover the 3 x 24 &.h batteries have been 
replaced by 4 x 40 A.h batteries to accommoda:te to a higher demand of 
power and to provide a reduced depth of discharge. 


The mechanical concept has beer changed by adding an “equipment case" 
between the platform and the payload. Such equipment case is needed to 
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accommodate the new tape recorders and to provide a mechanical support 
to the HRVIRsS and VGT payioads as well as to a number of passengers of 


the SPOT 4/5 platform. 


The in orbit configuration is shown on Figure l. 
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Fig. 1 SPOT 4 in Orbit Configuration 


6. THE SPOT-4 PASSENGERS 


The extended capabilities of the SPOT 4 satellite allow to embark the 
following additional experiments: 


--PASTEC (PASsager TECnologique) which is a group of experiments for 
studying orbital environment, 


--DORIS (Doppler Orbitograpnry and Radiopositioning Integrated by 
Satellite) which will be used for precise orbit determination and 


localization of ground platform, 


--a Radar Transponder which will be used to calibrate the ground radar 
of the Guyana Space Center and to train the operators, 


--PASTEL (PASsager Spot de TELecommunication Laser ) 
The PASTEL package is part of an European Space Agency program for the 
experimentation of an optical link between a low orbit spacecraft (SPOT 


4) and a geostationary satellite. This optical link will be used for 
transmitting the SPOT 4 payload data stream at arate of 2 x 25 Mb/sec. 
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These data will demonstrate the feasibility and performance of opiacal 
communication and will be used on a pre-operational basis to relay image 


data. 


Using these facilities will bring about a fundamental change in the 
Capabilities and management of the system. AS an example, Figure 2 
shows the visibility on earth between SPOT 4 and a geostationary 
satellite located at 19 degrees E. As seen in Figure 2 the visibility 
is restricted mainly by the Earth and the SPOT 4’s body which may 


obscure the optical link. 


The operational SPOT services will still use the tape recorder and 
ground station network, while the optical link will be used to extend 
image transmission capabilities in the event of lack of tape recorder 


capability or failure. 


























- a 

= , 

e 4 
* 

‘ . 4 
) de 4 
i ; 
- —— ‘ ' 
2 ee e. 140. a 


wa. 2. ee. we rr. 
-~E™ gg. O25 S= 
7 om tsmurie Gy mem 
Fig. 2 PASTEL Ground Coverage 
7. THE GROUND SEGMENT 
7.1. The image ground segment - Products 


Two stations--Toulouse and Kiruna--will be used to ensure a large direct 
and playback transmission capacity of the HRVIR data. Night passes and 
some day passes will be used to collect recorded VEGETATION data on the 


Toulouse station. 


Apart from these stations direct read-out X band stations can be 
implemented worldwide to acquire HRVIR and VGT data. 
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The Toulouse and Kiruna stations contain an HRVIR Archiving and 
Preprocessing Center which archives all the HRVIR data and processes 
level O and level 1 products. These centers are also capable to procet: 


images issued from the SPOT i1/2/3 archive. 


The level O product is an exchange level product with geographically 
arranged pixels and a raw radiometry. 


The level 1 is a radiometrically corrected (level 1A) and geometrically 
corrected (level 18) product. Only system corrections are applied to 
perform the geometric rectification. Note that oontrary to SPOT 1 the 
level 18 scenes can be mosaiqued along the track without further 
processing and centered upon request at any point along the track. Note 
also that the cloud cover notation at the inventory step (catalogue) of 
the image processing will be significantly improved relative to the SPOT 
1 system, by the addition of automatic and semi-automatic analysis 
methods (analysis of the reflectance and reliability, graphics 


facilities, etc...). 


The Toulouse station includes also a Preprocessing Center of VGT images 
which produces and archives level 1 data. The level 1 VGT product is 
characterized by radiometric and systematic geometric corrections, as 
well as a systematic correction of atmospheric effects. The VGT level 1 
product contains also reduced data such as the vegetation or turbidity 
indexes, and a simplified classification map (snow, cloud, water, 


ground). 


The products delivered by the Toulouse VGT Center will be available 
within 8 hours in the Regional Observation mode, 36 hours in the 
Worldwide Observation mode after their acquisition by the satellite, as 
soon as VGT enters its fully operational phase. 


7.2. Organization of the ground segment 


In the operational phase, after the post launch assessment, the 
responsibilities will be shared between the various entities as follows. 


The satellite operator (CNES) will be in charge of the Control Center 
and the Platform and Payload Management Center. These centers perform 
all satellite management functions, monitor the satellite operations, 
generate the necessary commands and determine the orbital parameters. 
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The commercial operator (SPOT IMAGE) has the exclusive license to 
organize the promotion, the distribution and sale of SPOT data. It will 
be in charge, moreover, of the Mission Programming Center and of the 
Toulouse image ground segment. The Mission Programming Center is 
responsible for mission management regarding the HRVIR and VGT payloads. 
It programs the on-board image acquisition aia storage in accordance 
with user requirements, and also sends the dat# necsssary for image 
acquisition and processing to the ground stations. 


8. CONCLUSION 


CNES has undertaken an ambitious earth observation and data distribution 
program. This ambition is obviously already apparent since SPOT 1 has 
been giving most satisfactory results for more than 2 years of 
operational life, and since the following satellites SPOT 2 and SPOT 3 
have been either built or fully decided. The aim now is to ensure a 
continuity and an improvement of the SPOT services and products up to 
the year 2000 and beyond. The SPOT 4/5 system has been designed to 
fulfill this requirement. This generation of satellites represents a 
technological challenge (in particular, the ca-slopment of CCDs in the 
MIR band, the wide field of view objectives with CCD arrays of 
Vegetation, the development of a new polar platform) and widens the 
field of potential applications of earth remote sensing, especially 
through the use of the new payload Vegetation. 
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[Text] Abstract 


The concepts in algorithmic design and operational use of the German 
Processing and Archiving Facility Geocoding System for ERS-1 are 
presented. The System will produce non-terrain corrected and terrain 
corrected products. The latter require complex algorithms in terrain 
height consideration and ground control pointing. 


The Geocoding System will be implemented on UNIX workstations and 
programmed within an image processing system. To guarantee operational 
throughput and flexibility, new software tools will be implemented. 


2. Introduction 


In 1990, the first European Remote Sensing Satellite (ERS-1) will be the 
first step towards a scenario of satellites carrying a Synthetic 
Aperture Radar (SAR) as a prime sensor (SIR-C, RADARSAT, J-ERS-1). For 
the 90’s and beyond, SAR-sensors--then upgraded in terme ~f wavelength 
and polarizatiorn--will be the basic source of image information besides 


optical sampling techniques. 


To cope with the large amount of data and to benefit from ite intrinsic 
information, the data procassing and data evaluation techniques must 
keep pace with the develorrant of space-segment hardware. 


Especially for SAR-data processing, powerful hardware and sophisticated 
software concepts will be realised, to generate digital raster image 
Gata out of the holographic information gathered by the sensor (Noack, 


1987). 
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The microwave image information is now basically interpretable by 
digital or visual means, but suffers from the lack of geographic 
context. Therefore, scientists, engaged in remote sensing, first try to 
rectify or geocode the satellite images (for optical images as well) to 
ease the interpretation and to utilize further existant geocoded 
information for daa synergism, covering several aspects (spectral 
wavelength, time, polarization, look-angle) of the investigation area. 


Rectification and geocoding algorithms are more or less a standard 
feature of many remote sensing image processing systems, but are time- 
consuming (by manual tiepoint search and raster data resampling) and 
produce often unsatifactory and uncomparable (with other systems) 


results. 


This is one reason, why many potential commercial and governmental users 
hesitate from using satellite borne images due to the complex 
proceuures, which must be performed, before getting useful information. 


For optical scanner images, a few national remote sensig organisations 
and commercial distributors are already offering geocoded satellite 
images. For the future microwave satellite scenarios mentioned above, 
such geocoded products should become a standard user service and is 
therefore planned by several facilities (JPL, CCRS, ESA). 


The German Aerospace Research Establishment (DFVLR) is currently 
building up--in contract with ESA--a Processing and Archiving Facility 


(PAF) for ERS-1 SAR images. 


One major sub-system besides the processing of the radar raw data by the 
new Intelligent SAR processor ISAR is the SAR Geocoding System GEOS. 
The major design aspects of the Geocoding sub-systems will be described 


in the following chapters. 


3. Algorithmic Needs for Future Satellite SAR Image Data Geocoding 
Systems 


Tne future remote sensing satellite scenario can be characterized by the 
amount and variety of data. The SAR processing systems must cope with 
this data by increasing the processing speed up to near real time. The 
variety of the data for different satellites and sensors iust be 
properly obeyed in the SAR cases, where orbit, wavelength and satellite 
technical parameters determine the SAR image generation process and its 
quality. Thus, the input data needs a flexible structure in hardware 


and software of the SAR processing system. 


The target where all SAR and also optical sensors can be referred to, is 
the imaged body itself--the earth. All seisor views of the same 
physical ground element and all sensor data can be geometrically re- 
ordered to a pre-defined cartographic presentation. The latter task is 


subject of Geocoding Systems. 
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Satellite image Geocoding Systems can be characterized by the following 
attributes: 


--kind of sensor and its imaging physics 

--parametric or tiepoint rectification (or intermediate) 
--consideration of terrain height displacements 
--support of different cartographic projections 


--data throughput capabilities 


The imaging physics (which intentionally substitutes “geometry” in this 
case) of Synthetic Aperture Radar sensors determines the geocoding 
algorithm to compute the proper geo-location of each output pixel. The 
orientation of one image pixel in a SAR image and its correspondence to 
an earth location can be established by solving three equations 


(Curlander, 1982): 


--The range timing (the absolute travel time between sensor to target) 


--The azimuth Doppler reference (depending on the SAR Processing 
algorithm applied) 


--The actual terrain shape (ellipsoid or Digital Elevation Model if 
terrain height is considered). 


In remarkable contrast to optical sensors, the pointing accuracy of SAR 
systems do not depend on angle measurements, but on time and frequency 
analysis of the radar echo. These parameters can be extracted with 
better accuracy than satellite angle control. 


Having the necessary parameters at hand, a geocoding approach with 
parametric rectification using the sensor/processor entities is obvious 
for SAR systems and gives sufficient results for many application areas. 


Unfortunately, SAR images and its geometry are affected by the local 
earth surface elevation more heavily, than it is the case for sensors. 
For ERS-1 (look angle 23 degrees at mid swath) 1000 m elevation come 
along with mis-registration of about 425 m. For areas with moderate to 
hilly terrain, these mis-registrations can cause pixel geocoding errors 
which are considerably higher than the possible system pointing accuracy 
mentioned above. The most common way to tackle with this problem is to 
introduce the terrain height of each earth point target by using a 
Digital Elevation Model (DEM). 


In principle, an absolute mapping accuracy could then be achieved within 
the bonds of the system parameter accuracy. However, the introducing of 
a "foreign data set like DEM, which has no direct relationship to the 
SAR mission data or even to remote sensing geometry, can cause “scaling” 
problems. These problems will come from: 
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--the remaining uncertainty of the SAR pixel location. 


--shifts and scalings in reference to the common ellipsoid. 


--errors and uncertainty in DEM generation. 


Having this in mind, it 1s obvious, that some ground control must be 
incorporated in the geocoding process, when precision mapping is 
aspired. In any case, the quality control of geocoded SAR products must 
use ground control to verify the achieved accuracy. 


Following the principies of geocoding stated so far, SAR geocoded 
products can roughly 66 divided into two categories: 


--Geocoded without terrain correction, using an ellipsoid as reference 
earth surface and referenced to the system initial geometric accuracy. 
These products are abbreviated here GEC (Geocoded with Ellipsoid 


Correction). 


--Geocoded with terrain correction, using a DEM as reference earth 
surface with more precise accuracy introduced by ground control points. 
These products are abbreviated GYC (Geocoded with Terrain Correction). 


The algorithms for generating these two basic products at the German-PAF 
are described in the following chapters. 


4. Geocoding Without Terrain Correction 


Due to the spare availability of DEM data, the GEC-data are considered 
to represent the major throughput of the Geocoding System. As already 
mentioned above, the SAR systems inherent pointing accuracy can be used 
to rectify SAR images without any foreign control or tiepoint 
information. Using this scheme, the geodetic location of each SAR image 
pixel could be estimated and resampled to the desired map projection. 

Of course, this approach would be rather time consuming and ineffective. 
Instead, the correct SAR image pixel to ground correspondence of four 
corner pixels of an appropriate rectangular frame is calculated 
precisely. The area and pixels in between the rectangle are transferred 
using a lower degree polynominal, estimated from the four corner pixels 


mapping function. 


The deviation of the geocoded pixel location using polynominal 
interpolations instead of the correct transformation formula depends on 
the mesh-size and the used interpolation formula. For SEASAT images 
this approach, with a 500 pixel mesh size gives only deviations of one 
pixel with respect to correct mapping. Due to the increased data 
throughput without considerable loss of accuracy, this rectification 
scheme is used in many geocoding algorithms. 


Considering that geocoding of satellite borne imagery consists mainly of 
a rotation towards grid nerth and some rectification of internal 
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distortions of the image, the rotation of points within the mesh frames 


ne (#) - (sors 223) (3) 


In re-arranging the transformation martrix and splitting in two steps, 
the formula reads (Curlander, 1987a). 


(3) = (ate sad) (og* #*)(3) 


A further initial step can be involved to avoid aliasing when rotating 
and resampling the image. This step oversamples the image in y- 
direction. The oversampling factor nf is derived from the rotation 


angle (Friedmann, 1981). 
nf = 1 + tane 


Although the initial l-step formula seems to be more convenient than the 
derived 3-step formulas, the re-arranged resampling scheme has important 
advantages. The step-by-step rectification touches only one coordinate 
direction per step and needs only i-dimensic‘%al resampling kernels. 
Considering a 4 x 4 cubic convolution 2-dimensional kernel, 16 points 
must be touched, whereas only 8 in the modified l-dimensional scheme 
(plus 4 additional in the initial oversampling step). Tie remarkably 
increases data throughput. The line/column access scheme fits to the 
Gata needs of array-processors and simplifies address calculation. 
Figure 1 shows a schematic sketch of the geocoding process. 
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Fig. 1 Geocoding Process Without DEM Correction 


Because no ground control information or control and merging of DEM data 
is required for this geacoding process, it can run without operators 


interaction. 


Therefore, the flexible architecture of the D-PAF Intelligent SAR 
Processor (ISAR), consisting of a huge hardware memory, a powerful 100 
MFlop array processor and a UNIX workstation host, is used as an 
implementation target for generating the SAR geocoded GEC product in one 
standard map projection (e.g.: UTM). 


The direct coupling of SAR image generation and SAR geocoding within one 
process will make it possible to use the complex intermediate data for 
geocoding input rather than the look summed and detected slant range 


image. 
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S. Geocoding With Terrain Correction 


When terrain correction is involved in the geocoding process, the 
computing algorithn becomes more complicated, due to additional raster 
data sets needed (DEM), which requires data handling and geometric model 


adjustment. 


Together with contracted institutes, DFVLR has performed a phase-8 study 
to identify the required algorithms and to introduce and test new 
approaches. According to the gained results, the operational processing 
scheme for GTS-products will be as follows: 


The terrain elevation data, which is necessary for GTS-product 
generation will be stored on-line in raster format ina large disk 
database. Before adding DEM data to this database, quality tests will 
be performed to ensure a minimum accuracy of the DEM data used. 
Although, it must be aspired that DEM generating organizations might 
perform a quality test, the results we have gained with different height 
data form different sources, unveiled a variety of errors--some severe 
and some neglectable for geocoding--which are partly derived from the 
generation process of the DEM data. 


From the DEM database, the required part of the DEM is resampled to the 
cartographic projection demanded by the geocoding process. This 
resameling of DEM data will use a comparable fast algorithm as outlined 
for the SAR GEC-Geocoding mentioned in the previous chapter. 


As stated earlier, the utilisation of the resampled DEM for terrain 
correction might give erroneous results, because the estimated SAR point 
location accuracy (ERS-1 about 100 m) and the DEM accuracy (FRG: approx. 
25 m raster size) might not fit on a pixel by pixel basis. Hence, both 
geometry models need to be adjusted and refined. Due to the origin of 
most DEMs from maps, the height model itself will be used as a geometry 


adjusting reference. 


This can be achieved by considering the strong influence of hilly 
topography to the backscatter of SAR images. Using the imaging geometry 
and simple backscatter assumptions, a simulated image can be generated 
out of hilly parts of the DEM, showing the backscatter behaviour of the 
SAR with an exact cartographic reference. After some pre-processing, 
the original SAR image can be correlated with the template, giving an 
“artificial” tiepoint. The phase-8 study results (Wiggenhagen, 1988) 
showed that the process can work automatically and gives better results 
for considerably good DEMs than matching SAR with Landsat TM images. 


The “artificial” tiepoint with DEM-derived information has several 
advantages: 


--No other external data set (like Landsat or Spot) is needed, which 
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might antroduce inconsistent geometric adjustment. Instead, the DEM 
data itself, which is used for terrain correction is the geometric 


reference. 


--The search for “candidate” hilly regions for simulation and matching 
can be done in advance. Using the pre-known geometry of ascending and 
descending ERS-1 orbit, candidate regions can be found and flagged in 
the DEM database. Operationally, it will oniy be necessary to simulate 
small regions, around hilly areas, rather than simulating the whole 


image. 


--Test show (Curlander, 1987b) that the DEM to SAR image geometry 
derivation 1s mainly driven by linear effects. Thus, only few tie- 
points (2-3) will be needed for the adjustment, which can be found even 


in moderate hilly terrain. 


By introducing these few tie-points, an adjustment calculation will fit 
the geometry estimation of the processor to the needs of the DEM data 

and cartographic accuracy (Raggam, 1987). The number of necessary tie- 
points and the Jegree of adjustment will depend on the sensor benaviour. 
For ERS-1 the orbit will be known very accurately and the SAR processor 
will produce considerable location estimates. Thus, only a few tiepoint 


might be necessary. 


Having the adjusted imaging model and the resampled DEM data, a backward 
resampling to generate the GTS-product can start. There, each DEM 
point, corresponding to the pixel size in the final output image, is 
transferred to a common cartesic system, as well as the satellite crrit. 
A slant-range vector from DEM point to an estimated orbit position is 
Gefined, giving a relationship between slant-range length and Doppier- 
azimuth reference. The orbit point 1s now iterated to fulfill the 
range-Doppler relation given by the system geometry. The calculated 
range and orbit position gives the ‘ine and column location in the SAR 
Slant range image which is used to resample the appropriate gray value 
of the output pixel of the GTS-product, located at the initial DEM point 
(Meier, 1985). Figure 2 shows a schematic flow chart of this process. 
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Fig. 2 Geocoding Process With DEM Correction 


In contrast to the GEC-product generation, the GTS-product generation 
must touch each DEM point for the exact location computation of the 
geocoded SAR pixel. Hence, the utilisation of fast geocoding 
algorithms, mentioned earlier, is more complicated to realise. A poly- 
nominal interpolation of the correct imaging geometry for small sub- 
frames must also treat the DEM geometry transfer and terrain height 
dependant adjustment. This resampling scheme will increase data 
throughput and is modified to meet the requirements of pixel location 
accuracy within the magnitude of the DEM raster size used. 


6. The Software Environment of the Geocoding System 


The sketch of the geocoding algorithms explained in the previous 
chapters reveals the diversification of the necessary data and 


input/output interfaces for this task. 


Ooviously, geocoding is tightly coupled to the output of the SAR 

processor and its computed parameters. As far as mors accuracy in 
location is aspired, geocoding has also to tackle with digital elevation 
data, ground control points and precise map projections. Finally, all 
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these data sets must be visualised, at least for quality control. For 
the D-PAF, geocoding is a set of modules to be integrated in a remote 
sensing image processing system, which is modified to the special needs 
of an ooerational production environment. Since several years, DFVLR, 
together with a software company, develops and uses the UPSTAIRS image 
processing software running on the VAX-VMS operating system. The 
UPSTAIRS software is dedicated to process remote sensing images and 
supports program development with pre-defined user/machine interface 
calls and simple access to a sophisticated image database. A magnitude 
of application programmes covers the area of radiometric and geometric 
processing of satellite borne raster images. 


For the geocoding system, the UPSTAIRS <ernel was vcorted to UNIX and Ae 
partly rewritten in C (see Figure 3). Along with the software transfer, 
some important modifications, which enhances operationality of the 
Geocoding System have been integrated. These upgrades support 

operational data throughput, new data like-DEM and GCP and a new user 
interface based on the bit-mapped screen of UNIX workstations. They are 


described in more detail as follows. 
Data Throughput Upgrades and Upgrades for Operational Use 


In contrast to an image processing systems as a “scientific tool,” the 
operational Geocoding System has to store only the SAR input and the 
geocoded output image as a minimum. This data must be accessed and 
transferred as fast as possible. fYfhersefore, software tools are 
integrated, which manage large pixel sets for fast data I/O by 
considering disk resources and IIW/SW capabilities, especially the data 
transfer to the attached array processor. 
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Fig. 3 Geocoding Image Processing System Software Layers 


To exchange data and parameters between different modules of the 

geocoding system, a “parameter pool” is integrated, which resides in the 
computer memory and delivers access via write and read functions. This 
feature will support the programming and the flexibility of the modules, 


rather than exchanging parameters by files. 


The Databases for DEM 


Although the Digital Elevation Model raster data could be stored like 
image data, its handling within the Geocoding System is organised ina 
DEM database. The database is structured according to a geographic 
binary reference tree. Each branch in the tree defines a reference 
area. The final data-tile on the earth surface is able to store DEM 
raster sizes of about 10 m x 10 m resolution plus additional annotation 
data. These include flags for candidate simulation areas for artificial 
tie-points stated earlier. All DEM data will be stored on-line on 
magnetic disks. The capacity of the DEM database is only limited by 
Gdisk capacity and--in principle--allows to store DEMs covering the whole 
globe. With subroutine calls, the DEM data can be accessed and 
resampled to the format required by the geocoding process. 
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The Operators Interface for the Geocoding System 


The geocoding system will run on a UNIX workstation with “bit-mapped" 
screen. This screen allicws to display concurrently graphics, images (up 
to 256 different colours out of 16 million) and alphanumerics in 
different “windows” and utilizes a “mouse” as a fast and versatile user 
input device. Whereas the normal menu screen is managed in 
alphanumerics and gets its input via kaydeard, the system will have some 
tasks, which directly benefit from the advanced user interface. Among 
these tools in the final image quality control, which will be performed 
using the normal colour screens, capable of displaying 900 x 1152 pixels 


in 256 shades of grey. 


There, features in the geocoded SAR image--displayed on the screen--will 
be compared with homologue points on map paper sheets on a digitizer 
table. Knowing the geocoded coordinates of the product the cursor con 
the screen will move, controlled by the movement of the digitizer pen. 


7. ##The Hardware Environment of the Geocoding System 


The revolutionary development in computer techniques gave smaller 
systems the capabilities of former mainframe hosts. Due to their 
reasonable costs, small workstations can be dedicated to one task and 
equipped with special hardware like bit-mapped screens, fast and large 


disks and array processor. 


Nevertheless, such dedicated systems are not isolated as far as they are 
coupled via communication links to other workstations, building up a 
computing system with one (or even more) computing units for each task, 
but with the transparency of a large mainframe. Together with the UNIX 
operating system, supported by more and more manufacturers, these 
computer networks can easily be adopted to different needs. 


The Geocoding System is part of a larger workstation network with 
Gedicated software and hardware for each function like SAR processing, 


data management and facility control. 


A sketch of the Geocoding System is given in Figure 4. The server 
workstation controls the Geocoding System and is the physical storage 
for the DEM and GCP database, as well as the I/0 interface to SAR data 


stored on optical disks. 
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Fig. 4. Hardware Configuration of the Geocoding System 


To cope with the variety of products, the expected data throughput and 
to perform manual quality control and experimental processing without 
interferring the operational throughput, three indentically equipped 
workstations are foreseen for the system. As stated earlier, the 
geocoded product without terrain correction will mainly be produced by 
the SAR processor. This processer is also based on the workstation 
concept but equipped with special hardware like huge corner turn memory 
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and a fast array processor (Noack, 1987). The Geocoding System 
mentioned above, will produce geocoded products without terrain 
correction for required special map projections and all geocoded 
products which need Digital Elevation Models for precise mapping. 


8. Conclusions 


Aspects for an operational system for geocoding SAR images with and 
without terrain rectification have been discussed and the design for the 
German Processing and Archiving Facility Geocoding System (D-PAF GEOS) 
for ERS-1 SAR images was presented. Though, this system is dedicated to 
ERS-1, it will also cope with SIR-C/X-SAR images and has the capability 
to process other satellite SAR images as well. 


Further, our aim was not only to build up a mission dedicated number 
cruncher, but to define a basis in terms of concepts in software, 
hardware and methodology for advanced user needs to geocoded SAR 
products for the polar platform area. 
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{Text} 1. Abstract 


The accuracy of photogrammetric mapping, as well as the reliability of 
image interpretation, strongly depend upon the quality of the images 
used. It is therefore important to assess or describe image quality. 
Different areas of application and various methods of evaluation imply 
quality demands that differ considerably. An investigation into the 
nature of these demands and into the differences between them is 
necessary as a basis for further investigations of image quality 
criteria. Image quality measures must be relevant with respect to 
subjectively experienced image quality or image usefulness. At the same 
time the measures should be based on objective measurements of image 


properties. 


In this survey, subjective as well as objective image quality criteria 
are studied. Methods for evaluation of image information are 
investigated from an image quality point of view. Special emphasis is 
put on the adjustment of objective assessment methods to subjectively 
experienced image quality with respect to certain applications. 


2. Background 


The results of an application of aerial photographs, either in the form 
of photogrammetric mapping accuracy or in the form of image 
interpretation reliability, is highly dependent on the image quality of 
the photographs used. An example of this i* the accuracy and 
r@liability improvement of photogrammetric mapping from low altitude 
aerial photographs, resulting from the introduction of aerial cameras 
using forward motion compensation (Ackerman 1986, Boberg 1986). Image 
quality affects not only the result of the mapping or interpretation 
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project, but also the time consumed for image evaluation and subsequent 
field completion, and therefore the costs involved in the project. It 
is therefore important, not only to produce aerial photographs of the 
highest possible quality, but also to be able to assess or describe the 


quality of the images produced. 


The main objective of the project, of which this is a status summary, is 
to find methods for assessment and description of the image quality of 
existing aerial photographic negatives. The quality of the image- 
forming process, or of the parts thereof, is thus not stressed. Also, 
the presence of imaged test objects cannot be presupposed. The study is 
restricted to near-vertical photographs, taken with 15/23x23 cm aerial 
mapping cameras on black-and-white film, disregarding, e.g., 
geometrical quality and color balance. In this paper, a number of image 
quality descriptors and assessment methods will be reviewed from the 
point of view of their relevance and practical applicability. The 
review is based upon a literature survey and upon a few preliminary 


experiments. 
3. The concept of image quality 


3.1. Image interpretability 


An imaging system translates an object’s spatial and spectral radiation 
characteristics into a density (or color) distribution over an image 
plane. The imaging process is an information channel, the end product 
of which is a memory for information. Information is retrieved from the 
memory by interpretation of the image. The amount and quality of 
information that can be retrieved, is dependent upon the quality of the 
image, expressed in a relevant way with respect to the human ability to 
interpret the image. Therefore, the interpretability of images will be 
the basis of a discussion of the relevance of different quality 
Gescriptors. It is at this stage important to bear a few things in mind: 


--Image interpretation is a very complex process, the result of which is 
dependent not only of the detail definition of the image and of the 
visual properties of the eye, but also of the psycho-physical 
interactions between the eye and the brain. 


--Image interpretation, 6@6.9g., mapping, is normally performed 
stereoscopically. The mechanism of stereoscopic vision is in this 
respect not fully known. According to Konecny, et al. (1982), features 
could be interpreted stereoscopically with only half the resolution 
required to interpret the same features on single images. Trinder 
(1986) claimed that the boundaries of the objects rather than the areas 
are the important properties used by the visual system to detect the 
stereoscopical disparity. These questions, as well as the effect of 
different image quality for the left and the right eye, should be 


further investigated. 
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--Most interpretation tasks are based upon a priori knowledge or 
expectation of the existence, location or shape of objects. This may 
affect object detectability and recognizability. 


The two main approaches to image quality description are the subjective 
and the objective one. Both are at the same time weak and strong. 
Combining the strengt:: of the methods while avoiding their weaknesses 
has been the more or less explicit aim of many image quality studies. 


ee s Image quality demands 


Aerial photographs are used for many different applications, utilizing 
many different types of instruments. This implies that image quality 
demands may differ considerably. One reason for this is that image 
quality may. mean different things for, and is described in different 
ways by, different categories of image users. For example, stero- 
plotting in analog or analytical plotters, photo interpretation in 
instruments with high image magnification, orthophoto production, image 
coordinate mesurement for block triangulation, image digitizing for 
automatic image matching measurements, and digital mono-plotting from 
image coordinates measured on a digitizing table all may put different 
demands on the quality of the images used. An investigation into the 
nature of these different quality demands is necessary. 


4. The influence of the photographic process 


4.1. Exposure and dodging 


The photographic process has a determining influence on image quality. 
The reason for mentioning this again is the availability of aerial 
mapping cameras with forward motion compensation (FMC), which have 
increased the possibilities to optimize film exposure and made fine- 
grained, low-sensitivity films possible to use also for low altitude 
aerial photography. As reported before (Boberg 1986), image quality as 
assessed by the National Land Survey of Sweden (NLS) has improved using 


this technique. 


Image resolution varies with exposure (Brock 1952, Corten 1960, Zeth 
1984). This restricts the useful log brightness scale and the exposure 
latitude. Exposure should bring also the darkest parts of the subject 
to the straight part of the D/logH curve. Corten used Dmin = Dfog + 
(0.1-0.3), while Graham and Read (1986) indicate best resolution at D =: 
Dfog + (0:4-0.6) and Istomin (1961) at D = Dfog + 0.6-1.0). Gerencser 
(1979) recommends Dfog < 0,2, Dmin = Dfog + (0,2-0,4) and AD =: (0,8- 
1,5). The NLS uses Dmin = 0.5 and Dmax = 1.4, giving a density range of 
OD = 0.9, according to practical experience. 


Istomin (1961) used the amount of exposure range, within which a 


photographic material provides for a pre-determined resolution, as an 
information capacity criterion. 
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Corten (1960) investigated the resolution of transparencies and paper 
prints of different gradation and exposure, as function of the copy 
density. With the aim to optimize the copy resolution, he recommanded 
the use of density compression (dodging) to an apparent density range of 
AD = 0.6 and a high contrast copying material. Since then, the silver 
content of photographic materials has decreased, and air pollution has 
changed the haze conditions. Tone reproduction studies of today’s 


conditions would therefore be of great interest. 


4.2. Photographic tone reproduction 


Tone reproduction diagrams show how densities and contrasts are 
transferred from the object to the positive copy. A limited, 
preliminary measurement and construction of photographic tone 
reproduction curves of aerial photographs of the NLS quality marks 7 
(which is the very best), 5 and 3 (which is barely accepted) has been 
made (see Figure 1). The DlogH curves were constructed with the help of 
the step tablets in the Zeiss Jena LMK camera frame. The solid parts of 
the curves represent the object range. From the graphs, the values in 


Table 1 were derived. 


The following, very preliminary conclusions can be drawn from the 
figures and the table: 


--The negatives are nicely linear, but only a limited part thereof is 
used--apparently the best part from resolution point of view, though! 


--Negative gradation is relatively high for film 2412, which could be 
expected. 


--Negative density difference 4Dn varies with NLS quality mark 
approximately as NLS# = 5 x 4 Dn. 


--Positive gradation is very high for copies from film 2405. As 
automatic dodging has been used, diminishing the positive exposure range 
Alogtp, this is in accordance with Corten’s (1960) recommendations. 


--The gamma product is very high, and inversely proportional to the 
quality mark. 


--The negative density range partly falls below the linear part of the 
(D/logH)p curve, resulting in a simila: result for the final tone 
reproduction curve. This is especially evident in the NLS#5 curve, and 


is probably the result of copy underexposure. 


--An increased exposure and a shorter development time of the copies 
seem preferrable. 
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Summing up, photographic tone reproduction curves seem to be an 
interesting source of information on reasons for image quality 
differences. The effect of dodging on them should be further 
investigated. To use them for image quality assessment, they have to be 


coupled to other quality descriptors. 
5. Subjective image quality descriptors 


5.1. Resolving power 


A resolution test means a subjective decision whether imaged test 
patterns, normally three-bar targets of increasing line frequency, are 
resolvable or not. The test leads to resolving power values in ip/mm 
for the imaging system. The subjectivity of the method is limited to 
the judgment whether the three-bar target is possible to resolve or not. 


Although resolving power measures the combined effect of all the 
information transmitting or degrading parameters in the imaging system, 
it has no absolute or permanent significance, derived from physical 
optics, but simply means the ability of a camera system to reproduce 
fine details of the type of the test object (Brock 1952,1969). 


Resolving power is an uncomplete description of image quality, although 
it is reported to show strong correlation to subjectively experienced 
aerial image quality. Resolution is different for objects of different 
shape, size, arrangement and contrast (Rosenberg 1971, Overington 1974). 
A high-contrast pattern limits the usefulness of the method, as the 
landscape normally imaged in aerial photographs consists of different, 
mainly low contrasts. S6rock (1952) is clearly in favor of an annulus 
target, with target contrast 4D = 0.2. Further, as the method 
presupposes the imaging of man-made targets, it is not applicable for 


routine airphoto production. 
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Thus, resolving power is useful for comparison of the performance of 
different imaging systems, but not for quality description of images in 


hand. 
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5.2. Quality mark assessment 


The National Land Survey of Sweden (the NLS) performs nearly all aerial 
photography for mapping purposes in Sweden. The assessment of image 
quality must not delay the production, and therefore must not be too 
@laborate. For this reason, a relatively simple, subjective image 
quality assessment method is used. Every photographed strip, as a unit, 
is given a subjective quality mark in a seven-step scale, expressing the 
mean photographic image quality of the negatives in the strip. The 
assessment is based upon criteria on illumination, gradation, haze and 
image motion, using a few density readings and ocular examination with a 
magnifying glass. The method was explained more in detail in (Boberg 


1986). 


The method works fairly well, according to experience on how the images 
are received by the users. However, a more thorough investigation into 
the relevance of the quality marks for different categories of image 
users is needed. E.g., Biedermann (1978) has experienced that quality 
judgment depends on object detail and contrast and on the attitude of 
the observer to the imaged subject. Also, the NLS method should be 
Calibrated in a more objective way than via a yearly review of the 
routines by the personnel responsible for the assessment. Measurements 
of objective parameters and calculation of quality descriptors of the © 
negatives would thus improve the method. 


6. Objective image quality descriptors 


6.1. Image acutance 


Acutance is a measure of image sharpness, using the gradient of a 
microdensitometer trace across in edge, which has an infinitely steep 
gradient in the cbject. Istomin (1961) reports of the acutance measure 
used by Jones and Higgins, namely the ratio of the RMS mean gradient of 
an edge response curve and the contrast of the edge, expressed as log 
exposure height of the edge. Contrast as information carrier for image 
interpretation is however better expressed as a density difference, QD - 
Dmax-Dmin. If density is used instead of log intensity, and the 
photographic process is considered practically linear, acutance could be 


expressed as 


1 § @p 


2 
Az= {€£ (==) AD 
§ as} ax 


in which the gradient dD/dx might be determined as the mean slope, like 
a gamma value. 
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Acutance 1s a complement to resolving power. The two might even be 
complementary, as high acutance often means low resolution and vice 
versa. The definition of acutance (in logH values) is related to edge 
gradient analysis for modulation transfer function determination. 


6.2. The modulation transfer function 


If an object 1s regarded to be composed of sinusoidal line patterns of 
different frequencies, an analysis of the degradation of the imaged 
object contrast, as a function of the spatial frequencies of the 
composing line patterns, leads to the modulation transfer function 
(MTF). This means, that if the line pattern frequencies are regarded as 
@ measure of the sizes of image details, the MTF will describe the 
reproduction quality of imaged object details. Although this is not 
generally true, the MTF has been shown to be a good starting point for 


Gerivation of image quality measures. 


MTFs could be determined for each step of the imaging process, as well 
as for the process as a whole. In our case, this is of certain 
interest, because the positive copy is the end user product, while the 
image quality assessment should refer to the negative. Knowing the MTF 
of the copying process, it could be reduced from the MTF of the copy, to 


get the function for the negative. 


MTFs of images in hand must be determined from microdensitometer traces 
of natural image details. Normally, Edge Gradient Analysis (EGA) of 
natural sharp edges has been used for this purpose. Gerencser (1972, 
1974, 1976) has gathered some experience from the EGA method, using the 
graphical method of Scott et al (1963) as well as an automated system. 
The accuracy requirements on the microdensitometer’s length measurement 
was found to be very high. Using known contrast panels, the MTF curves 
don’t have to be normalized to 100 percent modulation at O 1 p/mm 
frequency, which makes it possible to evaluate the contrast reduction of 
the atmosphere. This experience is of interest for this study, as an 
image quality descriptor should encompass the whole imaging chain. The 
possibilities to use natural, known contrasts should be investigated. 
Gerencser expected difficulties in finding suitable edges and to define 


the trace length. 


Especially in large scale aerial photography, the definition of the 
imaged object edge has to be considered carefully. Gerencser (1976) as 
well as Ackermann (1986) encountered great problems in finding suitable 
edges. In some cases, line objects for direct measurement of the Line 
Spread Function (LSF), may be easier to find, @.g., in built-up areas 
with road-markings. Figure 4 shows an example of this. Generally, 
image scale must be considered when suitable targets are chosen. 
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Fig. 4. An example of a microdensitometer trace across a white marking 
line on an airport runway. Image scale 1:15,000 


A basic prerequisite for MTF analysis is that the imaging system is 
linear. The photographic emulsion is principally not, however. 
Welander (1968) has put together results that show that within 
reasonable limits of contrast, gamma product, frequencies and negative 
exposure range, the photographic process may be considered linear for 
practical purposes. wHWendeberg (1967) showed that in the absence of 
adjacency effects, the MTF of the negative is independent of contrast 
and exposure within wide ranges. These results are of great importance 
if the total system result is to be assessed and related to subjective 
image quality. The negative quadrant of the tone reproduction curves 
shown above seems to fulfill the linearity well, Dut the positive 


quadrant does not, if copy exposure is insufficient. 


6.2.1. MTF condensation 


A MTF is normally expressed graphically. Moreover, many MTFs are needed 
to fully describe all quality aspects of an image. It is therefore 
natural, that much effort has been devoted to try to condense the MTF 
into a single figure of merit. Computing the logarithm of the 
integrated area beneath the curve, or of a part thereof, has been the 


most commonly used method. 


A Threshold Modulation (TM) curve for a specific emulsion and 
development describes the image modulation needed to resolve an object, 
as a function of object spatial frequency. The intersection of the 
appropriate MYF and TM curves gives a frequency threshold, which is a 
measure of the system resolving power. This threshold is often used as 
upper limit of the MTF condensation integration. The method is suited 
for prediction of the performance of imaging systems (Rosenbruch 1986). 


A more promising form to use with images in hand is b2sed upon 
integration of a frequency weighted MTF curve, as described below. 
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6.3. Granularity 


Image interpretability is strongly influenced by graininess, especially 
when detail contrast is low. This could be expressed as a signal-to- 
noise ratio (S/N), where the object image contrast is the signal, and 


the graininess is the noise. According to experience, S/N = 4 is needed 
for object detection, and S/N = 10 for recognition (Welander et al 


1972). 


Graininess is measured as granularity, defined as the standard deviation 
oD of a large number of density readings, each covering a very small 
area of the image. Typically, a 24um diameter aperture is used, but 
Gerencser (1972) used a 1,6x200um microdensitometer slit. The Kodak 
method for granularity measurement, leading to RMS granularity 

(10° op) ), 1S supposed to indicate the impression of graininess 
perceived at 12x magnification of the sample. The Selwyn granularity 
measure G is standardized to a unit measuring area a as 


G = oD /2a 


Granularity increases with the square root of the density value, and 
roughly in proportion to the contrast of the printing material (Welander 
1968). This reduces its value as quality measure for images at hand. 


Granularity could be investigated with the help of Fourier analysis, in 
which noise strength is expressed as a function of its spatial 
frequency. This form of description, a Wiener spectrum, is closely 
related to the MTF. Trinder (1980) has investigated this in detail, and 
found it important to consider the S/N in attempting to relate visual 


performance to image quality parameters. 
7. Relations between objective and subjective criteria 


Attempts have been made to correlate physically measurable figures of 
merit to subjectively perceived quality of photographic images. 
Biedermann (1967, 1978) tried in his thesis to relate subjective 
assessments of image quality via weighting functions to measurements of 
modulation transfer functions of the images. Biedermann’s study 
concerned mainly “normal” images like landscapes, portraits, etc., but 
also an a@rial photograph, at a viewing distance of 40cm. He found good 
correlation with subjective quality G (running from 1 “excellent” to 7 
“useless") using the following function of MTF, denoted M(v): 


G = Ci + C2logf 
0 


M(v )B(v)dv 
B(v)dv 


Different weighting functions B8(v) were tested, of the form 


1 
Biv) = TF vs 
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where F could be varied. Thus, image quality was related to the 
logarithm of the area under the weighted MTF. 


Rosenbruch (1986) proposes the agreement on a standard weighting 
function. However, frequency weighted MTFs, or derivatives thereof, as 
descriptors of aerial image quality, require weighting functions adapted 
to the image application and to the image evaluation situation at hand. 
As a base for the choice of frequency weighting functions, a study of 
the importance of different frequencies and contrast distributions for 
different applications is needed. The base for this should be an 
investigation into image quality demands of different categories of 
aerial image users, e.g., stereo operators, image interpreters, 
topographic map makers, forest taxators and physical planners. The 
judgment should be done in their real working situation, with the images 
oriented in evaluating instruments. 


The result of such a study might be a suggested calibration procedure, 
or a replacement method, for the NLS subjective image quality assessment 


method. 
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[Text] ABSTRACT 


The purpose of this paper is to show the present situation of the 
Brazilian SPOT Direct Receiving Station, involving both the SPOT 
Receiving and Recording Station and the SPOT Image Processing Center. 
One of the main functions of the Processing Center is to perform the 
radiometric and geometric treatments in order to eliminate system 
distortions inherent to the HRV-SPOT images generation process. 
Nowadays, the station is capable to process images in the 1A and 18 
levels, in both digital and analogical domains. In addition, new 
developments are being carried for the generation of HRV-SPOT images in 
the levels 2A, 2 and 3, where extra information obtained from ground 
control points and digital terrain elevation models are used to refine 


the geometric system corrections. 


1. INTRODUCTION 


The Brazilian SPOT Direct Receiving Station (8SDRS), set up by INPE, 
begun its operational phase in December, 1987. It comprises the SPOT 
Receiving and Recording Station (SRRS), located at Culaba, Mato 
Grosso State, and the SPOT Image Processing Center (SIPC), situated at 


Cachoeira Paulista, Sac Paulo State. 
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The antenna of SRRS is shared with the Landsat Receiving and Recording 
Station (LRRS), and its location allows full coverage of Brazil and 
most of South America. SIPC makes use of hardware facilities of INPE’s 
Image Generation Department (DGI) which were primarily configured in 
view of TM data processing, and are still used also to this end. 


Figure 1 shows the present hardware configuration botna for SRRS and 
SIPC, while Figure 2 presents the software already implemented. 
Together, they give an overview on the complete SPOT data manipulation 
cycle perfoi-med at BSDRS. 
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Fig. 1. S8SDRS Hardware Chain 


Nowadays, the station is capable of producing images in the 1A and 18 
preprocessing levels, and new efforts are being carried by INPE in order 
to enable the processing at higher levels. These efforts are oriented 
into two directions: 


--Development Project SPOT2, whose objective is to implement a 
processing subsystem corresponding to levels 2A and 2 before the end of 
1988. 


--Research Project SPOT3, where experiments with Digital Terrain 
Elevation Models are being made with the aim of giving subsidies toa 
future development project. 
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The differences between the various levels are concentrated on software 
functions GEO and COR, and so the following sections are dedicated to 
describe the different mathematical models associated with them. 


2. PREPROCESSING LEVELS 1A AND 18 


As stated earlier, these are the levels at which SPIC is presently 
capable to generate SPOT images. Thus, in the scope of the software 
package already implemented, the aim of GEO is to determine the 
geometric correction parameters for the image generation at level 18; 
for all scenes recorded during one pass of the satellite. COR is 
applied for the generation of each scene, both on level 1A and 18; if 
the level is 1A only radiometric corrections are done, while if it is 16 
the geometric corrections are also taken into account, based on the 


parameters computed by GEO. 


2.1. GEO FUNCTION 


On level 18, one-dimensional geometric corrections are used, which 
allows one to write the inverse mapping equations as follows: 


L =f (IT) (1) 
C = 9 (J) (2) 
Aj= hh (I) (3) 


where (L,C) and (I,J) correspond to (line, column) in raw and corrected 
images, respectively. 


The applied geometric correction model, known as photogrammetric model, 
relates each point in a raw image to another one in a tangent plane to 
the reference ellipsoid, with the tangency point corresponding to the 


raw image center. 


Due to the HRV-SPOT sensor characteristics, all the pixels belonging to 
a given raw line are acquired at the same time. Because of this, the 
sight directions relative to these pixels determine a plane, and thus 
the intersection of this plane with the projection plane defines a 
straight line, which validates the hypothesis related to equation 1. 
Nevertheless, it is important to observe that the direction of this line 
will be variable if the heading velocity and the yaw are not constant. 
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Fig. 2. SPIC Software Chain 
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Otherwise, the pixels belonging to a given raw column are sampled at 
different times, being submitted to different effects coming from 
attitude and heading velocity changes and, above all, from Earth’s 
rotation. These effects give rise to displacements between consecutive 
lines, causing discontinuities in the columns; these displacements are 


taken into account by equation 3. 


Functions f, 9g and h were specified as third order polynoimials, being 
left to GEO the determination of their coefficients {fi}, {gi} and {hi}, 
as well as the computation of geodetic coordinates (latitude 6 and 
longitudeA) of the scene center and of its corners. 


The polynomials coefficients are determined by applying the 
photogrammetric model to two sets of 61 breakpoints each, one regularly 
spaced along the center line ({gi}) and the other along the center 


column ({fi} and {hi}). 


The photogrammetric model can be depicted by a two equations system, tne 
first one corresponding to the straight line defined by the satellite 
position and some sight direction, and the other one corresponding to 


the reference surface. 


Initially, the position of the satellite and the reference surface are 
entailed to the Greenwich geocentric system, and the sight angles are 
related to the HRV imager reference system as a function of operation 
mode, mirror position and detector number. Transformation from the 

latter to the former system is done in three steps, through boresight 


and attitude angles, and ephemeris data. 


Attitude angles are determined by integrating attitude rate data 
broadcasted by the satellite. This integration is first carried out 
assuming a null attitude for the first received attitude rate record; 
later, when the geometric parameters for sach scene are to be computed, 
the attitude at scene center time is subtracted from the integrated data 
corresponding to that scene, thus zeroing the scene center attitude. 
Ephemeris data obtained via telex link are interpolated through a 


seventh order Lagrangian polynomial. 


The projection plane has its tangency point to the ellipsoid defined by 
the application of the photogrammetric model over the ellipsoid for the 
scene center. The photogrammetric model is then applied over this plane 
for every breakpoint, leading to a set of points with coordinates 
referred to the Greenwich geocentric system. 


Afterwards, a projection system is defined with origin in the tangency 
point to the ellipsoid, Y axis corresponding to the center line on the 
projection plane, eastward oriented, and X axis orthogonal to Y, 
southward oriented. Using the origin coordinates and the unitary 
vectors X and Y, both relat<d to the Greenwich geocentric system, it is 
possible to perform the transformation from this system to the 
projection system for all the breakpoints. 
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The lines and columns of the corrected image are oriented along the Y 
and X axis, respectively: for this reason, the conversion from metric 
projection coordinates to discrete coordinates (I,J) is easy and 

conditioned only to the nominal spatial resolution and the corrected 


dimensions. 


The coordinates »l,y and J, (J of the center column 
points are then related, and tne shetttoiefis {fi} and {hi} determined. 
In the same way, Ck and Jk of the center line points are related, and 


the coefficients {gi} computed. 


The geodetic coordinates of the scene corners are obtained by applying 
the photogrammetric model over the reference ellipsoid for the four raw 


image corners. 
2.2. COR FUNCTION 


This function 1s applied when an HRV-SPOT scene is to be generated, both 
on level 1A and 18. In the first situation, only a radiometric 
calibration is executed, corresponding to a normalization of detector 
responses in each spectral band. On level 18, besides the calibration, 
also radiometric deconvolutions along iines and columns and geometric 


corrections are applied. 


The gray level interpolation associated to the geometric corrections may 
be perfomed in two different ways: nearest neighbour for lines and 
columns, or nearest neighbour for lines and cubic convolution for 
columns. In the former option, polynomials f and g are used to 
determine respectively the lines and columns to be suppressed or 
repeated. In the latter, polynomial f is used in the same was as in the 
former option, and polynomial g defines, within a given line, the four 
neighbours to be considered in the convolution. In both options, 
polynomial h is used to compute the number of black pixels corresponding 


to the corrected image left skew. 


All operations executed by COR require the raw image to be read only 
once. 


3. PREPROCESSING LEVELS 2A AND 28 


The HRV-SPOT imagary generation subsystem for levels 2A and 2 uses a 
geometric correction model which is applied in three steps. 


In the first one, a direct mapping is determined through the 
photogrammetric model over the ellipsoid for a grid of 61 x 6l 
breakpoints regularly distributed over the raw image. A set 

is obtained for each breakpoint, and the geodetic coordinated are 
transformed into projection coordinates (UTM, SOM, etc.), which are then 
converted into discrete image coordinates (I,J). As a consequence, a 
direct relation (L,C)-(I,J) is set up. 











The second step begins with the direct mapping inversion, by using a 
grid of breakpoints regularly distributed over the corrected image. For 
this purpose, it is current to use from third up to fifth order 













polynomials, which are defined by the resulting grid in the first step. 
This model has some disadvantages, as polynomials cannot represent high 
frequency distortions. Moreover, they are not reliable when applied to 
points located outside of the original grid. INPE is now searching into 
the determination of the inverse grid breakpoints from the four direct 

grid nearest ones. 









The last step refers to the resampling process, where a gray level is 
determined for every point on the corrected image. The resampling 
includes a rotation over the images in order to orientate them 
northwards. Every line on the corrected image is related to 950 lines 
on the raw image in the P mode or 475 lines in the XS mode, so 
increasing the necessary computer memory, INPE intends to perform this 
resampling through the algorithms presented by Friedmann (1981) and by 
Craig and Green (1987). 




























The procedures above described are enough to produce a level 2A 
corrected image. In the case of level 2, ground control points are 
necessary to improve some geometric correction parameters, aiming at the 
use of the images as up-to-date planimetric maps. INPE has been 
concentrating its efforts in three possibilities: a polynomial model 
(Bahr, 1978), a Kalman filtering attitude correction model (Caron & 
Simon, 1975), and an orbit correction model, comprising an ephemeris 
data correction and some well-known numerical orbit integrators, with 
the purpose of producing a fully corrected image over a scene without 
ground control points (d’Alge et alli, 1988). 


4. LEVEL 3 PERSPECTIVES 


The SPOT satellite has an off-nadir viewing capability and, as a 
consequence, a stereo viewing possibility. It is feasible to achieve a 
base-to-he:eht ratio of around, being possible to use the images for 
topographic mapping. Moreover, the off-nadir viewing allows the same 
area to be observed on two or five successive days, depending on the 
required latitudy. It also permits the bidirectional reflectance 

di’ Sution function to be investigated, and covers an area which is 

l than in the near vertical viewing case. 


Tnis preprocessing geometric level rectifies the off-nadir viewing 
images within a geodetic accuracy better than 50 m. The merging between 
an off-nadir viewing image and another image over the same scene shall 
occur within a half-pixel accuracy. The planimetric displacements due 
to the existing relief are corrected and the resulting images are 
northward oriented. 











The relief displacement of a point on the Earth’s surface depends on its 
height and on the HRV viewing angle. Therefore, for height of about 
1.000m and a viewing angle of 272, the planimetric displacement due to 
the relief reaches 500 m (Guichard, 1986). 


The current methodology foresees the application of the same 
photogrammetric approach used for level 2. In the off-nadir viewing 
case, any of the sight directions intersects two ellipsoids, both 
related to the reference ellipsoid: one is defined by the maximum 
height, and the other by the minimum height over the scene to be 
corrected. As a consequence, a point on the corrected image is related 
to two different points on the raw image, and its true position is 
defined by some interpolation algorithm, considering the information of 
a digital terrain elevation model. 


The researches about this geometric correction level are being carried 
at INPE’s Geographic Information System, comprising four data sets: raw 
image, grid related to the maximum height ellipsoid, grid related to the 
minimum height ellipsoid, and digital terrain elevation model. 


S. CONCLUSION 


This year, BSDRS will receive and record 4,000 SPOT scenes. It is 
capable to process them only in levels 1A and 18. For the next year 
BSDRS is expected to generate images also in levels 2A, 2 and 3. 
Nowadays, the software system of the SIPC is being adapted to process 
these levels. The first results of level 18 geometric evaluation were 
satisfactory, leading to an intrinsic error of about 30 meters and an 
extrinsic error of around 500 amd 300 meters, in the across-track and 


along-track directions, respectively. 
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{Text}) Abstract 


Japanese Earth Resources Satellite-1(JERS-1) will be launched from 
Tanegashima Space Center, NASDA by H-I rocket in February 1992. The 
JERS-1 has two observation equipments: synthetic aperture radar (SAR) 
and optical sensor (OPS) which consists of visible, near-infrared and 
short-wave infrared radiometer. In order to evaluate parameters of SAR, 
to develop processing and calibration method and to evaluate the utility 
of SAR data in various fields of applications, NASDA conducted SAR-S80 
experiment by using variable parameter SAR mounted on Convair 580 owned 
by Canada Centre for Remote Sensing (CCRS) during October and November , 
1983 in collaboration with many research organizations. As a part of 
this experiment, NASDA conducted calibration experiment in Tsukuba where 
several corner reflectors were set up. NASDA processed about 100 scenes 
of SAR-S80 by newly developed autofocusing method and evaluated SAR data 
from various points of view. Moreover, Communications Research 
Laboratory (CRL) participated in SIR-B experiment planned and conducted 
by JPL in collaboration with several agencies including NASDA. 


In this paper, outline of data processing and calibration for SAR data 
in Japan are presented. 


1. Introduction 


Japanese Earth Resources Satellite-1 (JERS-1) is scheduled for launch in 
February 1992 and will be launched by H-I rocket from Tanegashima Space 
Center, NASDA. Space segment of the JERS-1 is jointly developed by 
NASDA and JAROS and ground segment is developed by Earth Observation 
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Center (EOC), NASDA. The JERS-1 has two observation equipments: 
synthetic aperture radar (SAR) and optical sensor (OPS). Main mission 
objectives of JERS-1 are to evaluate newly developed sensors and 
spacecraft and to collect data of the earth concerning resources 
exploration, geology, agriculture, forestry, land use, prevention of 
Gisaster, coastal monitoring, sea ice monitoring and others. 


In order to evaluate parameters of SAR, to develop processing and 
calibration metnods, to evaluate utility of SAR data in various 
application fields and to study future advanced SAR, NASDA conducted 
SAR-S80O experiment by using variable parameter SAR mounted on Convair 
580 owned by Canada Centre for Remote Sensing (CCRS) in October and 
November 1983 in collaboration with many research organizations. Asa 
part of SAR-S80 experiment NASDA conducted calibration experiment in 
Tsukuba Science City where several corner reflectors were set up. NASDA 
processed about 100 scenes of SAR-S80 data and provided SAR-S80 data to 
principal investigators who participated in SAR-S80 data evaluation 
program in 1984-1986. 


In this paper, data processing and calibration for SAR data in Japan are 
presented by focusing on SAR-S80 data and JERS-1 data. 


2. Outline of JERS-1 SAR data 
2.1. Outline of JERS-1 SAR 


Bread board model (88M) of mission instruments of JERS-1 was developed 
by NASDA. Engineering model (EM) of mission instruments of JERS-1 is 
developed by Resources Remote Sensing System (RRSS) and Proto-flight 
model (PFM) will be developed by Japan Resources Observation System 
Organization (JAROS). On the other hand, Bus instruments of JERS-1 are 
developed by NASDA. 


Table 1 shows the specifications of JERS-1 SAR. L band was selected by 
considering weight and power of payload and penetration through 
vegetation. It was confirmed in SAR-S80 experiment and other 
experiments that L band penetrates vegetation such as trees deeper than 
C and X band. This is desirable for resources exploration. 


2.2. JERS-1 data type and processing level 


Ground segments of JERS-1 are developed by EOC, NASDA. Ground system 
for JERS-1 is under conceptual design in 1988. Basic design will begin 
in 1989. In order to develop processing system, it is necessary to 
define processing data volume, data type and processing level for JERS-1 
data. In this case, both SAR data and OPS data should be considered. 
Table 2 shows tentative processing level of JERS-1 SAR data. Table 3 
shows tentative data type of JERS-1 SAR data. Since the velocity of 
SAR is about 7km/sec and PRF is 1505.8-1606.0, sampling interval is 
about 4.5m. So, pixel size of 1 and 3 looks are defined as 3.125m and 
12.5m, respectively. Spatial resolution of SAR is 18m (3 looks) and 6m 
(1 look). Though there are many window functions, typical window 











functions (Square, Riesz and Hamming) can be selected. Concerning 
output data type, both amplitude and logarithmic expression can be 
selected by considering dynamic range. Since the logarithmic 
expression provides quasi-normal distribution, this 1s desirable in 
case of land use classification by combining with optical sensor data 
and using maximum likelihood method. 


Table 1. Specifications of JERS-1 SAR 


Altitude 568 km 
Transmit Frequency L band (1.275 GHz) 
Polarization H-H 
Swath Width 75 km 
Off-nadir Angle 35 degree 
Incident Angle 38.7 degree (Center of Swath Width) 
Transmitted Pulse 

Type linear down chirp 

Bandwidth 1S MHz 

Initial Frequency 1282.5 MHz 

Pulse Length 35 sec 

Sampling Rate 17.1 MHz 

Pulse Repetition 

Frequency (PRF ) 1505-1605 in 25Hz incriments 

(1508.8, 1530.1, 1555.2, 1581.1, 1606.0) 

Noise Equivalent -20.5 dB 
Signal to Ambiguity >14 ab 
(S/A) 
Raw Data Bit Precision 31/30 
Raw Data Rate 30x2 Mbps 


Table 2. Tentative processing level of JERS-1SAR data 


Level 0: Raw data 
Level 1: Basic image (Range and Azimuth Compression) 
Level 2: System corrected (Bulk) image 
Radiometrically and geometrically corrected image) 
Level 3: Precisely corrected image by using GCP. 
Level 4: Precisely corrected image by using Digital Terrain 
Model (DTM) 


7 
Ma 
4, 


Remarks; Main product is level 2. 
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Table 3. Tentative JERS-1 SAR data processing types 


Map projection method: UTM (Universal Transverse Mercator) 
SOM (Space Oblique Mercator ) 
PS (Polar Stereo) 


Resampling method: Nearest Neighbor (NN) 

Cubic Convolution (CC) 
Frame Size 75km x 75km (level O: 7Skm x 9Okm) 
Pixel Size 3.125m (1 look), 12.5m (3 looks) 
Window Function Square, Reisz, Hamming 
Output Data Type Amplitude/Logarithmic expression 


Recording Format (CCT) Super structure 
(Based upon CEOS WGD SAR format 


standardization) 
Recording density(CCT) 6250 BPI, 1600 BPI 
Number of Tracks (CCT) 7 


3. Data processing flow of SAR data 


3.1. SAR-580 data 


NASDA collected SAR data over 8 test sites (Tsukuba, Hiratsuka, Off- 
shore of Hiratsuka, Mt. Fuji, Oshima, Miyakejima, Akita (Futatsui) and 
Shizukuishi/Kosaka) in 28 October-4 November 1983 by using Convair 580 
with variable parameter SAR which are chartered from CCRS based upon the 
arrangement between CCRS and NASDA which were concluded on 23 August 


1983. 


NASDA processed about 100 scene of SAR data in 1984 by newly developed 
auto-focusing method. Figure 1 shows the data processing flow of SAR- 
580 data. There are two methods for auto-focusing: (1) to employ the 
relationship between target positional shift and processing parameter 
error, (2) to employ the relationshp between image contrast in the 
processed image and processing parameter error. The first method was 
used for our program. Main feature of this program is to employ auto- 
focusing method both in range and in azimuth direction. The feed back 
loop of auto-focusing converges within 2-3 repetitions. 


S. 2. JERS-1 SAR data 


Since JERS-1 has mission data recorder (MDR), it is possible to collect 
SAR data and OPS data all over the world excluding some of arctic and 
antarctic region (SAR:86.015 degrees N-90 degrees N, 78.661 degrees S-90 
degrees south, O0OPS:82.338 degrees N-90 N, 82.338 degrees S-90 degrees S. 
If only land and coastal region is observed, about 30000 scenes are 
obtained. JERS-1 data will be received by Earth Observation Center 
(EOC) in Hatoyama, by Fairbanks ground station and Gatineau ground 
Station in Canada. Data will be dumped at Fairbanks ground station. 
There may be other foreigr ground stations to receive JERS-1 data within 
limited power and stored command number for spacecraft. 
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Minimum requirement of observation by JERS-1 is to observe land and 
coastal area in the world once during two years for SAR and twice during 


two years for OPS. 


Figure 2 shows the data processing flow for JERS-1 data. In order to 
get high quality image, it 1s necessary to introduce auto-focusing 
method even to spaceborne SAR data processing. In case of JERS-1 SAR 
data processing, this auto-focusing method will be applied to azimuth 
direction by using relationship between target positional shift and 
processing parameter error. In principle, it is impossible to apply 
this method to very flat area such as calm ocean. 


Processing level of JERS-1 SAR data will be explained as follows: 


(a) Level O: Raw data 


Radar hologram data including data and parameters required for 
processing such as orbital data, attitude data, on-board calibration 


data, PRF, Sensitive Time Control (STC). 


(bo) Level 1: Basic image 


Range and azimuth compressed image by using FFT accompanied by look 
processing (1 look and 3 looks) and auto-focusing. 1 look image is 
superior to 3 looks image in obtaining higher apesses resolution and 


detailed information. 
(c) Level 2: System corrected image (Bulk image) 


Geometrically and radiometrically corrected image. Geometric correction 
includes conversion from slant range coordinates to ground range 
coordinates, resampling. Radiometric correction includes STC 
compensation, antenna pattern compensation, incidence angle 
compensation, receiver gain compensation. (Position error: + over -lkm) 


(d) Level 3: Precisely corrected image 


Precisely corrected by using GCP. By using GCP and map coordinates, 
antenna position is estimated and by using orbital data and attitude 
data precisely corrected image is obtained. (Position error: spatial 


resolution) 
(e) Level 4: Geocoded image 


Position error due to terrain relief is eliminated by using digital 
terrain model (DTM). 
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Fig 1. Data Processing Flow of SAR-580 Data 
Note: 1. JERS-1 SAR simulation software is included. 


2. Geometric correction (A): The multi~-look image is converted from 
Slant range coordinates to ground range coordinates and resampled to 


adjust image pixel spacing. 


3. Platform location estimation: Using the GCP coordinates of the 
multi-look image and the map, platform position is estimated as an 
approximate function. 


4. Geometric and radiometric correction (8): The results of the 
Platform location estimation are used to correct for geometric 
distortion in the image data. The data is projected onto the map 
coordinates and radiometric distortion is corrected. 


5. Geometric correction (C): Geometric distortion is corrected for 
terrain relief, using digital terrain model. 
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Fig. 2 Data Processing Flow of JERS-1 SAR data 
4. Calibration flow of SAR data 


4.1. SAR-S8O data 
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Table 4. Type of Corner Reflection 


Type Size ! Number Testsite 
INTERA corner reflector 1725mm ukuba 
NASDA corner reflector A 1049mm : — Science City 
ditto B 822mm ‘« ditto 
ditto C 502mm 1 ditto 
eeee D 376mm 1 ditto 
, . ee E 14lmm 1 only for test 
NTERA ector 287mm 3 Ts Scien 
ERSDAC corner reflector 714mm 5 Oshima oo Ty 
rpm reflector 
al corner reflector) 150mm/ + 
ow Tsukuba Science City 
! pa 
4 2 
2 b 
6." 4ka 4(3A) 6.= 8Xa ory 


As a part of SAR-580 experiment, NASDA conducted calibration experiment 
in Tsukuba Science City on 28 October and 4 November 1985. NASDA 
developed several corner reflectors. Size 1s selected depending upon 
wavelength for L band, C band and X band. Figure 3 shows one of corner 
reflectors. Figure 4 shows example of the dependency of back scattering 
coefficient on incidence angle by using calibrated SAR-S580 data. 
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Fig. 3 Example of corner reflector Fig. 4 Dependency of on incidence 
in case of field in 


Tsukuba 


4.2. JERS-1 SAR data 
Calibration of JERS-1 SAR data is essential in various fields of 


application such as land use, monitoring of vegetation. Figure 5 shows 
calibration flow of JERS-1 SAR data. 
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Fig. 5 Calibration Flow of JERS-1 SAR data 


In Japan, there are very few testsites for space borne SAR calibration 
because very few flat and uniform areas exist. Therefore, it is 
necessary to develop special method for calibration which can be 


utilized in Japan. 


If transponder set near lake or pond has delay time such as 5-6 sec 
between receiving and transmitting, point target appear in the lake or 
the pond. Since back scattering coefficient of the lake or pond is 
sufficiently small, it is easy to calibrate SAR data. 


S. Evaluation program 


Before SAR-580 experiment, committee was held 4 times and flight course 
and content of experiment were studied. As part of SAR-S80O experiment, 
evaluation program was made and processed SAR-580 data were provided to 
princical investigators who participated in this program. All fruitful 
results were presented at SAR-580 experiment symposium in Tokyo on 25 


September 1986. 


The evaluation theme of SAR-580 data is as follows: 


(a) Joint Research Organization 


National Research Center for On detection of ocean waves 
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Disaster Prevention 


National Institute for En- 
vironmental Studies 
Forestry & Forest Product 


Tokai University Research 
and Information Center 


Earth Resources Satellite Data 
Analysis Center (ERSDAC) 


(b) Cooperative organization 


Institute of Space and 
Astronautical Science (ISAS) 


Meteorological Research Institute 


Geographical Survey Institute 


Department of Computer 
Science, Iwate University 


Science Information Processing 
Center, University of Tsukuba 
Remote Sensing & Image Research 
Center of Chiba University 
Kanazawa Institute of Technology 


College of Engineering, 
Hosei University 


Toba National College of 
Maritime Technology 
Remote Sensing Technology 
Center of Japan (RESTEC) 
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Analysis of urban land use 
utilizing SAR-S80 image 
Possibility for interpretation 
on landslide landform 

Soil moisture measurement using 
SAR-580 data 

Evaluation of SAR-580 data on 
forest type classification 

An optical image display 

method and effectivity for land 
cover classification using SAR 
data 

Geologic analysis of SAR-580 
data 


Studies on microwave signature 


Investigation on reflectivities 
of radar antenna reflection 
measured by SAR 

SAR imageries obtained in the 
sea areas around the Oshima and 
Miyake islands 

Inspection of ground surface 


state of Oshima and Miyakejima 


volcanoes with SAR-S80 images 
Feasibility study of SAR image 
interpretation for planimetric 
objects 

Investigation of corrsepondency 
between Images 1n SAR-S80 data 
and characteristics of land 
cover objects 

Image processing of synthetic 
aperture radar (SAR-580) data 
Evaluation of SAR-S5S80 imagery 


The classification of the ground 
cover types on the rugged 
terrain by the SAR-S80 data 
Fundamental study for verifica- 
tion of classifying SAR images 
by checking the objects in field 
Sea surface phenomena investi- 
gation using the SAR data 
Detection of terrain height 

from stereo SAR image and 
optically processed images 
Accuracy of land cover classi- 
fication by SAR image 








CRL participated in SIR-B experiments planned and conducted by JPL in 
collaboration with several agencies including NASDA in October 1984 and 
lt was found that a simulated o11 slick area over sea was clearly 


detected (Ref.8). 
6. Conclusion 


EOC, NASDA is developing JERS-1 ground system (Conceptual design). 
Processing level and format was studied in the committee on JERS-1 
ground system. Considering user requirement and capability of 
spacecraft, basic and detailed design will begin from 1989. 
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{Text} Abstract 


The paper investigates the capability of object reconstruction from 
space for two specific projects. Both of them are based on the use of 
three line scanner systems allowing for simultaneous digital image 
recording from three different angles. An efficient model for combined 
point determination is proposed, considering the digital image data as 
well as orbit information, GPS and INS observations. The resulting 
accuracy 1S estimated from project oriented simulations. 


Finally the possibilities for DTM generation, orthophoto production ana 
stereo compilation from three line scanner imagery are discussed. 


1. Introduction and Review 


A line scanner or opto-electronic camera allows for digital image 
recording line by line according to the push-broom principle. A three 
line scanner records the object simultaneously in forward, downward and 
backward looking direction. The three-fold coverage makes a rigorous 
reconstruction of the object possible (Hofmann et al., 1982). 


Three line scanner systems shall be used from space in conjunction with 
the MEOSS experiment (Lanzl, 1986) AND THE MOMS-02/D2 project (Ackermann 


et al., 1988). 
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The paper first describes the proposed mathematical model for combined 
point determination using three line scanner imagery and general contro! 


information. 


Then accuracy estimations for the space projects MEOSS and MOMS-02/D2 
are performed by means of computer simulations. 


Finally the possibilities for DTM generation, orthophoto production and 
stereo compilation from three line scanner imagery are discussed. 


2. Combined Point Determination Using Three Line Scanner Imagery 


A three line scanner contains in the focal plane of the objective three 
linear CCD-sensors, which are oriented perpendicularly to the direction 
of flgiht. During the flight the shutter of the objective 1s open and 
the sensors are read out with constant frequency. Within each read 
cycle a digital image is recorded, consisting of only three lines. If 
read cycle frequency, speac of flight, image scale and size of sensor 
elements are in proper relation to each cther the terrain is scanned 
with nearly square shaped object surface elements (Figure 1). Generally 
each object point 1s imaged three times, first by the forward looking 
sensor a, then by the downward looking sensor b and finally dy the 
backward looking sensor c. So, classical frame photography with proper 
longitudinal overlap is replaced by a composition of a large number of 
successive three line images. Basically, each of them has its own 


exterior orientation. 


Consequently the mathematical model for the reconstruction of the 
exterior orientation should use 6 unknown parameters for each three line 
image Ij . In practice, however, there is not enough information to 
determine such a large number of unknowns. A solution to this problem 
can be found by the introduction of a proper functional model for the 
variation of the exterior orientation parameters. In the simplest case 
piecewise linear functions can be used. Then so-called orientation 
images |, are defined for which *%e parameters of the exterior 
orientation are unknown and the 6 parameters of an arbitrary image 

are represented by linear functions of the 12 unknown parameters of the 
two nearest orientation images I, and Ie . For details see (Ebner 
and Muller, 1986). 
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Fig. 1 Image Recording Using a Three Line Scanner 


In that way point determination using three line scanner imagery can be 
performed similar to classical photogrammetric point determination. 
Based on the bundle method collinearity equations are formulated for tne 
image coordinates Xi4 of the object points Py in the images Ij 
The exterior srcanibiel parameters Ry 94 a5 64 43 ey) , Nowever. 
are replaced by the unknown parameters! 

Bx Sie 3x Re Ker Picea 2k+1 Sk+r bk+l Ek+2 

of the neighbouring efiatiatial images 1, “and k+1° 





All relevant image coordinates are considered as 
observations ina least squares » SR Further ~- observation 
equations are formulated for the control point coordinates. Unknowns of 
the adjustment are the coordinates of the object points P; and the 
parameters of the orientation images I, - This concept can easily be 
generalized to combined point determination, if instead of or in 
addition to control points general control information is available, 
e.g., direct geodetic observations, measured orientation parameters or 


DTM data. 
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In conjunction with the planned space project MOMS-02,;02, GPS and INS 
observations are of particular interest. Wath the MEOSS experiment 
where such data are not avallable orbit information can be used. 


a Accuracy Estimations for che Space Projects MEOSS and MOMS-02/D2 


is wie MEOSS Experiment 


The Monocular Electro-Optical Stereo Scanner (MEOSS), whaiach has been 
designed and manufactured by the Institute of Optoelectronics at the 
‘Deutsche Forschungs’ und Versuchsanstalt fur Luftund Raumfanrt e.V. ° 
(DFVLR) will be the first three line camera in space (Lanzl, 1986). Tne 
mission is planned for 1988 on the SROSS satellite provided by the 
Indian Space Research Organisation (ISRO). 


Besides thematic applications relating to cartography, meteorcloay and 
geosciences, MEOSS offers tha first possibility of three dimensional 
photogrammetric object reconstruction based on spaceborne three line 
imagery. Although the ground resolution is limited to 52m by 79 m, th2 
MEOSS image data may be used for basic investigations and tests of 
photogrammetric compilation methods for threefold stereo scanner data. 


The intention of this paper 1s to give an a priori estimation of the 
accuracy of point determination using MEOSS imagery. By means of 
computer simulations with project oriented camera- and flight 


parameters, the on. ie standard deviations Oxi» oyi? C2 of the 
unknown coordinates Xi» i> 24 of object points have been 
calcuated. In (Hofmann and Muller, 1988) the principle of this method 


1s described in detail. 


The simulations point out: 


--the advantages of a simultaneous adjustment of two overlapping strips 
versus single strip adjustment and 


--the possibility of the integration of orbit models, whaich provide 
information about the satellite position with limited accuracy. 


The first item has to be seen 1N conjunction with the important feature 
of the MEOSS misssion, that areas of interest can be recorded twice from 
different orbits, providing sixfold overlap. The advantages of such a 
block adjustment of three line imagery have been demonstrated in (Ebner 


and Muller, 1987). 


Figure 2 shows an example for 5 consecutive orbits in the course of the 
MEOSS experiment representing the coverage of one day. 


The rastered areas have been arbitrarily selected to investigate the 
accuracy of the object reconstruction from MEOSS imagery: 


--area A represents a single strip, 
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--area 8 represents a block, consisting of two strips, which intersect 
at an angle of approximately 22 grad, 


--area C represents a block, consisting of two strips, which intersect 
at an angle of approximately 61 grad. 


The characteristic data for the MEOSS simulations are listed in Table l. 


Concerning the chosen distance between orientation images, a remark is 
necessary: 16 km 1s probably a too pessimistic value, because the 
Variations of the parameters of the exterior orientation might be 
described using longer intervals without increasing the interpolation 


errors. 























Fig. 2: Example for 5 consecutive orbits in the course of the MEOSS 
experiment (taken from Lanzl, 1986) 
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camera parameters 
calibrated focal senge 
spacing forward/vertical sensor 

and vertical/backward sensor 
standard deviation of image coordinates do 











flying height 

baselength between forward/vertical look 
and between vertical/backward look 

strip length 

strip width 

distance between orientation images 


ebiect point distribution 
distance in flight direction 
Gistance across flight direction 
height of all object points 





soy gyey e BEE 





Table 1. Characteristic Data for the MEOSS Simulations 
The following simulations have been performed: 


--strip / block adjustment for A, 8 and C assuming error-free parameters 
of the exterior orientation (accuracy limits), 


--strip / block adjustment for A, B and C using four error-free ground 
control points (GCP) for each strip, 


--strip / block adjustment for A, B and C assuming given position 
parameters for the orientation images ' with Oxk - Svk = Os; = 
100 m (orbit information) in addition to the GCPs- 7 


Table 2 shows the results of the simulations. Given are the rms values 
BQ» Ad» ug of the theoretical etendersg deviations OR» OFi» Od, 
for the adjusted object point coordinates i> 94> 2,. 
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Table 2. Results of the MEOSS Simulations 


accuracy limits 4 GCPs/strip 4 GCPs/strip 
given pos.paran. 








(m) |ue(m)Juetm) |wolm) |ue(m)fuetm) |uotm) |uetm)} 


152 268 309 





30 30 91 





35 | 35 | 109. 





























The following conclusions can be drawn from Table 2. 


--The single strip shows poor geometry. The accuracy could be improved 
by using a greater distance between the orientation images, if the 
satellite dynamics allows for this. 


--The use of two overlapping strips in a simultaneous block adjustment 
improves the accuracy of point determination drastically. Especially 
the achieved planimetric accuracy is close to the accuracy limit. 


--The introduction of given position parameters (orbit information) 
further improves the accuracy, particularly in height. 


3.2. MOMS-02/D2 Project 


The MOMS-02/D2 project, which is also directed by the DFVLR, 1s based on 
a panchromatic three line scanner with high ground resolution in 
connection with multispectral channels (Ackermann et al., 1988). The 
use of the system is planned for 1991 during the German D2 mission on a 


NASA Space Shuttle flight. 
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One of the main tasks of the MOMS-02/D2 project is the rigorous and 
precise geometric evaluation of threefold stereo data. To secure, that 
the high requirements of this mission concerning the geometry of object 
reconstruction are met, comprehensive studies have to be performed. 


The simulation results given in this paper show a general trend, but 
should not be interpreted as final. 


There is one characteristic feature of the project, which changes the 
geometry of photogrammetric point determination for the worse: the 
ratio of flying height and strip width is approximately 9:1. 


Improvement of the geometry may be obtained by using information from 
navigation systems for the position and attitude parameters. 


The following Simulations have been performed: 


--adjustment of a strip assuming error-free parameters of the exterior 
orientation (accuracy limits), 


--adjustment of a strip using four error-free ground control points, 


--adjustment of a strip assuming measured position and/or attitude 
parameters with different accuracies additionally to four error-free 
ground control points. 


In practice attitude and position information can be derived from 
Inertial Navigation Systems (INS) and/or from the Global Positioning 
System (GPS). 


Table 3 lists the characteristic data for the MOMS-02/D2 simulations. 





In Table 4 the results of the investigations are summarized. Again. the 


rms values MR, w¢, Ue of the theoretical standard deviations 
Oxi» go OF for the adjusted object point coordinates 
Ri, § $, are listed. 


In this case, nowever, not all points have been used to calculate the 
rms values, but only the ones, which are located in the threefoldly 
covered area of the strip. ; 


The obtained results may be analysed as follows 


--Compared to the accuracy limits (Sim. 11) the accuracy of the strip 
adjustment based on only 4 control points (Sim. 1) is poor. 


--The use of measured attitude parameters (Sim. 2 to Sim. 4) improves 
the height accuracy to some extent. There is no difference, whether 
1.0, 0.5 or 0.2 mgrad is assumed. 











--Measured position parameters (Sim. § to Sim. 7) yield smaller values 
for pe and pp, but there are still significant differences to the 
accuracy limits. 


--Really satisfactory results can only be obtained by the simultaneous 
use of accurate attitude and position data (Sim. 9 and Sim. 10). 


Table 3. Characteristic Data for the MOMS-02/D2 Simulations 


camera parameters 
calibrated focal length 660.00 mm 
effective spacing forward/vertical sensor 
and vertical/backward sensor 300.42 mm 

standard deviation of image coordinates 09 5.0 um 
flight parameters 
flying height 334 im 
baselength between forward/vertical look 

and between vertical/backward look 152 ka 
strip length 606 kn 
strip width 36 kn 
distance between orientation images 16 ka 
object point distribution 
distance in flight direction 2 kn 
distance across flight direction 9 kn 
height of all object points 0 km 








The described studies will be continued. The final goal 1s a 
comprehensive analysis of the accuracy properties of combined point 
determination within the MOMS-02/D2 project. 
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Table 4. Results of the MOMS-02/D2 Simulations 


estimated accuracy 


attitude data position data 


















































“ ‘tmgtaa) ue (m)} wp (m)] we (m) 
ra none 21.8 8.0 64.1 
2 1.0 none 21.2 6.7 35.7 
| 3 0.5 none 21.2 6.6 35.7 
4| 0.2 none 21.2 6.6 35.7 
5 none 25.0 11.1 7.4 24.3 
. none 5.0 7.1 7.3 18.1 
7 none 1.0 6.7 7.3 16.6 
~ 1.0 25.0 6.7 4.4 14.1 
9 0.5 5.0 2.6 2.1 6.2 
10 0.2 1.0 1.6 1.6 4.2 
fal error-tree orient, parameters | is | a3 | 39 


4. Steps of Object Reconstruction and Conclusion 


The process of object reconstruction based on digital imagery of three 
line scanners includes the following tasks. 


--Image recording analogous to the push-broom principle (see Chapter 2). 


--Determination of homologous points by a human operator or by methods 
of digital image matching. 


--Preprocessing of additional control information, e.g., from INS or 
GPS. 


--The main task consists in the reconstruction of the exterior 
orientation based on the mathematical model described in Chapter 2. 


--After the reconstruction of the exterior orientation, DIM image points 
can be determined with methods of digital image matching, yielding DTM 
object points. Then the DTM surface can be mathematically described by 
proper models. 
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--Using the image data of line a, 6b or ¢ orthophotos and central 
perspectives may be derived by means of digital image transformation 


methods. 


--If an analogue film is produced from the original or transformed 
Gigital data, analytical plotters can be used for the evaluation of 


three line imagery. 





--Finally digital mono plotting (compilation of the images generated by 
sensor a. bor c) and digital stereo plotting (stereoscopic evaluation 
of the imagery using lines a/v, a/c or b/c) is possible. These methods 
will play an important role in future photogrammetry. 
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{Text} Introduction 


In 1972 the first earth observation satellite ERTS (LANDSAT-1) was 
launched by NASA. A number of investigators from Japan participated in 
this ERTS program for the satellite data analysis. This was an initial 
step of space remote sensing activities of Japan. NASDA established the 
Earth Observation Center (EOC) at Hatoyama near Tokyo in 1978 and 
started direct reception of LANDS&AT data in 1979. Since then, 

various organizations and universities in Japan have been involved ina 
wide range of remote sensing activities. 


In 1976, NASDA, supported by various groups of Japanese experts, started 
the investigation study of the Japanese needs and eventual contributions 
in earth observation program. Following these activities, in 1978, the 
Space Activities Commission of Japan (SAC), which supervise all space 
activities in Japan, recommended in the “Outline of Japanese Space 
Policy” that “marine and land observation satellite series" should be 
developed in order to establish the earth observation technology and to 
step forward to operational satellites. This recommendation resulted in 
developing two remote sensing satellites, namely Marine Observation 
Satellite-1(MOS-1) and Japanese Earth Resources Satellite-1(ERS-1). 
NASDA started the development of MOS-1 in 1979 and was launched on 19 
February 1987. NASDA also started the R&D of ERS-1 in 1980 and which is 
now in the critical design phase. Development task of ERS-1 is shared 
by MITI and NASDA. MITI is developing the mission equipment (sensors) 
while NASDA the satellite system. ERS-1 is currently planned to be 
launched in 1992. 
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Following the fulfillment of these programs, Japan is to pursue its next 
step earth observation programs. As a next generation earth observation 
Satellite, NASDA plans to launch Advanced Earth Observing Satellite 
(ADEOS) in early 1994. Moreover, a new earth observation program with 
Polar Orbiting Platform(POP) 1s now planned and coordinated among the 
Space Station Partners, namely, NASA, NOAA, ESA, Canada and Japan. 
Currently, NASA, ESA and Japan have their own platform launch plans in 
late 1990’s. The overall schedule of Japanese earth observation 
satellite programs is shown in Table l. 


2 MOS-1 Program 


The Marine Observation Satellite-1(MOS-1), the Japanese first 
experimental earth observation satellite, was launched on 19 February 
1987 successfully by NASDA. The aims of MOS-1 is to establish the basic 
technology for earth observation system, to carry out practical 
observation of the earth (primary the ocean) using on-board sensors and 
to verify the performance of the sensors. 


(1) Satellite System and Mission Instrument 


The profile and description of MOS-1 are shown in Fig. 1 and Table 2. 

MOS-1 has three sensors and a Data Collection System Transponder (DCST). 
These sensors are the Multispectral Electronic Self Scanning Radiometer 
(MESSR), the Visible Thermal Infrared Radiometer (VTIR), and the Micro- 
wave Scanning Radiometer (MSR). The characte~istics of the sensors are 


shown on Table 3. 
(2) MOS-1 Operation 


After the launch of MOS-1 on 19 February 1987, with spacecraft check out 
for three months, images data from the three sensors were flawlessly 
acquired. NASDA is now performing $1ix month experimental operation to 
check the compatibility between spacecraft and ground station for data 
reception and processing. After this test operation, MOS-1 shifts to 
routine operation status. 


The ground support system for MOS-1 is shown in Figure 2. Eart- 
Observation Center (EOC) is the main ground station responsible for 
mission management, scheduling, data acquisition and processing. In 
addition to EOC, MOS-1 data are also planned to be received at nine 
other receiving stations in the world. The world coverage of MOS-1 is 


shown in Figure 3. 
(3) MOS-1 Verification Program (MVP) 


MOS-1 Satellite Program includes the MOS-1 Verification Program (MVP) in 
which the system parameters of mission instrument, function and 
characteristics of sensors and satellite systems, and its data 
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usefulness will be evaluated. In order to obtain fruitful result, NASDA 
publicly announced the opportunity to participate in the MVP. In 
response to this announcement, NASDA received a total of 114 proposals 
from 11 countries. After evaluation 93 proposals were selected. NASDA 
1s now conducting the MVP by using MOS-1 data, airborne data, and 
ground/sea truth data in collaboration with the proposal selected 


organizations. 
(4) MOS-1 Data Distribution 


All the MOS-1 data acquired and archived by any ground station will be 
made available on public and non-discriminatory basis. On 10 August 
1987, NASDA started to distribute MOS-1 data through Remote Sensing 
Technology Center of Japan (RESTEC). The concept of MOS-1 data 
distribution flow is shown in Figure 4. 


3. ERS-1 Program 


Japanese Earth Resources Satellite-1 (ERS-1), with its launching 
scheduled in early 1992, is under development as a joint program of 
Ministry of International Trade and Industry (MITI) and NASDA. The main 
objectives of the ERS-1 are to establish the technology of the active 
microwave sensor, namely the Synthetic Aperture Radar (SAR), and the 
high resolution Optical Sensor (OPS), and to examine the terrestrial 
resources and environment, primarily focusing on the geological and 


topographical survey. 


The profile and description of ERS-1 is shown in Figure S and Table 4. 
ERS-1 is planned to carry two sensors which are the SAR and OPS. The 
characteristics of these sensors are shown in Table S. Their primary 
design was completed in 1987. The critical design has been started. 


4. ADEOS Program 


Following the fulfillment of MOS-1 and ERS-1 programs, NASDA plans to 
launch Advanced Earth Observing Satellite (ADEOS) in early 1994 as a 
next generation earth observation satellite. 


The main objectives of ADEOS are as follows: 
(1) The development of advanced earth observation sensors. 


(2) The development of the modular satellite that is the prototype of 
the future platform. 


(3) The experiments on earth observation data relay using data relay 
satellites to form a global observation network. 


(4) The contribution to domestic and international cooperation by 
installing sensors developed by domestic and/or foreign organizations. 
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The profile and description of ADEOS are shown in Figure 6 and Table 6. 
ADEOS will carry two core sensors which are Ocean Color and Temperature 
Scanner (OCTS) and Advanced Visible and Near Infrared Radiometer 
(AVNIR). In addition to the core sensors, NASDA issued an Announcement 
of Opportunity (AO) to anstall the AO sensors as mentioned above. The 
characteristics of these sensors are shown in Table 7. 





5. Polar Orbiting Platform Program 


Earth observation program by means of Polar Orbiting Platform 
(hereinafter referred to as POP Program) is being established among 
NASA, ESA, STA/NASDA and Canada who are the partners of the Space 
Station Program. The main objectives of POP Program are to guarantee 
the continuity of operational meteorological observation which is 
currently operated by NOAA satellite and to make effective observation 
of the area of earth science. According to the current schedule, NASA 
plans to launch one platform in 1995 and ESA plans to launch one in 
1997. Japan also have a plan to launch its own platform in late 1990’s. 
The countries who participate in POP Program within the frame work of 
the space Station Program will provide on-board sensor instruments. In 
Japan, NASDA’s AMSR (Advanced Micro-wave Scanning Radiometer) and MITI’s 
ITIR (Intermediate Thermal Infrared Radiometer) are planned to be 
provided to NASA POP as Core sensors. The characteristics of AMSR and 


ITIR are shown in Table 8. 


To make this POP program more effective and fruitful, announcement of 
opportunity (AO) for the participation in this program was issued from 
NASA, ESA and STA respectively to the organizations and/or individual 


researchers concerned. 
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Table 1. Schedule of Japanese Earth Observation Satellite Programs 
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Table 2. Description of MOS-1 
Ites Paraseter 
Orbit Hear-Polar (808ka) 
Sua-Syachreseus 
(Lecal Neaa Tise 10:0S48) 
17 days Ceverage Cycle Buraties 
Weight 740k 
Peover Selar Arrey 640(BOL) 
Stabilizaties 3-Anis, Barth-Peistiag 
Lasach Vebicle a- 8 free THSC 
Desiga Life 2 Years 
Iastreneats NESSR (Nultispectral Blectreaic Self 
Scassiag Radieneter) 
VIIR (Visible aed Theraal Iafrared 
Radieseter) 
BSE (Hicrevave Scasniag Radieseter) 
BCS (Data Cellecties Systea) 
Laveach date Peb. 18, 18987 
States Plight segseat PR Coapleted 
PA Coupleted 
aed is flight 
Table 3. Characteristics of MOS-1 Sensors 
Boe Ba 
agsse vTie asa 
Ites 
Ceveleagts (a8) 0.61 - 6.88 6.0- 7.8 
@.61 ° 6.0 6.8 bad 6.7 16 * e 11.8 a 
6.72 ° 6.08 a8 7 12.8 
mY a 1.3 
Prequeeecy (68s) —— 23.8 31.4 
Geosetric Peseleticn oo 8.8 2.7 32 23 
(19069 ite ke) 
Radiesetric Beseletics (Shae) ses" ou is Tt | 
(Ale. 2868) 
Seath eidth (ha) 108 (ene optical Bee ae 
eleseat) « 2 
Seseaiag Bethed electric eek aelen! woabasiocn! 
Ressrts @ Signal te selee ratic esclediang quastisetice seise. 
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Fig. 4 World Coverage Map of MOS-1 
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Fig. S MOS-1 Data Distribution Flow 
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Fig. 6 


Table 4. 





Profile of ERS 


Description of ERS-1 











Recerrest peried 
Leesa] sese tise 
(descending) 
Heaber of orbit per day 
Weight 
Lasach Vebicle 
Besiga Life 
Lavach Sehedele 





1 





ltes Parsseter 
Orbit Sea-syachreosess 
bltitede Appres. $66 te 
Iacliaatios Appres. 97.7° 


dA4daye(Vesteard) 
10:00-11:0048 


Appres. 1§ 
Appres. 1400 tg 
B-1 Becket 

2 Tears 

Barly 1882 
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Table 5S. 






























































Table 6. 








Specification of ERS-1 Sensors 
Sensor 
SAR OPS 
Item 
head | Center Favelength | bead sidth hots 
1 0. 56 0. 08 
2 0. 66 0. 06 
Wavelength (2) L-band 3 0. 81 0. 10 | Bendé is fer off asdia 
4 0. 81 0. 10 | viewiag. 
5 1. 685 O. 11 | Gandd and 4 ashes 
6 2.065 0. 11 | steree-pair 
7 2.19 0.12 
8 2. 4368 0.13 
Pelarizatica H-H 
haw 1 8 (appreainetely) 
Spatiai reselutica 18 (appresiastely) 18 X 24 (appresiaately) 
Off Radia dagle 35 15. 33 (Gead4 ealy) 
moe 
Swath Width (Ka) 78 (appresiastely) 7S (appreniastely) 
Steresscopic Imaging Yes 
Capability 
Description of ADEOS 
ites Pareseter 
° * 
- tapoenas eaatah 
Veight Adppres. 3000 tg 
Peover Aperes. 3.8 be 





Leeseh Vebiele 


Besige Life 
Leesch Sebedele 





8-8 Becket 
(deel Lovech) 


3 Teere 
eT 1883 
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Fig 7 Profile of ADEOS 








Table 7. 


Characteristics of ADEOS Sensors 
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Characteristics of AMSR and ITIR 





CanOiBATE SEASOR 


AUER 
ADVARCED AICROUAVE SCARBING RADIGAETES 


ITIR 
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JAPAN 


Estimation of the Simulated JERS-1 Data 


43070010G Kyoto Selections from INTERNATIONAL ARCHIVES OF PHOTOGRAMMETRY 
AND REMOTE SENSING in English Vol 27 1988 Part 811 pp 770-777 


[Article by T. Takemura, I. Kohno, I. Kakimoto (MSS), Mitsui Mining & 
Smelting Co., Ltd., System Development Office, 2-1-1, Nihonbashi- 
Muromachi, Chuo-ku, Tekyo, 103 Japan; KK. Arai, K. Tonoike (ERSDAC), 
Earth Resources Satellite Data Analysis Center, No 39 Mori Bldg. 2-4-5, 
Azabudai, Minato-ku, Tokyo, 106 Japan) 


[Text] 1. INTRODUCTION 


OPS borne in JERS-1, which will be launched in 1992, is paid attention 
from the world, because it has four bands in short wave infrared region. 
The JERS-1 data were simulated by using the GER (GEOPHYSICAL 
ENVIRONMENTAL RESEARCH CORP.) 64 channels airborne imaging spectrometer 
(AIS) data and evaluated as mineral exploration tool. The area studied 
is Goldfield, Nsvada of U.S.A. 


2. DESCRIPTION OF THE TEST SITE 


Goldfield is located in south central Nevada (Fig. 1) and prospered in 
the gold rush the early years of this century. Gold deposits are 
centered at 3km NE of the town and are distributed in N-S direction. 

The ore bodies are irregular platy veins in the hydrothermally altered 
host rocks, which mainly consist of porphyritic rhyodacite. The altered 
zone is in a circular shape with a diameter of 4km (Fig. 1). Over 130 
tons of gold were mined from this area and a small mining activity 
persists to now. 


3. THE AIRCRAFT DATA 


GER AIS data consist of 63 channels data from visible through infrared 
region (Table 1). Scan parameters are in Table 2. The central 
wavelength of AIS channels and JRES-1 bands are in Table 3. Each pixei 
1s rectangular, 20m in azimuth direction (NS), and 15m in range 
direction (EW). The imagery has been used without geometric correction. 
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The original data 1S provided in digital form and 1s converted to 
radiance by the conversion functions supplied by GER. 


4. JERS-1 DATA SIMULATION 


Energy P i of elecrtromagnetic waves which streamed into the JERS-1 band 
1 detector can be written as follows: 


Pie=- k J Tr(a) © Fi( a) © Hola) da 


where b = constant value iacludiag f auaber of optical systes 
Tr(A)= transeission factor of the condenser systes at 


wavelength A 
Fi(A)e spectral transeission factor of the spectral systes 


of band i at wavelength A 
Ho(A)= radiance gotten by AIS at wavelength A 


Gaussian noise is added to this energy P i to correspond to S/N of band 
1, and P i is digitized into 6 bits. The data thus digitized was used 
as the simulated JERS-1 band i data for study. Flow chart of JERS-1 
data simulation is shown in Figure 2. Main properties of JERS-1 are 


shown in Table 4. 


5. EVALUATION OF THE IMAGES MADE FROM SIMULATED JERS-1 DATA 


The images made from simulated JERS-1 data were evaluated as a mineral 
exploration tool. Simulated images were black-and-white images, 
falsecolor images, ratio images and log residual images. Black-and- 
white images of band 1, 3, 6 and 8 are shown in Photo 1 to 4. These 
images are excellent qualities, and their brightness differences due to 
spectral characteristics are recognized very well. Simulated images 
were evaluated from the standpoint of mapping altered zones. We tried 
to extract and classify the altered zone by studying spectral absorption 
features from simulated OPS data. Spectral reflectance of important 
altered minerals and JERS-1 bands in region of 2.0Oum to 2.Sym are shown 


in Figure 3. (Photos 1 to 4 omitted) 


From these simulated images and ground proof, several characteristics 
were recognized as follows: 


(1) Falsecolor image (band 1,2,3=8,G,R) 


The brightest zones correspond to mine dumps or dried waste pond. 
Goldfield town showed reddish color and it suggests the existance of 


vegetation. 


(2) Falsecolor image (band 5,7,8=R,G,8) 
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Pinkish zones show the altered zone where alterection clay minerals are 
alunite, kaolinite, montmorillonite and sericite. These minerals have 
typical spectrum absorption near 2.2m, so the brightness of band 7 is 
relatively darkened and the altered zone became pink. However, the 
color difference among these minerals could not be recognized. 


(3) Ratio image (bandS/band7, band6/band7, band7/band8 = 8,R,G) 


Violet zones correspond to the altered zone as the pinkish colored zones 
of the preceding image. 


(4) Log Residual image (band 1,2,3 = 8,G,R) 


The log residual technique is the method to remove effects of 
atmospheric scattering and absorption, and brightness difference due to 
Slope orientation. This technique is effective to be applied in the 
area with many kinds of clay minerals and sparsely existing vegetation 
like Goldfield area because the minute change in reflectance is 


emphasized. 


Image (4) has more variety of colors than image (1). Greenish zone 
corresponds to the iron oxide zone. This greenish zone coincides with 
the area extracted as the hematic zone by Chebyshev waveform analysis of 
GER that used 63 channels data. This fact shows the validity of JERS-1 


data. 
(5) Log Residual image (band 5,7,8 = R,G,8) 


In image (5), the altered zone is in dark reddish, reddish and pinkish 
zones. The change in color suggests the differences of minerals. 


6. CONCLUSION 


The results of this JERS-1 data simulation over Goldfield area are as 
follows: 


(1) The area covered with vegetation can be extracted by using 
falsecolor images, ratio images and log residual images. 


(2) The altered zone where alterection clay minerals are 
montmorillonite, alunite, sericite, kaolinite, etc., can be extracted by 
using falsecolor images and ratio images. 


(3) The iron oxide zone can be extracted by using log residual images. 


(4) Clay minerals in the altered zone could be classified in log 
residual images. 


(S) We could not evaluate the possibility for extraction of carbonate 
rocks like limestone that have spectrum absorption in JERS-1 band 8 
because of no carbonate rock in the study area. 
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Map of the Test Site 
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able 4? 1 
— 
RESOLUTION 18.38 xX 24.28 
SWATH WIDTH 75ke mn 
| BAND . NGTH BAND WIDTH S/N 
BAND | 0.56( 4 8) 80 (ne) 46 (db) 
BAND 2 0.66( ue) 60 (ne) 46 (db) 
BAND 3 0.81( 48) 100 (ne) 46 (db) 
BAND 4(STEREOSCOPE) 0.81( ue) 100(n8) 46 (db) 
BAND § 1.65( 4 @) 110(ne) 42 (db) 
BAND 6 2.06( ue) 110(as) $2 (db) 
- BAND 7 2.19( 4 8) 120(ns) $3 (db) 
BAND 8 .3 $0 7(db 
BASE HEIGHT RATIO IN 0.3 
| STEREOSCOPIC IMAGEING | 
UANTI ZATION 6 BITS 
bps x 208. 














Table 3. 


CaGNVd 


CHANNEL WAVELENGTH 


'tanvad 








tEanvgd 





— fee pee pee fee fee pee 
eon Wwe © oO ee -3 BO cr & CrP — 


The Central Wavelength of AIS Channels and JERS-1 Bands 


a — i — oo — 2 — 2 — 2 — 2 — 2 — 2 — 2 — 2 — 2 — 2 2 — 2 — 2 2 2 2 2 — 2 — 2 — 2 2) 


(um) 


. 4336 
4570 
4804 
. 5038 
5272 
. 5506 
9740 
5974 
6208 
6442 
. 6676 
6910 
7144 
. 1378 
7612 
. 1846 
. 8080 
8314 
8548 
8782 
9016 
. 9250 
9484 
9718 
0800 


SGNV@ 


LaNnNvV@ 


CHANNEL WAVELENGTH 








S9aGNV@ 
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(um) 


1.2000 
1.3200 
1.4400 
1.5600 
1.6800 
1.8000 
1.9860 
2.0024 
2.0189 
2.0358 
2.0517 
2.0682 
2.0846 
2.1010 
2.1174 
2.1389 
2.1508 
2.1667 
2.1882 
2.1996 
2.2160 
2.2825 
2.2489 
2.2658 
2.2817 


CHANNEL WAVELENGTH 





= 


SGNVd 


(um) 


2.2982 
2.3146 
2.3310 
2.3475 
2.3639 
2.3808 
2.3968 
2.4182 
2.4296 
2.4460 
2.4625 
2.4789 
2.4953 
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SIMULATED JERS-1 DATA | 


Fig. 2 Flow Chart of JERS-1 Data Simulation 
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Fig. 3 Spectral Reflection Factor of Altered Minerals 
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